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ABSTRACT 


Two  recent  computer  programs,  WRCOMP  and  ENGINE,  by  A. 
Mathur,  were  used  to  examine  the  performance  to  be  gained  by 
incorporating  a  wave  rotor  component  in  a  turbofan  engine  with 
mixed  exhausts.  The  programs  were  transferred  to  a  VAX-?oon 
computer,  extended,  and  test  cases  reported  by  A.  Mathur  were 
successfully  reproduced.  A  comparison  was  made  between 
ENGINE,  in  which  real  gas  effects  are  accounted  for,  and  ONX 
(by  J.  Mattingly)  in  which  constant  specific  heats  are  used. 
The  inclusion  of  real  gas  effects  proved  to  have  a  significant 
impact  on  the  predicted  performance.  An  extension  of  Mathur ' s 
results  was  made  by  varying  the  overall  pressure  ratio  in  the 
t  wave-turbofan  engine.  Further  cycle  studies  and  experiments 

to  measure  wave  rotor  component  performance  are 
recommended . 
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I. 


INTRODUCTION 


Interest  by  DARPA  (Defense  Advance  Research  Projects 
Agency)  and  NAVAIRSYSCOM  (Naval  Air  Systems  Command)  in 
obtaining  increased  range  for  unmanned  air  vehicles  (UAV's) 
prompted  the  present  study  of  wave  rotor  applications  in 
engines.  The  wave  rotor  offers  potential  to  achieve  lower 
specific  fuel  consumption  (SFC)  and  higher  specific  thrust  by 
allowing  increased  cycle  temperatures  without  increasing 
turbine  inlet  tem.perature .  The  goal  of  this  study  was  to 
examine  the  potential  benefits  of  incorporating  a  wave  rotor 
in  a  turbofan  engine  with  mixed  exhausts.  The  results  showed 
that  significant  benefits  could  be  obtained  in  both  SFC  and 
specific  thrust  depending  upon  the  selection  of  other  cycle 
parameters . 

Chapter  II  describes  the  wave  rotor  concept,  presents  a 
general  background  on  trends  in  gas  turbine  technology,  gives 
a  description  of  a  wave  rotor  engine  component  and  introduces 
the  two  computer  programs  (ENGINE  and  ONX)  which  were  used  in 
the  study.  Chapter  III  gives  a  detailed  description  of  the 
ENGINE  code  and  presents  a  comparison  of  results  obtained  with 
ENGINE  and  ONX.  Chapter  IV  gives  a  description  of  the  WRCOMP 
program,  which  analyses  the  unsteady  wave  rotor  flow,  and  its 
interface  with  the  ENGINE  program.  Chapter  V  gives  results 
obtained  using  the  codes  for  the  performance  of  a  turbofan 
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with  mixed  exhausts,  with  and  without  a  wave  rotor  component. 
Conclusions  and  recommendations  for  further  work  are  given  in 
Chapter  VI.  Details  are  given  in  the  appendices.  Appendices 
A  and  B  contain  procedural  guides,  sample  results  and  program 
listings  for  ENGINE  and  WRCOMP,  respectively,  while  Appendix 
C  describes  file  transfer  and  graphics  procedures  which  were 
used  in  the  course  of  the  study. 
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II.  BACKGROUND 


A.  POTENTIAL  FOR  WAVE  ROTOR  APPLICATIONS  IN  ENGINES 

A  wave  rotor  is  a  partial  admission  device  which  causes 
one  gas  to  compress  another  by  wave  propagation  [Ref.  l:p.  1] . 
A  series  of  straight  compression  tubes  are  aligned  on  a  rotor 
encased  in  a  drum.  The  rotor  rotates  on  a  fixed  axis.  Flow 
in  and  out  of  the  wave  rotor  is  steady  while  flow  inside  the 
wave  rotor  tubes  is  unsteady.  Figure  1  shows  an  example  of  a 
wave  rot^r  component.  The  single  wave  rotor  component  serves 
the  same  function  as  the  combination  of  a  compressor  and  a 
turbine  mounted  on  the  same  shaft. 


Figure  1.  Wave  Rotor  Component  [Ref.  l:p.  294] 

Figure  2  shows  a  turbofan  engine  configured  with  an 
axially  mounted  wave  rotor  component. 
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Figure  2.  Wave  Rotor  Turbofan  Engine  [Ref.  l:p.  293] 


Development  of  wave  rotor  devices  for  engines  began  in  the 
early  1950 's.  By  the  early  1960 's,  gains  were  achieved  in 
performance  through  improved  seal  design  and  improved  timing 
of  the  waves  arriving  at  inlet  and  outlet  ports  ("tuning") 
[Ref.  l:p  174].  Few  developments  were  attempted  after  the 
1960 's,  except  by  Brown-Boveri  in  Switzerland.  The  Brown- 
Boveri  "Comprex"  is  presently  being  produced  as  a  commercial 
supercharger  and  is  available  in  trucks,  heavy  equipment  and 
Mazda  diesel  automobiles.  Reference  1  provides  a  comprehen¬ 
sive  review  of  wave  rotor  technology. 

Wave  rotor  components  have  promising  applications  in  small 
engines  in  the  600  to  1000  lb.  thrust  range.  Such  engines 
have  applications  in  cruise  missiles,  remotely  piloted 
vehicles  (RPV*s),  helicopters  and  small  thrust  turbofans.  The 
more  conventional  approach,  however,  has  been  to  pursue 
increased  turbine  inlet  temperatures  in  conventional  gas 


turbine  engines  to  achieve  the  desired  improved  performance. 
However,  in  the  process,  the  requirement  for  turbine  cooling 
also  increased.  The  wave  rotor  offers  the  potential  to 
increase  performance  without  the  attendant  necessity  for 
cooling. 

Studies  of  engine  performance  indicate  that  a  wave  rotor 
component  must  be  included  in  the  initial  cycle  optimization 
rather  than  be  added  to  an  existing  engine  design.  Wave  rotor 
components  require,  optimally,  different  compressor  and  fan 
pressure  ratios.  The  inclusion  of  a  wave  rotor  reduces  the 
number  of  stages  required  for  compression  to  a  constant  cycle 
pressure  ratio  and  allows  higher  maximum  temperatures  in  the 
cycle  for  the  same  turbine  inlet  temperature.  Performance 
predictions  for  gas  turbine  engines  incorporating  wave  rotor 
components  show  decreased  specific  fuel  consumption  and 
increased  specific  thrust  for  any  fixed  value  of  turbine  inlet 
temperatures . 

Figure  3  shows  the  gains  that  are  thought  to  be  achievable 
in  the  performance  of  small  engines  [Ref.  l:p.  292].  In  Figure 
3,  "near  term  improvements"  include  higher  rpm  (for  a  given 
engine  size)  and  higher  turbine  temperatures.  "New  engines" 
will  include  composites,  ceramics,  or  new  materials.  Some 
current  "revolutionary  engines"  include  the  compound  cycle, 
the  eccentric,  the  recuperated  and  the  wave  rotor  engine  [Ref. 
l:p.  292].  The  conventional,  eccentric  and  the  recuperated 
engines  are  fundamentally  limited  by  the  maximum  allowable 
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Figure  3.  Trends  in  Engine  Development  [Ref.  l:p.  292] 

turbine  inlet  temperature  (TIT)  [Ref.  2;p.  60).  The  wave 

rotor  engine  is  the  only  engine  that  promises  significantly 
increased  performance  without  requiring  an  increase  in  TIT 
[Ref.  l:p.  293]. 

B.  ENGINE  PERFORMANCE  CODES 

Two  engine  performance  codes  were  available  for  the 
present  study,  ENGINE  and  ONX.  There  were  some  differences 
between  the  two  programs.  ENGINE  calculated  real  gas  effects 


based  on 

fuel-air 

ratio 

and 

temperature,  while  ONX 

used 

constant 

specific 

heats 

at 

the  compressor  outlet 

and 

downstream  of  the  turbine  inlet  station.  In  both  programs, 
the  engine  performance  is  obtained  by  equating  power  balances 
across  the  turbine  and  compressor  sections,  and  describing  an 
expression  for  the  enthalpy  rise  a-  ’-oss  the  burner  section. 

Program  ENGINE  was  developed  by  A.  Mathur  [Refs.  3;  4], 
and  uses  notation  developed  by  Vavra  [Ref.  5],  and  Vanco 
[Ref.  6].  The  program  calculates  the  performance  of  a 
turbofan  engine  with  mixed  exhausts  with  and  without  the  wave 
rotor  component. 

Program  ONX  was  developed  by  J.D.  Mattingly  [Ref.  7],  and 
it  uses  methods  and  notation  developed  by  G.  C.  Oates  [Ref. 
8].  The  program  calculates  the  performance  for  a  variety  of 
engine  configurations  including  a  turbofan  with  mixed 
exhausts . 

C.  WAVE  ROTOR  INTERNAL  FLOW  DESCRIPTION 

In  1984,  A.  Mathur  developed  a  FORTRAN  program  (WRCOMP)  to 
analyze  the  unsteady  processes  in  wave  rotors  [Ref.  9].  A 
detailed  description  of  the  analysis  and  code  are  provided  in 
Reference  9.  WRCOMP  uses  a  1-D,  random  choice  method  (RCM)  to 
solve  the  hyperbolic,  nonlinear  conservative  system  of 
equations  by  solving  localized  Riemann  problems.  Shocks  have 
sharp  resolution.  Discontinuities  are  modelled  as  jumps.  The 
code  is  first  order  accurate  in  time  and  it  is  unconditionally 
stable.  Further  explanations  and  details  on  the  Random  Choice 
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Method  are  given  by  Glimm  [Ref.  10],  Chorin  [Ref.  11],  and  Sod 
[Ref.  12]. 

WRCOMP  calculates  the  unsteady  process  inside  the  wave 
rotor,  outputs  the  inlet  and  outlet  port  opening  and  closing 
times  required  for  correct  design,  and  computes  the  perfor¬ 
mance  of  the  wave  rotor  given  the  total  pressure  ratio  (TPR) , 
the  static  pressure  ratio  (SPR) ,  the  exit  Mach  number  (AMEXIT) 
and  the  (initial)  reference  density  inside  the  wave  rotor. 
The  performance  calculations  depend  on  a  favorable  comparison 
for  the  mass  flow  rate  into  and  out  of  the  wave  rotor 
component.  The  outputs  from  WRCOMP  are  used  as  inputs  to 
engine  performance  calculations  performed  using  the  ENGINE 
code . 
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III.  ENGINE  PERFORMANCE  CODE 


A.  DESCRIPTION  OF  THE  CODE 

A  complete  listing,  procedural  guide  and  sample  results 
for  ENGINE  are  provided  in  Appendix  A.  The  code  incorporates 
modular  subroutines  to  calculate  each  section  of  the  engine 
cycle.  A  complete  derivation  of  the  code  can  be  found  in 
Reference  4. 

The  following  changes  were  made  to  improve  the  utility  of 
the  code: 

(1)  The  original  data  statement  format  was  modified  for  a 
user  input  interface; 

(2)  A  pop-up  menu  depicting  the  configurations  was  added 
(see  Appendix  C  for  details) ; 

(3)  An  altitude  table  was  incorporated  to  automate  ambient 
conditions  using  Reference  13; 

(4)  A  graphics  plotting  routine  was  added  using  GRAFkit 
software,  Reference  14; 

(5)  A  loop  was  constructed  to  increment  for: 

a)  bypass  ratio, 

b)  fan  pressure  ratio, 

c)  compressor  pressure  ratio, 

d)  LPT  exit  Mach. 

A  complete  listing  of  the  input  parameters  and  a  sample  output 
result  appears  in  Appendix  A. 

ENGINE  (and  ONX)  computes  engine  performance  in  terms  of 
uninstalled  specific  thrust  and  specific  fuel  consumption. 
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The  uninstalled  specific  thrust  is  defined  as  the  thrust  per 
unit  mass  flux  through  the  engine  and  determines  the  engine 
size.  The  specific  fuel  consumption  is  the  mass  flow  rate  of 
fuel  per  unit  thrust  and  determines  the  range.  A  power 
balance  is  written  for  each  turbine  and  compressor  spool.  An 
enthalpy  rise  is  written  for  the  flow  through  the  combustor. 
ENGINE  uses  loss  coefficients  and  efficiencies  to  account  for 
component  performance.  The  program  also  uses  real  gas  effects 
based  on  fuel  air  ratios  and  temperature.  In  contrast,  ONX 
uses  pressure  ratios  and  polytropic  efficiencies  to  account 
for  component  performance  and  assumes  constant  specific  heats 
across  the  compressor  (Cp^)  and  downstream  of  the  combustor 
(Cpt)  .  A  necessary  step,  prior  to  the  use  of  ENGINE,  or 
comparison  between  ENGINE  and  ONX,  was  to  reproduce  the 
results  reported  by  A.  Mathur  in  Reference  4.  These  results 
were  successfully  duplicated  using  the  inputs  given  in 
Reference  4  and  Appendix  A. 

To  achieve  an  accurate  comparison  between  ENGINE  and  ONX, 
the  polytropic  efficiencies  req[uired  in  ONX  were  computed  from 
the  outputs  of  the  ENGINE  code  using  the  component  efficien¬ 
cies,  the  actual  and  isentropic  pressures  and  temperatures, 
and  gamma.  A  sample  calculation  is  provided  in  Appendix  A. 

B.  COMPARISON  WITH  ONX  CODE 

Since  engine  cycle  calculations  are  greatly  affected  by 
the  inputs,  especially  the  component  efficiencies,  an  accurate 
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comparison  between  ENGINE  and  ONX  could  only  be  performed  if 
the  flight  parameters,  the  design  choices,  the  design 
limitations  and  the  component  efficiencies  were  the  same. 
Flight  parameters  include  the  flight  Mach  number  and  the 
altitude.  Design  choices  include  the  maximum  temperature  in 
the  burner,  the  lower  fuel  heating  value,  the  maximum  turbine 
inlet  temperature,  the  low  pressure  turbine  exit  Mach  number, 
the  compressor  pressure  ratio,  the  fan  pressure  ratio,  the 
bypass  ratio,  etc.  Fans,  compressors,  burners,  turbines  and 
nozzles  are  characterized  by  component  efficiencies,  which  are 
the  zc...plimei.t£.  to  the  loss  coefficients.  The  simplest  test 
case  for  comparison  of  ENGINE  and  ONX  was  the  non-afterburning 
turbojet.  ENGINE  calculated  the  turbojet  case  by  setting  the 
bypass  ratio  equal  to  zero. 

1.  TURBOJET 

With  all  input  and  component  performance  parameters 
matched,  the  results  of  ENGINE  were  compared  to  those  of  ONX 
for  the  case  of  the  turbojet.  The  results  are  shown  in  Table 
I.  For  the  ONX  (1)  case,  the  Cp^  and  Cpt  were  selected  to  be 
in  the  range  that  Mattingly  referred  to  as  "typical”  values 
[Ref.  7:pp.  116,  438].  The  ONX  (2)  values  for  Cp^  and  Cpt  were 
the  result  of  numerous  iterations  using  ONX  to  try  to  achieve 
agreement  with  ENGINE.  The  ONX  (3)  values  for  Cp^  and  Cpt  were 
obtained  by  computing  the  average  values  for  Cp^  and  Cp^  from 
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TABLE  I 


TURBOJET  PERFORMANCE  COMPARISON 
Mo  =  0.79,  h=35000  ft. 


ENGINE:  F/m  =  41.88  Ibf/ (Ibm/sec) 

SFC  =  0.899  (Ibm  fuel/hr) /Ibf  thrust 


ONX 

ONX 

ONX 

(1) 

(2) 

(3) 

Typical 

Specific 

heats 

Values  req'd 
for 

agreement 

Average 
from  ENGINE 
results 

CPc 

0.235 

0.235 

0.2434 

CPt 

0.295 

0.251 

0.2508 

F/m 

55.18 

40.04 

35.94 

SFC 

0.992 

0.954 

1.001 

% 

Difference 

F/m 

31.76 

4.38 

14.18 

% 

Difference 

SFC 

10.34 

6.10 

11.34 

the  ENGINE  program  output  data.  The  differences  were  seen  to 
be  surprisingly  large 

A  comparison  table  was  constructed  of  the  pressure  and 
temperature  into  each  component  to  trace  where  the 
calculations  began  to  diverge.  The  divergence  was  found  to 
occur  in  the  enthalpy  rise  calculation  across  the  combustor. 
For  ONX,  the  enthalpy  rise  gives  the  fuel-air  ratio  (f)  as 
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(1) 


[Cp^yT^f+h^xT\  -Cp^xT^^ 


and  the  test  case  in  Table  I  gave  f  =  0.0152.  For  ENGINE,  the 
fuel-air  ratio  equation  was: 

[Cp fXT h^xT\  -B{  Tj.^) 


where , 


and 


c  c  c  c 

A{T)  ^r  xT-  —  xT^  +  —  xT^-—^xT^  +  —^xT^ 
1  2  3  4  5 


,  >  ,  P,  ,  D.  .  D.  . 

BiT)  =r>  xT+—^xT^ - lxT^  +  -±xT^ - ^xT^ 

2  3  4  5 


(3) 


(4) 


With  fi  =  0.0  (f  prior  to  the  burner), the  test  case  designated 
ONX  (1)  in  le  I  gave  f  =  0.0103.  Thus  ONX  calculated  a 
32.24%  higher  value  for  the  fuel  air  ratio  than  did  ENGINE. 
The  question  then  was,  can  agreement  be  obtained  between  ONX 
and  ENGINE  using  any  reasonable  values  of  the  Cp's  in  ONX. 
The  result  of  an  attempt  to  achieve  such  agreement  is  shown  as 
ONX  (2)  in  Table  I.  Using  Cp^  =  0.235  and  Cp^  =  0.251,  the 
disagreement  in  specific  thrust  was  reduced  to  4.4%  and  in  SFC 
to  6.1%.  It  is  noted  here  that  the  Cp  value  for  air  at  the 
temperature  at  the  exit  of  the  combustor  (1860°R)  in  the 
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present  example  is  Cp  =  0.271,  considerably  higher  than  the 
value  required  to  approach  agreement. 

One  further  attempt  was  made  to  compare  the  code 
predictions  for  a  turbojet.  Values  of  Cp  were  calculated  at 
the  inlet  and  outlet  stations  of  the  compressors  and  turbines 
using  the  enthalpies  calculated  at  the  stations  by  ENGINE  in 
the  expression,  h  =  CpT.  An  average  of  the  in  and  out  values 
was  then  used  as  an  input  to  ONX.  The  results  are  shown  as 
ONX  (3)  in  Table  I.  The  strong  sensitivity  of  the  performance 
to  the  values  used  for  the  specific  heats  is  evident  in  the 
comparison  of  the  results  of  ONX  (3)  with  ONX  (2). 

The  selection  of  Cp  values,  therefore,  had  a 
significant  affect  on  the  performance  calculations  for 
specific  thrust  and  specific  fuel  consumption.  In  order  to 
examine  the  sensitivity  to  specific  heats,  ONX  was  run  with 
varying  Cp^  and  varying  Cp^.  Figures  4  and  5  show  the  affect 
that  Cpt  and  Cp^  had  on  the  specific  thrust  and  SFC  given  by 
the  ONX  code.  The  sensitivity  to  the  selection  of  the  input 
values  for  ONX  is  clearly  evident  in  these  figures. 

A  complete  deviation  of  the  real  gas  equations 
programmed  in  ENGINE  can  be  found  in  Reference  3 .  Further 
details  can  be  found  in  References  5  and  6.  In  short,  the 
stoichiometric  combustion  equation  for  a  general  hydrocarbon 
fuel  (CH2)n  is  expressed  as  a  polynomial  expansion.  The 
constants  in  the  expansion  to  C5  and  Di  to  D5  in  Equations 
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Figure  4.  Effect  of  Cpt^  on  Specific  Thrust  and  SFC 
(  Yc=  Yt=  1-3  50,  Cp^  =  0.2434) 

(3)  and  (4)  are  deiived  from  the  constituents  of  the  combus¬ 
tion  gases.  The  final  form  of  the  enthalpy  rise  for  the 
combustor  is  given  as  Equation  (2).  Thus  Cp  is  effectively, 
and  properly,  a  function  of  fuel-air  ratio  and  temperature. 
Since  the  ENGINE  code  calculated  the  enthalpy  using  the  fuel- 
air  ratio  and  temperature  rather  than  maintaining  constant 
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Figure  5.  Effect  of  Cp^  on  Specific  Thrust  and  SFC 
{  =  1.400,  Yt  =  1.350,  Cp^  =  0.2508) 

values  for  specific  heats,  the  ENGINE  code  results  were  now 
considered  to  be  more  reliable  than  the  ONX  results. 

2 .  TURBOFAN 

After  successfully  reproducing  A.  Mathur's  results 
[Ref.  4:p.7]  for  a  turbofan  with  mixed  exhausts,  a  comparison 
between  ENGINE  and  ONX  was  carried  out  for  the  same  engine. 


The  average  values  of  Cp's  given  by  the  ENGINE  results  were 
input  into  ONX.  The  results  are  presented  in  Table  II.  The 
percent  differences  in  the  specific  thrust  and  SFC,  predicted 
by  ONX  in  comparison  to  ENGINE  were  11%  aid  23%  respectively. 
[Ref.  4:p.  7] 

TABLE  II 

TURBOFAN  PERFORMANCE  COMPARISON 


Mo=  0.79,  h=35000  ft. 

ENGINE  [Ref.  4:P.7]: 

F/m  =  30.185  Ibf/ (Ibm/sec) 

SFC  =  0.832  (Ibm  fuel/hr) /Ibf  thrust 


Difference 

F/m 

% 

Difference 

SFC 


ENGINE 

ONX 

Cp,=0.2434 

Cpt=0. 2508 

30.18 

33.837 

0.832 

1.0926 

0 

10.8 

0 

23.8 

The  bypass  ratio  was  varied  keeping  other  parameters 
fixed.  Using  ENGINE,  a  very  narrow  range  of  bypass  ratio  was 
found  for  which  the  static  pressures  of  the  bypass  stream  and 
the  core  flow  could  be  equal  at  the  constant  pressure  splitter 
plate  aft  of  the  low  pressure  turbine.  The  results  in  Figu.'e 
6  show  that  common  trends  were  obtained  with  ENGINE  and  ONX. 
As  bypass  ratio  was  increased,  both  SFC  and  specific 
thrust  decreased.  However,  the  effect  on  specific  thrust  was 
the  largest.  In  fact,  only  a  5%  decrease  in  SFC  could  be 
obtained  at  the  cost  of  a  20-30%  decrease  in  specific  thrust. 


Figure  6.  Engine  Performance  Trends  for  ONX  and 
ENGINE  Varying  Bypass  Ratio 

(Mo  =  0.79,  h  =  35000  ft,  Ttj  =  2.3,  =  10.87,  TIT  =  1860  °R) 
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To  complete  the  comparison  of  ENGINE  and  ONX  and  to 
examine  the  predictions  of  the  ENGINE  program,  the  following 
parameters  were  varied  and  the  results  obtained  for  specific 
thrust  and  SFC  were  plotted  in  the  indicated  figures:  1)  exit 
Mach  number  from  the  LPT  (Figure  7)  ;  2)  maximum  temperature  in 
the  combustor  (Figure  8) ;  3)  fan  pressure  ratio  (Figure  9) ; 
and  4)  compressor  pressure  ratio  (Figure  10) . 


e  .1  .2  .3  .4  .5  .b  .7  .e  .<»  i.« 

4M4CH8 


Figure  7.  Effect  of  LPT  Exit  Mach  Number  (AMACH8) 
on  Specific  Thrust  and  SFC 

(M„  =  0.79,  h  =  35,000  ft,  tIj  =  2.3,  n,  =  10.87,  TIT  =  1860°R) 
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As  shown  in  Figure  7,  the  LPT  exit  Mach  number  had 
little  effect  on  the  specific  thrust  or  on  the  SFC. 

In  Figure  8,  results  are  shown  for  varying  turbine 
inlet  temperature  obtained  using  ENGINE  with  bypass  ratio 
allowed  to  vary  and  with  ONX  with  bypass  ratio  fixed.  In  the 
cases  shown  in  Figure  8,  the  bypass  ratio  of  ENGINE  was 
allowed  to  vary  such  that  the  static  pressure  at  the  splitter 
plate,  where  the  bypass  air  is  mixed  with  core  air,  was 
matched.  The  mixed  exhaust  forced  higher  values  of  bypass 
ratio  as  the  TIT  was  increased.  The  higher  bypass  ratio 
allowed  the  specific  thrust  to  remain  fairly  constant.  Since 
significantly  more  air  was  bypassed,  less  fuel  was  required  to 
produce  the  same  specific  thrust  as  the  TIT  was  increased. 
Since  the  bypass  ratio  increased,  the  SFC  decreased  as  the  TIT 
increased. 

The  results  from  ONX  plotted  on  Figure  8  show  the 
effect  of  increasing  TIT  while  keeping  bypass  ratio  constant. 
Fixing  the  bypass  ratio  resulted  in  increased  specific  thrust 
but  at  the  price  of  increased  SFC. 
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Figure  8.  Effect  of  Varying  TIT  on  Specific  Thrust  and  SFC 
(Mo  =  0.79,  h  =  35,000  ft,  tif  =  2.3,  71^  =  10.87,  TIT  =  ISeO^R) 

The  range  of  fan  pressure  ratio  was  taken  to  be 
between  1.0  and  4.0.  Figure  9  shows  the  performance 
calculated  as  the  fan  pressure  ratio  was  varied  between  1.0 
and  3.4.  ENGINE  required  that  the  bypass  ratio  vary  to 
satisfy  the  equal-static-pressure  condition  at  the  splitter 
plate.  A  very  narrow  band  of  fan  pressure  ratios  was  obtained 
if  the  bypass  ratio  was  also  specified.  ONX  produced  a 
broader  range  of  values  for  a  specified  bypass  ratio. 
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PRLC 


Figure  9.  Effect  of  Varying  PRLC  on  Specific  Thrust  and  SFC 
(Mo  =  0.79,  h  ==  35,000  ft,  *  10.87,  TIT  =  1860°R) 


The  result  of  varying  compressor  pressure  ratio  is 
shown  in  Figure  10.  The  bypass  ratio  was  fixed  at  P  =  0.5. 
ENGINE  predicted  a  minimum  SFC  for  compressor  pressure  ratio 
range  between  11.5  and  12.5  (whereas  the  minimum  predicted  by 
ONX  was  greater  that  20)  .  It  is  clear  from  the  curves  in 
Figure  10  that  a  pressure  ratio  less  than  that  for  optimum  SFC 
would  be  chosen  to  obtain  a  significant  improvement  in 
specific  thrust. 
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Figure  10.  Effect  of  Varying  PRHC  on  Specific  Thrust  and  SFC 


(Mo  =  0.79,  h  =  35,000  ft,  Ti^  =  2.3,  TIT  =  1860'’R) 


IV.  WAVE  ROTOR  ANALYSIS 


A.  SUMMARY  DESCRIPTION 

The  wave  rotor  component  and  the  combustion  chamber  can  be 
considered  to  be  combined  to  form  a  two  port  component  called 
a  gas  generator.  Other  models  of  the  wave  rotor  component  are 
possible,  but  the  gas  generator  model  is  the  most  straight¬ 
forward.  The  equations  and  code  for  the  wave  rotor  flow 
modeling  used  in  the  WRCOMP  program  are  described  in 
Reference  9.  The  equations  for  the  wave  rotor  component  as  a 
gas  generator  included  in  the  ENGINE  code  are  derived  in 
Reference  4.  Figure  11  shows  an  unwrapped  wave  rotor  and  a 
simplified  wave  diagram  [Ref.  l:p.  38],  Reference  1  describes 
the  gas  generator  process  as  a  filling,  emptying,  filling  and 
emptying  process. 

The  process  can  equally  be  described  as  two  scavenging 
processes  separated  by  periods  of  stationary  flow  within  the 
rotor.  There  is  a  low  pressure  scavenge  at  station  (4)  and  a 
high  pressure  scavenge  at  station  (2).  The  combustion  exhaust 
gas  initially  compresses  the  air  inside  the  wave  rotor  and 
then  is  scavenged  to  (4)  by  the  incoming  air  at  the  low 
pressure  port  (1)  .  When  the  high  pressure  exhaust  port  opens, 
the  high  pressure  gas  inside  the  wave  rotor  exits  to  the  HPT 
and  effectivaly  scavenges  the  wave  rotor  [Ref.  15:p.  63]. 
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(USUALLY  COMBUSTOR  FOR  GAS  GEN.  OR 


Figure  11.  Unwrapped  Wave  Rotor  Component  [Ref.  l:p.  38] 

In  the  ENGINE  program,  the  gas  generator  is  modelled  as  an 
internal  combustion  process  in  which  a  gain  in  total  pressure 
and  a  rise  in  the  total  temperature  is  achieved.  The 
performance  of  the  gas  generator  is  expressed  as  a  function  of 
the  static  pressure  ratio,  the  total  pressure  ratio  and  the 
total  temperature  ratio  across  the  inlet  ports  [Ref.  3;p.  14]. 

B.  INTERFACE  WITH  ENGINE  PERFORMANCE  CODE 

WRCOMP  can  be  run  alone  or  in  conjunction  with  the  ENGINE 
code.  Since  the  VAX  system  operates  in  windows,  both  programs 
can  be  run  simultaneously.  Both  procedures  for  running  WRCOMP 
are  outlined  in  Appendix  B.  The  original  code  used  data 
statements  in  SI  units.  It  was  modified  to  accept  input 
pressures  in  psi,  temperatures  in  °R,  density  in  Ibm  per  ft^. 
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gamma,  fuel  air  ratio  (f) ,  and  wave  rotor  exit  Mach  number. 
In  the  coding,  the  analysis  is  nondimensionalized  to  the  input 
parameters.  An  output  file,  WRPERFl.OUT,  was  incorporated  to 
store  the  input  parameters  and  to  convert  the  English  units  to 
SI  units.  Routines  were  added  to  compute  the  procedures 
outlined  in  Reference  4  and  to  send  the  output  to  WRPERFl.OUT. 
The  DISSPLA  plotting  routines  were  disabled.  In  practice, 
WRCOMP  is  run  until  the  mass  flow  into  the  wave  rotor  is  equal 
to  the  mass  flow  out.  Adjustments  are  made  progressively  to 
RREF,  the  reference  density  inside  the  wave  rotor,  until  the 
mass  flows  are  equal.  For  1681  iterations,  a  CFL  number  = 
0.6,  a  grid  cell  width  of  0.01,  and  performing  2  0  Riemann 
iterations,  the  total  run  time  was  about  14  minutes  on  the 
VAX.  An  example  output  is  presented  in  Appendix  B. 

From  experience  with  WRCOMP  and  ENGINE,  it  is  easier  to 
keep  the  overall  pressure  ratio  constant  for  a  series  of  test 
cases.  By  keeping  the  overall  pressure  ratio  constant,  the 
outlet  pressure  from  the  HPC  stays  the  same.  If  the  outlet 
pressure  from  the  HPC  is  the  same  and  the  performance 
parameters  for  the  wave  rotor  are  kept  the  same,  then  the 
value  for  the  reference  density  which  satisfies  the  mass  flow 
solution  for  WRCOMP  also  stays  the  same. 

In  practice,  it  is  better  to  vary  the  overall  pressure 
ratio  and  maintain  the  same  turbine  inlet  temperature.  This 
allows  optimizing  the  fan,  the  compressor  and  the  wave  rotor 
based  on  realistic  cycle  constraints. 
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V.  LOW  BYPASS  TURBOFAN  CYCLE  CALCULATIONS 


A.  GENERAL  CONSIDERATIONS 

The  present  work  has  been  concentrated  on  the  low  bypass 
turbofan  engine,  and  what  performance  improvements  can  be 
obtained  by  incorporating  a  wave  rotor  component.  The  initial 
calculations  have  been  made  without  reference  to  a  specific 
vehicle  and  mission.  Considerations  which  underlie  the 
calculations  and  which  limit  their  scope  are  discussed  in  the 
following  paragraphs - 

The  characteristics  of  low  bypass  turbofans  provide  a 
reasonable  match  to  the  performance  requirements  for  cruise 
missiles  and  RPV's.  However,  the  design  limitations  for  the 
components  of  a  cruise  missile  may  vary  significantly  from 
those  of  an  RPV,  even  though  the  two  engines  may  have  very 
similar  thrust  recfuirements .  A  cruise  missile  may  require 
compact  packaging  and  "low  cost  design"  as  it  is  considered  a 
"throw-away"  or  one-time-use  engine,  thereby  relaxing  the  TIT 
and  cooling  requirements.  An  RPV  may  require  a  lower  TIT  to 
prolong  its  service  life. 

Some  limitations  are  inherent  in  the  bypass  engine  type. 
A  turbofan  with  mixed  exhausts  may  be  bypass  ratio  limited. 
The  bypass  air  is  mixed  with  the  core  air  at  the  splitter 
plate,  where  the  static  pressures  must  match.  If  the  Mach 
number  of  the  bypass  air  is  too  low,  the  losses  are  too  great. 
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Reasons  for  mixing  the  exhaust  streams  are  to  reduce  engine 
noise  and  to  reduce  the  IR  signature  by  reducing  the 
temperature  of  the  core  air  by  mixing  v’ith  the  cold  bypass 
air.  There  is  an  improvement  in  performance  ideally  if  mixing 
occurs  internally  when  there  is  a  large  temperature  difference 
between  the  two  streams.  [Ref.  8:pp.  165,  172] 

The  present  representation  of  the  wave  rotor  component  is 
somewhat  limited.  The  modeling  associated  with  instantaneous 
port  opening  and  closing  may  require  further  investigation. 
Slow  port  openings  may  provide  more  efficient  wave  rotor 
operation  at  length  to  width  ratios  below  three  and  stagnation 
pressure  losses  are  predicted  to  increase  as  the  length  to 
width  ratio  increases  from  three  to  six  [Ref.  l:p.  236]. 

In  performing  cycle  calculations,  care  must  be  taken  to 
use  realistic  inputs  for  the  wave  rotor  pressure  ratio. 
Reference  1,  p.  18,  suggests  some  practical  limits.  As  the 
static  pressure  ratio  (SPR)  increases,  the  efficiency  of  the 
wave  rotor  decreases.  For  good  efficiency,  total  pressure 
ratios  should  be  about  2.0.  In  the  gas  generator  model,  or 
pressure  gain  combustor,  the  total  pressure  ratio  (TPR)  should 
range  from  1.1  to  2.0,  corresponding  to  a  total  temperature 
(TTR)  range  between  1.5  to  3.8  [Refs  14;  15:p.  707]. 

The  treatment  of  the  wave  rotor  as  an  ideal  component  in 
the  engine  cycle  is  a  limitation  which  eventually  must  be 
removed.  The  engine  SFC  and  specific  thrust  depend  greatly  on 
the  component  efficiencies  used  in  the  performance 
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calculations.  Their  affects  are  not  always  obvious.  A  key 
parameter  in  turbofan  engine  design  is  the  bypass  ratio.  As 
the  bypass  ratio  increases,  specific  thrust  decreases  but  SFC 
decreases  to  a  minimum  and  then  increases  [Ref.  8:p.  163]. 
The  optimum  arrangement  becomes  dependent  on  the  efficiencies 
attained  by  the  individual  components. 

B.  RESULTS 

A  study  was  conducted  using  WRCOMP  and  ENGINE  to  compare 
a  turbofan  with  mixed  exhausts  and  incorporating  a  wave  rotor, 
with  the  baseline  engine  operating  with  the  same  turbine  inlet 
temperatures.  Table  III  lists  the  input  parameters  used  in 
the  comparison.  A  list  of  the  component  efficiencies  is  given 
in  Appendix  A. 

In  maintaining  a  constant  TIT,  the  effect  of  adding  a  wave 
rotor  component  to  the  turbofan  was  to  change  both  the  SFC  and 
the  specific  thrust.  The  results  obtained  for  an  engine 
designed  with  a  wave  rotor  component  operating  in  a  gas 
generator  mode  are  shown  in  Figure  12. 

For  each  run  the  fan  pressure  ratio  (Tif)  was  fixed  and  the 
compressor  pressure  ratio  (it,.)  was  varied.  Initially,  the 
bypass  ratio  was  allowed  to  vary  to  satisfy  the  equal-static- 
pressure  requirement  at  the  splitter  plate,  for  a  mixed 
exhaust.  The  calculated  bypass  ratio  for  all  the  test  cases 
ranged  broadly  from  zero  to  three.  The  overall  pressure  ratio 
was  varied  by  fixing  the  fan  and  the  wave  rotor  pressure 
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TABLE  III 


RANGES  OF  INPUT  PARAMETERS  IN  ENGINE  COMPARISON 


Mo  =0.79 

h  =  35000  ft 

TIT  =  1860 

M^P-  =  0.4 

BASELINE 

TURBOFAN 

TURBOFAN  WITH 
WAVE  ROTOR 

p  (bypass  ratio) 

0.5 

varied 

2.3 

2.3  -  3.8 

10.87 

mtam 

^OPF. 

25.0 

15  -  32 

— 

1.6 

static  Press. 
Ratio 

— 

0.4 

rw 


ratios  and  varying  the  compressor  pressure  ratio.  The  fan 
pressure  ratio  was  then  increased.  At  each  fan  pressure 
ratio,  a  specific  range  of  values  was  obtained  for  the 
compressor  pressure  ratio  for  which  valid  solutions  existed. 
As  the  fan  pressure  ratio  was  increased,  the  feasible  range  of 
values  for  the  compressor  pressure  ratio  decreased  (the  best 
performance  range  for  the  overall  pressure  ratio  was  from  22 
to  26)  ,  and  the  bypass  ratio  decreased.  For  example,  for 
approximately  the  same  TIT  =  1860  and  ^opj,  =  23.6  and  for  = 
2.3  then  =  6.4  and  p  =  1.802;  for  =  2.8  then  =  5.3  and 
P  ^  1.144.  With  the  wave  rotor  incorporated,  the  SFC  was 
consistently  lower  than  the  baseline  turbofan.  The  specific 
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0  1 


Figure  12.  Effect  of  Varying  PRHC  and  PRLC  for  a 
Turbofan  with  a  Wave  Rotor 

(M,,  =  0.79,  h  =  35,000  ft,  AMACH8  =  0.4,  P  varies) 


thrust  was  also  higher  than  the  baseline  engine  except  for  the 
case  of  rvj  =  2.3. 

Since  a  large  bypass  ratio  may  not  be  acceptable  for  a 
cruise  missile  (a  larger  bypass  ratio  would  require  a  larger 
fan  size  which  may  not  meet  the  packaging  requirements) ,  the 
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bypass  ratio  was  set  to  p  =  0.5,  and  the  runs  repeated.  The 
results  are  given  in  Figure  13.  A  much  smaller  range  of 
feasible  compressor  ratios  resulted. 


Table  IV  shows  the  trade-off  between  %  change  in  SFC  and 
specific  thrust  (referred  to  the  baseline  turbofan  engine)  as 


the  fan  pressure  ratio  was  changed. 


TABLE  IV 

PERCENT  DIFFERENCE  FROM  BASELINE  ENGINE 


FAN  PRESSURE 
RATIO 

%  DIFFERENCE 
SFC 

%  DIFFERENCE 
F/m 

2.3 

-15.62 

-16.83 

2.8 

-10.22 

+  9.01 

3.0 

-  8.05 

+19.38 

3.2 

-  7.33 

+23.79 

3.4 

-  5.65 

+32 . 07 

3.6 

-  4.09 

+39.76 

3 . 8 

-  3.12 

+44.86 

Since  it  was  not  obvious  that  the  improved  performance  was 
due  to  the  wave  rotor,  the  baseline  engine  performance  was 
calculated  as  the  fan  pressure  ratio  was  varied  from  2.0  to 
3.5.  The  results  are  shown  in  Figure  14.  It  can  be  seen  in 
Figure  14  that,  as  the  fan  pressure  ratio  was  increased,  the 
SFC  increased  very  slightly  and  the  specific  thrust  was  also 
almost  constant.  Thus  the  improved  performance  in  Table  IV 
was  a  direct  result  of  the  wave  rotor. 
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Figure  14.  Effect  of  Varying  the  Fan  Pressure  Ratio 
for  the  Baseline  Turbofan 

Table  '  summarizes  the  results  of  the  comparison  of  the 
baseline  turbofan  with  an  engine  that  incorporates  a  wave 
rotor  component.  The  baseline  performance  is  shown  compared 
to  the  wave  rotor  turbofan  with  the  same  turbine  inlet 
temperature  (TIT)  ,  the  same  bypass  ratio  (P)  ,  or  the  best 
specific  thrust  (ST) .  Of  interest  to  a  missile  application, 
a  9%  increase  in  specific  thrust  and  a  7.3%  decrease  in  SFC 
was  obtained  by  incorporating  a  wave  rotor  into  the  baseline 
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turbofan  and  increasing  the  fan  pressure  ratio  to  2.8,  while 
maintaining  P  =  0.5. 

Table  V  also  infers  other  options  available  in  particular 
applications.  For  example,  if  the  capture  area  was  not  a 
design  limitation,  a  higher  bypass  ratio  (P  =  0.91)  and  fan 
pressure  ratio  (itf  =  3.0)  could  improve  the  specific  thrust 
(ST)  by  16.5%  and  decrease  the  SFC  by  6.2%  while  operating  at 
the  same  TIT  =  1860.  Similar  improvements  in  specific  thrust 
(ST)  and  SFC  were  obtained  for  lower  values  of  TIT  at  the 
specified  Jif,  as  shown  in  the  "best  ST"  entry  in  Table  V. 
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TABLE  V 


COMPARISON  OF  BASELINE  TURBOFAN  WITH  WAVE  ROTOR  TURBOFAN 


CONFIG 

■ 

^opr 

P 

TIT 

°R 

B 

SFC 

Ibm 

hr 

Ibf 

% 

DIFF 

ST 

% 

DIFF 

SFC 

BASE¬ 

LINE 

2.3 

10.8 

25.0 

.50 

1860 

30.2 

.832 

— 

B 

WR: 

1 

SAME 

TIT 

2 . 3 

D 

23.6 

1.8 

1856 

25.6 

.727 

-15 

-12  1 

SAME 

P 

2.3 

M 

18.4 

.50 

1443 

B 

.744 

-13 

-11 

BEST 

ST 

2.3 

B 

16.2 

.03 

1267 

26.2 

.744 

-13 

-11 

SAME 

TIT 

2.8 

5.3 

B 

B 

1876 

32.6 

.762 

8.0 

-8.4 

SAME 

P 

2.8 

B 

20.6 

.53 

1624 

32.9 

.771 

9.0 

B 

BEST 

ST 

2.8 

4.0 

17.9 

.06 

1409 

33.1 

.773 

9.6 

B 

SAME 

TIT 

3.0 

m 

23.5 

.91 

1859 

35.2 

.780 

16.5 

-6.2 

SAME 

P 

3.0 

B 

MNj 

.48 

1665 

B 

.782 

17.3 

-6.0 

BEST 

ST 

3.0 

4.0 

19.2 

.17 

1511 

35.6 

.780 

17.8 

-6.2 

SAME 

TIT 

3.8 

3.9 

B 

1874 

44.3 

.832 

0 

SAME 

P 

3.8 

B 

24.3 

.53 

1923 

44.2 

.832 

46.5 

0 

BEST 

ST 

3.8 

3.8 

23.1 

.36 

1825 

B 

.832 

B 

0 

VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


In  the  present  study,  the  computer  codes  ENGINE  and  WRCOMP 
were  installed  and  successfully  operated  on  a  VAX  computer 
system.  Modifications  were  made  to  enhance  the  utility  of  the 
programs  and  the  results  reported  by  the  codes'  author,  A. 
Mathur,  were  reproduced. 

Comparisons  between  the  predictions  of  ENGINE  and  ONX  (by 
J.  Mattingly)  led  to  the  conclusion  that  the  use  of  constant 
specific  heats  in  ONX  leads  to  less  accurate  results  than  can 
be  obtained  with  ENGINE,  which  includes  real  gas  behavior 
throughout.  Indeed,  the  predictions  of  ONX  were  found  to  be 
highly  sensitive  to  the  input  values  of  specific  heats. 

In  first  attempts  to  use  the  codes  to  examine  the 
potential  benefits  of  incorporating  a  wave  rotor  component  in 
a  mixed-exhaust  by-pass  fan  engine,  it  was  found  that  benefits 
were  available  which  could  not  be  obtained  by  varying  the 
conventional  engine  components.  For  example,  while 
maintaining  the  same  by-pass  ratio,  a  9%  increase  in  specific 
thrust  and  7.3%  decrease  in  SFC  were  predicted  if,  in  addition 
to  incorporating  the  wave  rotor,  the  fan  pressure  ratio  was 
increased  from  2.3  to  2.8.  Increases  in  specific  thrust  of 
almost  47%  could  be  obtained  with  no  change  in  SFC  by 
increasing  the  fan  pressure  ratio  to  3.8. 
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These  results  show  that  engines  incorporating  wave  rotors 
can  have  improved  performance  over  conventional  gas  turbines 
which  are  limited  to  the  same  turbine  inlet  temperature 
levels.  Further  extensions  of  the  ENGINE  code 

(1)  to  include  other  engine  types 

(2)  to  include  friction  and  heat  transfer  in  the  wave  rotor 
simulation 

are  recommended.  Also,  studies  to  examine  the  effects  of 
varying  static  pressure  ratio  and  stagnation  pressure  ratio 
across  the  wave  rotor,  and  of  increasing  turbine  inlet 
temperature,  should  be  carried  out. 

Finally,  an  experimental  program  to  validate  the  wave 
rotor  flow  predictions  and  performance  levels  used  in  cycle 
studies,  is  requirec. 
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APPENDIX  A 


PROGRAM  _ENGJENE 


A.  PROCEDURES 

1 .  Log  on  to  the  VAX 

2 .  At  the  DCL  prompt  type : 


$gksetup . initialize  GRAFkit  software 

$uis  . sends  plots  to  screen 

$run  ENGINEIA . runs  ENGINE  program 


. select  options  from  screen 

3.  To  make  changes  in  ENGINEIA. FOR: 


$gksetup  . only  required  once 

$uis  . 

$edit  ENGINEIA. FOR  ....to  enter  edit  mode 


*C  . to  enter  full  page  editing 

4  «  «  •  • 

<CNTRL/Z>  . exits  edit  mode 

♦exit  . saves  changes,  DCL  prompt  returns 

$FORTRAN  ENGINEIA . compiles  FORTRAN  code 


$link  ENGINEIA, 'GKL*  ..links  ENGINEIA  and  GRAFkit 
$run  ENGINEIA 

. follow  screen  instructions 

4.  To  send  your  results  to  a  laser  printer: 


$gksetup  . initialize  GRAFkit 

. environment 

$ln03  . sends  output  to  laser  printer 

$define  gk_out  plot.dat  .sends  output  to  <fn>.<ft> 

. where  plot.dat  is  your  output 

. file  name 

$ln03s_sp  . output  to  laser  for  square 

. plot 

$run  ENGINLiA  . follow  screen  instructions 


$print/passall/que=ln03  plot.dat 

. plot.dat  is  the  same  file  above 

Note:  Since  the  results  are  sent  to  an  output  file, 
"plot.dat",  they  will  not  appear  on  the  screen. 
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B.  POLYTROPIC  EFFICIENCY  SAMPLE  CALCULATIONS 


The  following  sample  calculations  use  the  ENGINE  output 
sample  results  from  Section  C  in  Appendix  A. 

Compressor  efficiency: 

In  terms  of  cycle  parameters, 

ideal  worx  interaction  A(l) 

' actual  work  interaction 

for  a  given  pressure  ratio.  In  terms  of  cycle  parameters. 


1f  c-l 


Conicresocr  polvtropic  efficiency; 

For  a  given  differential  pressure  ratio, 

^  ideal  work  interaction 
actual  work  interaction 

The  relationship  with  compressor  efficiency  is., 

Yc-l 


so  that 


A(3) 


A(4) 


In  X 


Yc-l 
Yc  • 


•  In  71  p 


A(5) 


or 


Yc-1  .  In  7t^ 

Y  c  In  X  ^ 


A(6) 
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For  the  fan  (LPC) : 

From  ENGINE1.DAT: 


V 


c 


71 


f 


1 .4004 
11.7 
5.1 
582 
444 


2.2941 

1.3108 


giving,  with  Eg.  A(6), 
e,  =  0.8772 


For  the  high  pressure  compressor  (HPC) : 


From  ENGINE1.DAT: 


Yc 


7t 


C 


I. 3870 
127.0 

II. 7 


10 . 8547 


giving,  with  Eg.  A(6) , 
e,  =  0.9132 


Turbine  efficiency; 

For  a  given  pressure  ratio, 

actual  work  interaction 
^  ideal  work  interaction 

In  terms  of  cycle  parameters. 
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n  t  = 


1  -  T 


1  - 


Turbine  polvtropic  efficiency: 

For  a  given  differential  pressure  ratio. 


A(8) 


actual  work  interaction 

Q  =  I  .Ill  — I.  — , 

ideal  work  interaction 


A(9) 


The  relationship  to  turbine  efficiency  is 


t 


gt  •  <Y 

Vt 


A(10) 


so  that 


In  T  j.  =  Sc 


In  71 


A(ll) 


^ 

Y  t-1  In  71  . 


A(12) 


For  the  high  pressure  turbine  (HPT) ; 
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From  ENGINE1.DAT: 


Ye 


Tt 


t 


T 


t 


1 . 3436 
22.8 
123.2 
1295 
1860 


=  0.1851 
=  0.6962 


giving,  with  Eg.  A(12), 

=  0.8392 

For  the  low  pressure  turbine  (LPT) : 


Y  t 


n 


t 


T 


t 


1.3612 
11.6 
22 . 8 
1104 
1295 


0.5088 

0.8525 


giving,  with  Eg.  A(12)  , 
e^  =  0.8899 


ENGINE  was  modified  to  calculate  the  polytropic  efficiency  for 
comparison  with  ONX. 
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SAMPLE  RESULTS 


NOTATION  COMPARISON  FOR  ONX  AND  ENCINE 


3 


4  4b  4c  t' 6  6  /  9 


Layout  of  baseline  Turbofan 


Station  l^cnti^ioation : 


0-t  Inlet  Section 

1~2  Low  Trossurc  Compressor 

2~3  fressure  Compressor 

3- 4  Msln  Burner 

4- 5  High  Pressure  turbine 

5- 6  Constant  Pressure  Miser 

6- 8  Low  Pressure  turbine 

i-7  Bypass  0»ict 

7.0-9  Constant  Area  Mixer 
9-10  Afterburner 
10-12  Hoztie 


Nave  Rotor  as  Gas  Generator 

6 1  a t ion  X^ant 1 ^ioB  t ion : 

0-1  Inlet  section 

1- 2  Low  Pressure  Compressor 
f  -9  MIfh  Pressure  Compressor 
0-4  V»ve  Potor  s«4  Mtin  Burner 
4“9  Nl^h  Pressure  t»nb1r» 

5- B  Constant  Pressure  Mixer 

6- 0  Low  Pressure  Turbine 

2- 7  Bypass  Duct 

7,  0-9  Constant  Aree  Mixer 
•-|0  Afterburner 
10-12  Mottle 


2.  BASELINE  ENGINE  CALCULATIONS 
a .  DATA  INPUT 
Fiji ght  parameters: 

Mo  =  0.79 
h  =  35000  ft. 

Pso=  3.468  psi....\  (calculated  within 

I  the  program) 

Tso=  394.1  °R _ / 

Design  Choices: 

Tf.ei  =  520  °R 
FLHV  =  18000  BTU/lbm-R 
LPT  exit  Mach  number  (AMACH8)  =04 
Tif  =  PRLC  =2.3 
Tto  =  PRHC  =  10.8696 


Component  Figures  of  Merit; 

r\  f  =  ETALC  =0.86 
Ti  ^  =  ETAHC  =0.88 

T)  f.pj.  =  ETAHT  =0.87 
T]  Ipp  =  ETALT  =0.90 
r\t>  =  ETACC  =0.96 


Maximum  Component  Temperatures; 

T„„oc  =  I8  60  °R 
=  34  00  °R 
T„3^,  =  208  5  °R 
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Aircraft  System  Parameters; 

XI  =  HPT  coolant  air  =  0.0 
Bypass  ratio  =  P  =  0.5 


Component  losses: 

Xgp  =  BPLOSS  =0.02 
\cc  =  CCLOSS  =0.03 
Xj^  =  ABLOSS  =  0.02  if  AB  off 
=  0 . 06  if  AB  on 


Miscellaneous : 

=  PIMAX  =0.97 
=  PIDMAX  =  0.97 
VCOEFF  =0.97 
With  the  Wave  Rotor; 

’‘SFR  =  SPR  =  0.4 
^TPR  ~  TPR  =  1.6 

RREF=  reference  density  in  WR  =  3.6  kg/m^ 
Wr  exit  Mach  number  =  AMEX  =  0.9 
Ttf  =  PRLC  =  3.8 

Ti,  =  PRHC  =  4.1 
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b 


ONX  OUTPUT 


ON -DESIGN  CALCULATIONS 
rURBOFAN  ENGINE  WITH  MIXED  EXHAUST 

*••***•******♦♦•*»*•*•  INFUT  DATA  ********•**•*■»**•»**•* 


MACH  NO  =  .790 

AL  r  HA 

= 

.  00 

ALT  (FT)  =  35000. 

ri  c 

= 

8 . 30 

TO  (R)  “  394.10 

FI  D  (MAX) 

=■ 

.97 

FO  (ESI  A)  =  3.46Bi;i 

FI  P 

= 

.97 

DENSITY  =».  00073024 

FI  N 

-97 

(SLUG/CLIFT  ) 

EFnCIENCY 

CF  C  =  .043  PTU/LPM-R 

PURNCR 

= 

.96 

CF  T  -  .265  BTU/LDM-R 

MF.(  H  HI  FR 

1  -00 

GAHMA  C  -  1.400 

MECH  10  FR 

1  .  OO 

GAMMA  T  1 . 350 

LF  CONI  R  (FArj) 

= 

.GO 

(EC  ) 

TT4  MAX  -  106';).  (F;> 

HF  COMEF 

- 

.91 

(ECU) 

H  -  FL'EL  (EiTU/LPM)  IPOOO. 

HF  TUrjPlMr 

= 

.04 

(  E  T  H  ) 

CTO  =  .Cii'OO 

LI  TURAINE 

= 

.R7 

(ETL) 

COOL  I  MET  AIR  HI  =  .00  '/. 

FUR  mec:h  eff 

^  1 

.  00 

COOLIN(-i  AIR  «2  -  .OO  V. 

BLEED  AIR 

= 

.  00 

Fo/pf?  t-  1.00 

***  MIXER  PI  MIXER  MAX  -  1.00 

********  *^*«*»****»»-ini»-»*  RESUL,  rS 


TAU  R 

1.125 

A/'  =  973. 

1  Fi/Br.c: 

FI  R 

1 . 509 

Vo  =  76R . 

7  FI  /F:FC 

FI  U 

.970 

MASS  FLOW 

®  1 

LEn-l/SEU 

TAU  L 

5.138 

AREA  rCRO 

=  10. 95''^ 

SOE-T 

FUR  TO 

.'■0  KN 

At  EA  rERO« 

=■  lO.fX'7 

SOFT 

TAU  C' 

- 

1.312 

PT5'  /F  A  = 

3,367 

EtA  C  ' 

e 

.8604 

T15  /I0 

1 .476 

FI  C 

25.00 

TAU  Ml  = 

1  . K> 

FI  C' 

= 

8.300 

TAU  M2  = 

1  ,  !;u  *00 

lAU  C' 

r 

1.312 

TAU  M  =  1 

•  (.>'  u  >0 

F.TA  C 

ST 

.B604 

FI  M  =  1 

•  C'00'.‘ 

FI  CH 

= 

10.870 

M5 

•  4(,u  >0 

TAU  0-1 

=- 

2. 108 

M5  = 

•  or''.'0 

EtA  CH 

.0824 

M6  = 

•  0‘'*or> 

FI  TH 

=: 

.  1753 

AS'  /A5i  - 

•  00'  ”  ’ 

TAU  TH 

s- 

.  6055 

GAMMA  M-  1 

,  37(0(:' 

ETA  TH 

=- 

.0669 

CF  M 

-867)0 

FI  TL 

= 

.6?9i;i 

ETA  TL  = 

.0  793 

lAU  IL 

s- 

.9014 

PT9/FO  =  3 

.0694 

UI  I  HO' IT 

r 

= 

.  Tl  11  6 

FO 

= 

.  IT  1  1  6 

FXM 

s 

42.415  LPr/LBM/S 

S 

=S 

.9834  1/HR 

T9/T0 

ST 

2 . 0535 

V9/V0 

= 

2.7434 

M7/M0 

= 

1  .9615 

t 

A7/A0 

ts 

.7400 

A9  /  AS 

- 

I  .2571 

THRUST 

= 

0403.  LPF 
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ENGINE  OUTPUT 


rFR>2.3  OFR-25.0  TtT-ie60.  ABTEnF-3400.  MACH  H0‘ 

STATION  IDENTt  ENGINE  INLET,  STATIONS  0-1 

TI-  394.  TO-  444.  PI-  3.5  PO-  5.1  FI-O.OOO 

STATION  IDENT:  LP  CONPRESSOR,  STATIONS  1-2 

TI-  444  .  TO-  582.  PI-  5.1  PO-  11.7  PI-. 000 

ETALC-O.eeoo  GBAR-1.4004 

ELPC-0.e757 

STAION  IDENT:  HP  COMPRESSOR,  STATIONS  2-3 

TI-  582.  TO-1206.  PI-  11.7  PO-127.0  PI-. 000 

ETAHC-0.0800  GBAR-1.3e70 

EHPC-0 .9144 

STATION  IDENT:  MAIN  BURNER,  STATIONS  3-4 

TI-1206.  TO-1860.  PI-127.0  PO-123.2  PI-.OOO 

STATION  IDENT:  HP  TURBINE,  STATIONS  4-5 

Tl-1860.  TO-1295.  PI-123.2  PO-  22.8  Pl-.OlO 

ETA-0.8700  GUAR-1. 3436 

ET-0.8379 

STATION  IDENT;  PRE  LPT  MIXING,  STATIONS  3,5-6 

Tl-1206.  T2-1295.  TO-1295.  PO-  22.8  Pl-.OOO 

STATION  IDENT:  BYPASS  DUCT,  STATIONS  2-7 

TI-  582.  TO-  582.  PI-  11.7  PO-  11.5  PI-.OOO 

STATION  IDEM;  LP  TURBINE,  STATIONS  6-8 

TI-1295.  Tr-M04.  PI-  22.8  PO-  11.6  Pl-.OlO 

ETA.0.9000  OtAR-1.3612 

ET-0.8901 

STA'^ION  IDENT;  CONST.  AREA  MIXER,  STATIONS  7,8-9 

Tl-1104.  t2-  582.  TO-  939.  PI-  11.6  P2-  11.5 

PI-. 010  F2-.000  PO-.007  Ml-0.40  H2-0.37  MO- 

STATION  IDENT:  AFTERBURNER,  STATIONS  9-10 

TI-  939  TO-  939.  PI-  11.2  PO-  11.0  PI-. 007 

STATION  IDENT:  EXHAUST  NOZZLE,  STATIONS  11-12 

TI-  939.  TO-  697.  PI-  11.0  PO-  3.5  PI-. 007 

GLOCO-1.3948  PRN-. 31661  AEXIT-1294.1 

PERFORMANCE  PARAMETERS 

AMACNS-  0.400 


0.79  ALT-35000. 

PO-.OOO  VO-  769. 

PO-.OOO 

FO-.OOO 

PO-.OlO 

FO-.OlO 

P2-.010  PO-.OlO 

FO-.OOO 

PO-.OlO 

PO-  11.2 
1.4  0  ARAT-0.38 

FO-.007 

VO-1728.  RO-1.34 


6P. THRUST  -  39.18  SFC  -  0.832 
(LBF/LBM/S)  <LBM/HR/LBF) 


BETA-  0.500 


EXMACH-  1.336 


3 


TURBOFAN  WITH  WAVE  ROTOR:  ENGINE  OUTPUT 


rPR-3.8  OPR-24.9  tit-1860.  ABTEMP-3400.  EACH  NO-0.79  ALT-35000. 

STATION  IDENTi  ENGINE  INLET,  STATIONS  0-1 

TI-  394.  TO-  444.  PI-  3.5  po-  5.1  PI-0.000  PO-.OOO  VO-  769. 

STATION  IDENT!  LP  COMPRESSOR,  STATIONS  1-2 

TO-  682.  PI-  5.1  PO-  19.3  PI-.OOO  PO-  000 

ETALC-0.8600  GBAR-1.3992  ro-.OOO 

ELPC-0.8841 

STAION  IDENT:  HP  COMPRESSOR,  STATIONS  2-3 

TI-  682.  TO-1056.  pi-  19. 3  po-  79.2  Pl-.OOO  PO-  000 

ETAHC-O.aeoO  GBAR.1.3870 

EHPC-0.9022 

STATION  IDENT:  NAVE  ROTOR,  STATIONS  3-4 

TCCIN-  1199.0  TO-  1972.7  TREF-  2197.7 
PIN-  79.2  PO-  122.9  PREF-  183.1 
PlN-0.0000  FO-0.0165  GLOCO-1.371  IDBl-1 

STATION  IDENT:  HP  TURBINE,  STATIONS  4-5 

PO-  52.4  F1-.016  PO-.016 

ETA-0.6700  ..BAR-1.3263 

ET-0.B551 

STATION  IDENT:  PRE  LPT  MIXING,  STATIONS  3,5-6 

Tl-1056.  T2-1649.  TO-1649.  PO-  52.4  Pl-.OOO  P2-.016  PO-.016 

STATION  IDENT;  BYPASS  DUCT,  STATIONS  2-7 
TI-  682.  TO-  682.  PI-  19.3  PO-  18.9  Pl-.OOO  PO-.OOO 

STATION  IDENT:  LP  TURBINE,  STATIONS  6-8 

PI-.016  PO-.016 

CTA-0.9000  GBAR-1.3416 

ET-0.8e66 

STATION  IDENT:  CONST.  AR  MIXER,  STATIONS  7,8-9 

T2-  682.  TO-1133.  PI-  20.6  P2-  18.9  Po-  19.4 

F1-.016  r2-.000  PO-.Oll  mi-0.40  H2-0.19  HO-0.31  ARAT-0.74 

STATION  IDENT:  AFTERBURNER,  STATIONS  9-10 

Tl-1133.  TO-1133.  PI-  19.4  PO-  19.0  PI-.Oll  PO-.Oll 

STATION  IDENT:  EXHAUST  N022LE,  STATIONS  11-12 

lO-  737.  PI-  19.0  PO-  3.5  PI— .011  VO-2232  MO-l  6fl 

GLOCO-1.3e90  PRN.. 18289  AEXIT-1329.8  VO-2232.  MO-1.68 

PERFORMANCE  PARAMETERS 


AMACH8-  0.400 

SP. THRUST  -  46.22  SFC  -  0.856  BETA-  0.500  EXMACH-  1.678 
(LBF/LBM/S)  (LBM/HR/LBF) 
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D.  ENGINE  PROGRAM  LISTING 


PROGRAn  ENGINE 

DIHENSION  XARAYdlS)  ,¥ARA¥(12S)  .XARAy(12S)  ,BPR(125) 

REAL  XP1(50),TP1(S0),¥P2(50I 
REAL  ST.SPTHR.SrC.PRAT 
INTEGER  CONFIG, nCONPIG, lit 22, 11, J1 
CHARACTER*6  XILABEL 
CRARACTER*15  YlLABEL 
CHARACTER*!  y2LABEL 

COHNON/ETAS/ETALC.ETAHC, ETACC, ETAHT, ETALT, ETAAS, VCOEFP 

COHNON/LOSSES/DI FLOS , CCLOSS , ABL055 , BFLOSS 

COHHON/BURN/nFUEL , FLHV , TNAXCC , THAXAB , THAXHR 

COHMON/NORRl/HORKLC.KORKHC 

COMMON/CONST/C2 , G 

COnnON  CONFIG 

**••  DEFINITION  OF  PROGRAN  SHITCRES 

**••  CONFIG-Ol CONVENTIONAL  ENGINE^  CONFIG-1 iNAVE-ROTOR  ENGINE. 

*•••  KODE-Ol NON-AFTERBURNING  CASE)  RODE-1 l AFTERBURNER  LIT  CASE. 

•*••  RABET-OtBYPABS  RATIO  (BETA)  PRESCRIBED.  KABET-llBETA  CALCULATED. 
DATA  IDC, IDS, IDT, IDN/0. 0,0,0/ 

DATA  XPl/50*0.0/ 

DATA  YPl/ 50*0.0/ 

DATA  ¥P2/50*0.0/ 

I  FLAG-0 

INPUT  DATA  (GENERAL) 


INITIALIZE  THE  DATA  BASE : *************************************** • 
WRITE) 6 , * ) 'HOULD  you  LIRE  TO  BEE  THE  POSSIBLE  CONFIGURATIONS 

>  FOR  ENGINE17  '.1-YES  0-NO)' 

WRITE) 6 , * ) ’ HI t  <Return>  for  next  picture  and  to  continue.' 

READ) 5, ' (I) * )1DESIGN 
IF  )IDEEIG*’.E0.0)  GOTO  11 

CALL  LIB$SPAWN  ('RENDER  ENGl.UIS') 

CALL  LIB$EPAHN  ('RENDER  ENG2.UIS') 

CONTINUE 

11  HRITE(6,*)'1NPUT  DESIRED  CONFIGURATION:’ 

WR1IE(«, * ) 'CONFIG-(0-CONVENTlONAL  l-WAVE  ROTOR)' 

READ(5, ' ( I ) ' )CONFIG 
BCO'-TIG-CONFIG 

WFITC(6, *J 'EODE-(0-NON  A/B  1-A/B  ON)' 

BEAD(5,  'O'  )R0CE 

WRITE ( b, *)' SASCT-(0-ByFASS  RATIO  PRESCRIBED  1-BETA  CALCULATED)' 
READ(5,'(1)')RASET 

WRITE (6,*) 'WOULD  YOU  LIRE  TO  USE  YOUR  OWN  PARAHETERS  OR  IRE 

>  DEFAULT' 

WRITE(6,*)'  PABAHETERS7  ( 1-OWN, 0-DEFAULT ) ' 

READ) 5, ' ( I ) ' )IQ1 
IF  (IQl.EQ.l)  THEN 

CALL  FLIGHT) FNACH, ALT, PSO.TSO) 

C  WBITE(6,*)'CONFIG-', CONFIG 

IF  ;-0V"G.E3.1)  THEN 

CALL  NR8FECE(EFR,TPR.RREF,AHEX) 

END  IF 

C  WRITE) 6,*) 'CONFIG-', CONFIG 

CALL  DESLin(TFUCl ,FLRV,AnACH8,HFUEL,PRLC,PRHC,PRAT,TPR, 
>CONriG) 

CALL  CONFFOn ( ETALC , ETAHC , ETAHT , ETALT , ETACC , ETAAB ) 

CALL  TKAX)THAXCC, THAXAB, TNAXWR) 

CALL  ACSyS(CTO,THAXCC,Xl,BETA) 

CALL  COHLOSSjBPLOSb, CCLOSS, ABLOSS) 

CALL  COnFETA(FinAX,PIDnAX.VCOEFF) 


50 


o  n 


ELSEir  (IQl.EQ.O)  THEN 

CALL  BASEl(rHACH,ALT,PSO.TSO,TrUEL,rLnV,AHACnB,nrUEL, 
>PRLC,PRHC,  PRAT ,  ETALC . ETAHC , BTAHT , ETALT , ETACC , ETAAB  ,  EATDB , 
>THAXCC,TnAXAB,THAXHR,CTO,Xl.BETA,8PLOSS,CCLOSS,ABLOS8, 
>PinAX,PtDMAX,VCOErr,SPR,TPR,RREP,AMEX,CONriG,EODE> 

ENDir 

c - SET  INTIAL  CONSTANTS - 

NN-0 

CJ-778. IS 
G-32.174 
ro-D.o 

- END  OP  INPUT  DATA 

OPEN  (UNlT-9,riLE-'EN0INEl.DAT' , STATUS- ' NEW ' ) 

CALLING  SEQUENCE  OP  COMPONENT  SUBROUTINES  POR  CHOSEN  ENGINE  *** 
- CONPIGURATION 

MR1TE(6,»)'W0ULD  TOO  LIRE  TO  ITERATE  A  VARIABLE  ( 1-TES  O-NOl' 
READ(5,M1)'  )IANS1 
IP  (lANSl.EQ.OI  THEN 
GO  TO  6 
ELSE 

writ:;* . » ; ' rNP"T  the  desired  itsp-rtion  parameter  (!■  l)' 

HRITE(S,*J'  1.  BYPASS  RATIO' 

WBITE(6,»)'  2.  PAN  PRESSURE  RATIO' 

NRITE(6,»)'  3.  COMPRESSOR  RATIO' 

WRITE(6,*)'  4.  LPT  EXIT  MACH  NUMBER' 

READ(5,»)IANS4 

1ANSH4-IANS4 

WfiITE(6,») 'INPUT  MIN  VALUE  -' 

READ(S,»)8MIN 

WR1TE(6,»)' INPUT  MAX  VALUE 
REA0(5,«)SHAX 

HRITE(6,*)'INPUT  STEP  SHE 
READ (5,*) STEP 

0PEN(UN1T-14,PILE-'PL0T1.DAT' ,STATUS-'HEH’ ) 

IP  (IANS4.EQ.il  then 
X1LABEL-'  BETA  ' 


C 

WRITE(14, ' ( 2X,A,2X,A,2X,A)' )' 
EL5EIP  (IANS4.EQ.2)  THEN 
XILABEL-'  PRLC  ’ 

BETA' , ' 

ST 

p  r 
t 

SPC  ' 

c 

WRITE! 14, ' ( 2X,A,2X,A,2X,A) ' ) ' 
ELSEIP  (IAN54.EQ.3I  THEN 
X1LABEL-'  PRHC  ' 

PRLC' , ’ 

ST 

f  9 
* 

SPC  ' 

c 

WRITE! 14, ' (2X,A,2X,A,2X,A|' )' 
ELSEIP  (IANS4.EQ.4I  THEN 
XILABEL- 'AHACH8’ 

PRHC' , ’ 

ST  '  ,  ' 

SPC  ' 

c 

WRITE! 14, '(2X,A,2X,A,3X,A)' ) ' 
ENDIF 

AHACH8' , ' 

ST  '  ,  ' 

SPC  ' 

WRITE(14, ' (2X,A,2X,A,2X,A,2X,A,2X, 

A,2X,A6,2X 

,A,2X,A) ' 

) '  BETA 

>’  PRLC  PRHC  PRAT  TIT  ' ,X1LABEL, '  ST  8PC  • 

IC2-1 

C - START  LARGE  LOOP  ITERATION  POR  ITERATION  VARIABLE - 

DO  41  SVALUE-SHIN.SMAX.STEP 
IP  (IANS4.E0.1)  THEN 
BETA-SVALUE 

EL5EIP  (IANS4.E0.2)  THEN 
PRLC-SVALUE 
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BLSEir  (IANS4.EQ.3)  THEN 
PRHC-SVALUE 

ELSEir  (IANS4.EQ.4)  THEN 
AMACH8-8VALUE 
ENDir 

ir  (CONFIG. EQ.l)  THEN 
FRAT-FRLC*PERC*TFR 
ELSE 

PRAT-PFLC»PRHC 
END  IP 

- KRITE  STATION  RESULTS  TO  DEVICE  9  ON  PILE  ENGINE1.DAT - 

- ITERATION  PLOT  RESULTS  ARE  SENT  TO  14  AND  PLOTl.DAT - 

6  NRITE(9,1)FRLC,PRAT, TNAXCC , THAXAB , PHACR , ALT 

1  FORMAT ( 2X, 'PPR.' ,P3.l,2X, 'OPR- ’ , P4 . 1 , 2X , ' TIT-' , PS . 0 , 2X, 'ABTEMP-' , 
>P5.0,2X, 'MACH  NO-' .P4.2,2X,'ALT-' ,P6.0/) 

CALL  INLET(FSO,TSO,FO.PI,Tl,Pl,VO,FKACn,PIDHAX) 

CALL  C0MPR(PRLC,P1,T1 ,P2,T2,T2IS,F1,P2,IDC) 

CALL  C0MPH(PRHC,P2,T2.P3,T3.T3IB,P2,P3,IDC) 

IP(CONriG.EQ.l )  GOTO  10 

CALL  BURNER(P3,T3,P4,T4,P3,P4,1DB, RODE, CONFIG) 

GOTO  20 

10  CALL  ROTOR ( P3, P4,T3,T4, PREP, RREP,TREP,F3,P4,SPR,TPR, AMEX, TCCIN, 

>IDB, CONFIG) 

10  CALL  ROTOR ( PRAT, ?1, P3,P4,T3,T4, PREP, RREP,TPPr,P3,P4,SPR,TPR, AMEX, 

>TCCIN,1DB, CONFIG) 

20  CALL  TURBIN(P4,T4,P5,T5,TS18,P4,P5,B5,PRAT,PRLC,PRHT,PRLT,X1,BETA, 
>KASET, IDT) 

CALL  MIXCP(F3,PS,P3,T3,P5,T5,P6,P6,T6,XI,IDH) 

CALL  DUCT(P2,T2,P2,P7,T7,F7,H7) 

CALL  TURBIN(P6,T6,P8,T8,T8IS,P6,F8,H8,PRAT,PRLC,PRHT,PRLT,X1,BETA, 
>RASET, IDT) 

IP(RASET.EQ.l)  GOTO  30 

CALL  HIXCA(P7,P8,T7,T8,R7,Re,P7,P8,BBTA,AHACBB,PlHAX,T9,P9,P9, 
>AMACH9) 

GOTO  40 

30  CALL  I1IXCP(P7,P8,P7,T7,P8,TB,P9,P9,T9,BETA,1DM) 

40  CALL  BUBNEB(P9, T9, P10,T10,P9,P10,IDB, RODE, CONFIG) 

CALL  NOZZLC(P10,T10,P12,T12,T12IE,P10,PRLC,VEXIT,BXHACB,PSO,IFLAG) 

IP(IPLAG.EO.l)  GOTO  50 

CALL  PERP(AMACH8,P10,VO,VEXIT,BETA,SPTHR,SPC,EXMACH) 

WRITE ( 14, 3014 )  BETA, PRLC, PRHC, PRAT, T4 , SVALUE, SPTHR, SFC 
WRITE! 14, ' (2X,P6. 3,2X,P6.3,2X,F6. 3)' )8VALUE, SFTRR, SFC 
WFITE(6,*)'IC2-',IC2 
WRITE (6,*) 'ST-', SPTHR, 'SFC-', SFC 
XPMIC2)  -  SVALUE 
yPl(IC2)  -  SPTHR 
yP2(IC2)  -  SFC 

WRITE(6,»)'XPl-',XPl(IC2),'yPl-',yPl(IC2),’yP2-',yP2(IC2) 

IF  ( SVALUE. LE.SHAX)  IC2-IC2^1 

- BESET  THE  COUNTER  SWITCHES  EACH  LOOP 

IDC-0 
IDB-0 
IDT-0 
IDH-0 

CONFIG-HCONFIC 

- CHECR  PROPER  COUNTER  OPERATION - 

MRITE(6, •) 'IDC-' ,IDC, 'IDB-' ,IDB, 'IDT-' ,IDT, 'IDM-’ ,IDH 
NBITE(6,«)'CONFIG-',CONFIO 

41  CONTINUE 

c - end  large  loop  variation  FOR  ITERATION  VARIABLE - 
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END  IF 
ENDFILE  14 
REWIND  14 
CLOSE  (UNIT-14) 


CALL  GRAFKIT  PLOT  ROUTINE 


NR1TE(6,*)'N0ULD  YOU  LIRE  TO  PLOT  THE  RESULTS?  (l-YES  O-NO)' 
READ( 5, ’ (I ) ' )IC5 
IF  (IC5.E0.1)  THEN 
NOW  CALL  GRAPHICS  ROUTINE 

CALL  GRSPLOTl ( XPl . YPl , YPl , 1ANSN4 , IC2-1 ) 

ENDIF 


•**END  OF  USER-DEFINED  CALLING  SEQUENCE 


OPEN(UNIT-10,FILE-'ENGINE2.DAT' , STATUS-' NEW ' ) 

WRITE( 10, 3000)P6,P7,P8,Pll 
CLOSE(UNIT-IO) 

00  FORNATCSX, ’P8-' ,F5.1,5X, 'P?-' , FB . 1 , 5X, ' P8- ' , F5 . 1 , 5X, 'P11-' , 
>F5.1//) 

14  FORnAT(2X,F6.3,2X,F6.3,2X.P6.3,2X,F6.a,2X,F6.1,2X,Ffi.3,2X, 

>F6.3,2X,F6.3) 

GOTO  70 

50  WRITE(6,60) 

60  FORMAT(5X, 'NOZZLE  PRESSURE  RATIO  18  NEGATIVE'/) 

70  CL0SE(UNIT-9) 

STOP 

END 


SUBROUTINE  BASEl ( FHACH , ALT, PSO,TSO, TFUEL, FLHV, AHACHB , HFUEL , 
>PRLC,PRHC. PRAT , ETALC , ETAHC , ETAHT , ETALT , ETACC , ETAAB , EATDB , 
>TMAXCC,Th<.r.;.”  ,  TTO.XI  ,BETA,  oPLOSS  , CCLOSS  ,  ABLOSS  , 

>Pl(1AX,FIDnAX,VCOEFF,8PR,TPR,RREr,AnEX,  1  CONFIG ,  IRODE ) 

COMMON  CONFIG 
-FLIGHT  PARAMETERS 
FMACH-0.79 
ALT-35000. 

PSO-3.46e 
TSO-394 . 1 
-WAVE  ROTOR  SPECS 
EPR-0.40 
TPR-1 .60 
RREF-3 .6 
AMEX-0.9 

- TCCIN  IS  SOLVED  FOR  IN  ROTOR  ROUTINE - 

TCCIN-1213.3 
-DESIGN/FUEL  CHOICES 
TFUCL-520. 

FLHV-ieOOO. 

AMACH8-0 . 4 
HrUEL-TFUEL»0.5 
IF  (ICONFIO.EQ.l)  THEN 
PRLC-3.8 
PRHC-4.1 
ELSE 

PRLC-2.3 

PRHC-10.8696 

ENDIF 

IF  (ICONFIG.EQ.l)  THEN 
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PRAT-l>Rt.C»PRHC*TPR 

C  MR1IE(6,*) 'PRAT-' , PRAT, 'TPR-' ,TPR, ’ CONFIG-' , ICONPIO 

ELSE 

PRAT-PRLC*PRHC 

ENDir 

C - COMPONENT  POM 

ETALC-O.ee 

ETARC-0.88 

ETAHT-0.87 

ETALT-0.90 

ETACC-0.96 

IP  (IRODB.CQ.O)  THEN 
ETAAB-1 .00 
ELSE 

ETAAB-0.96 
END  IP 


C - MAX  COMPONENT  TEMPERATURES  (R) 

C  NOTE!  TMAXCC  AND  TMAXWB  GET  RESET  TO  VALUES  CORRESPONDING  TO 

C  THE  NAVE  ROTOR  ANALYSIS  I  TMAXCC-TREF,  TMAXNR-TRBP 

TMAXCC-ieSO.O 
TNAXAB-3400 . 

TNAXWR-2085 

C - AIRCRAFT  SYSTEM  PARAMETERS 


CTO-0.0 

WBITE( 6 , • ) ’ INPUT  t  COOLING  OH  BASE  IT  ON  TMAXCC:' 

WRITE (6,*)'  (I-INPUT  FOR  XI  0-BASE  XI  ON  TMAXCC)' 

BEAD(5,*)ICOOLl 

IP  ( ICOOLl .EO.O)  THEN 

IP  ( TMAXCC. GT. 2 40u)  THEN 

XI-2.C*(TMAXCC-2400.0)/l6000.0 

ELSE 

XI-0.0 

ENOIP 

ELSE 

WRTTE(S, «' 'INPUT  XI' 

READ(5,*)XI 
END  IP 
BETA-O.S 

C - COMPONENT  LOSSES 

BPLOSS-0. C2 
CCLOSS-0.03 
, ’  ( IRODE.EQ.O)  THEN 
ABLOSS-0.02 
ELSE 

ABLOES-0.06 
END  IP 

C - COMPONENT  EPPICIENCIES 

PIMAX-0.97 
PIOMAX-0.97 
VCOEf  P-fj  .  , 

RETURN 

END 

C— 

C  INPUT  DATA  BASE  SUBROUTINES  .............. 

C - 

SUBROUTINE  PLIGHT) FHACH,ALT,PSO,TSO) 

DIMENSION  TT(Oi50) ,PP(0i50) 

DATA  TT/1 .0000,0.9931,0.9863,0.9794,0.9725,0.9656,0.9583,0.9519, 
>0.9540,0.9381,0.9313,0.9244,0.9175,0.9107,0.9038,0.8969,0.8901, 
>0.8832,0.8764,0.8695,0.8626,0.8558,0.8489,0.8420,0.8352,0.8233, 
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>0.6215,0.6146,0.8077,0.8009,0.7940,0.7672.0.7803,0.7735,0.7666, 

>0.7596,0.7529,0.7519,0.7519,0.7519,0.7519,0.7519,0.7519,0.7519, 

>0.7519,0.7519,0.7519,0.7519,0.7519,0.7519,0.7519/ 

DATA  PP/1. 0000, 0.9644, 0.9296, 0.6963, 0.6637, 0.8321, 0.8014, 0.7717, 
>0.7429,0.7149,0.6878,0.6616,0.6362,0.6115,0.5877,0.5646,0.5422, 
>0.5206,0.4997.0.4795,0.4599,0.4410.0.4227,0.4051,0.3860,0.3716, 
>0.3557,0.3404,0.3256,0.3113,0.2975,0.2643,0.2715,0.2592,0.2474, 
>0.2360,0.2250,0.2145,0.2044,0.1949,0.1858,0.1771,0.1688,0.1609, 
>0.1534,0.1462,0.1046,0.1329,0.1267,0.1200,0.1151/ 

WRITE(6,*) 'INPUT  PLIGHT  CONDITIONS' 

NRITE(6,*)'  INPUT  PLIGHT  MACH  NUMBBRi' 

READ(5,*)PHACH 

MR1TE(6,»)'  INPUT  ALTITUDE  (0-50000  PT)' 

READ(5,*)ALT 

IALT-ALT/1000 

TSO-( 459. 67+59. 0)*TT(IALT) 

PSO-14.696*PP(IALT) 

RETURN 

END 


SUBROUTINE  DE8LIH( TPUEL, PLRV,AMACHeA, HFUEL, PRLC , PRBC , PRAT.TPR, 

>icoNriG) 

COHHON  CONFIG 


MR1TE(6,+ 
WRITE( 6 , • 


READ(5,*’iPUEL 


WRITE( 6 , + 


READ(5,*)PLHV 


NRITE(6,* 


NRITFie  * 


INPUT  DESIGN/PUEL  LIMITATIONS i ' 
TPOEL-(TEMP  OP  PUEL  IN  R)' 


PLHV«(PUEL  LONER  HEATING  VALUE)' 


READ(5,*)ANACHeA 


AMACHe-(EXIT  MACH  AT  LPT)' 


PRLC-' 


REAO( 5, * )PRLC 
KRITE(6,*)'  PRHC»' 
BEAD(5,*)PRHC 
IP  (ICONPIG.EQ.O)  THEN 
TPR-1.0 
ENDir 

PRAT-PRHC+PRLC*TPR 

HFUEL-TPUEL+O.S 

RETURN 

END 


SUBROUTINE  COHPFOM ( ETALC , ETAHC , ETAHT, ETALT, ETACC , ETAAP ! 

KR1TE(6,*)'INPU1  COMPONENT  FIGURES  OF  MERIT! (1-OWN  0-STD)' 
KRITE(6,*)'(STDI  BTALC-0.86  ETAHC-0.88  ETAHT-0.87  BTALT-0.90 
>  ETACC-0.96  ETAAB-0.96)' 

READ( 5 , ' ( I ) ' ) ISUB2 
IP  (ISUB2.EQ.0)  THEN 
ETALC-0.6S 
ETAHC-O.eS 
ETAHT-0.e7 
ETALT-0.90 
ETACC-0.96 
ETAAB-0.96 
ELSE 

NRITE(6,»)'  input  ETALC- ' 

READ(5,»)ETALC 
NRITE(6,*)'  input  ETAHC-' 

BEAD(5,‘)ETAHC 
NRITE(6,»)'  input  ETAHT-' 
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ri;ad(5,»)«ta«t 
MRITE(6,»>'  INPUT  BTALT«' 

*EA0<5,»)BTAI.T 
NRITe(6,«)'  INPUT  BTACC-' 

READ(S,*)BTACC 
WRITE (6,*)'  INPUT  BTAAB>' 

RBAO(S,*>BTAAB 

BNt;iP 

RETURN 

END 

C 

SUBROUTINE  TNAX ( THAXCC , TNAXAB , TNAXWR ) 

WRITE* 6,*) 'INPUT  MAX  COMPONENT  TEMPERATURES' 

WRITE (6,*)'  INPUT  THAXCC- ' 

WRITE(6,»)'  THAXCC  IS  RESET  TO  THAXWR  IN  ROTOR  CONPIG' 

REAO(S,*)THAXCC 
WRITECS,*)'  INPUT  THAXAB-' 

READ(S,*)THAXA8 
WRITE* 6,»)'  INPUT  THAXWR- ' 

WRITE(6,«)'  THAXWR  IS  RESET  TO  TREE  IN  ROTOR  CONFIG’ 

READ* 5.*) THAXWR 
RETURN 
END 
C 

SUBROUTINE  ACSYS( CTO, THAXCC, XI ,BETA) 

WRITE(6,»)'VEH1CLE  system  PARAMETERS’ 

WRITE(6,»)’  input  CTO-  (TAaE-OPF  POWER)' 

READ<5,»)CT0 

WRITE(6,*)'  input  %  COOLING  OR  BASE  IT  ON  THAXCC:' 
MRITE(6,*)'  (1-INPUT  for  XI  0-BASB  XI  ON  THAXCC)' 

READ(5,*)ICOOLl 
IF  (ICOOLl.EQ.C)  THEN 

IF  ( THAXCC. GT. 2400)  THEN 

XI-2.0«(TMAXCC-2400.0)/16000.0 
ELSE 
Xl-0.0 
ENOir 
ELSE 

WRITE(6,«)’ INPUT  XI-» 

READ(S,* )XI 
ENDIF 

KRITE{6,*)'  INPUT  BETA-  (BY-PASS  RATIO)' 

RE  At)  (5,*)  BETA 
RETUHN 
END 

SUBROUTINE  COHLOSS (BPL055 , CCLOSS , ABLOS5 ) 

WRITE* 5, ♦) 'COMPONENT  LOSSES  ( PERCENTS )!< 1-OWN  0-STD)' 
WB1TE(6,«)'8TDi  BPLOSS-0.02  CCLOSS-0.03  ABLOSS-0.02  )' 
BEAD(5,’<I)')I05 
IF  (IQS.EQ.O)  THEN 
BFLOSS-O.Oa 
CCLOSS-0.03 
ABLOS6-0.02 
ELSE 

write<6,»)'  input  BPLOSS-' 

REAP(5,*)BPL0SS 
WRITE(6,*)'  input  CCLOSS-' 

REAO(5,*)CCLOSS 
WRITE<6,»)'  input  ABLOSS-' 
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I(EAD(S,*)ABLOSS 

ENDir 

RETURN 

END 

C 

SUBROUTINE  COHFETA( PINAX, PIDHAX , VCOEPP ) 

write(6,*)'hiscellaneous  component  choices :(1-0WN  O-STD)' 

NBITE(6,»|'(STDi  PIHAX-0.97  lIDKAX-0.97  VCOEPP-O . 97 ) ’ 

REAO(5,'tI)')IQ« 

IP  (I06.E0.0)  THEN 
PIHAX-0.97 
PIDHAX-0.97 
vcoErr-0.97 
ELSE 

NRITE(6,«)'  input  PIHAX-' 

REAP(5,*)PIMAX 
NBITE(6,»J'  INPUT  P1DMAX-' 

REA0(5,»)P1DMAX 

NRITE(6,*J'  INPUT  VCOEPP-  (NOZILE  EPP)' 

READ (S,») VCOEPP 
ENDIP 
RETURN 
END 
C 

SUBROUTINE  NRSPECS( SPR.TPR, RRET, AMEX ) 
nrite(6,*)'nave  rotor  SPECIPICATIONS' 

KRITE(6,*)'  INPUT  6PR-' 

READ( S,« )SPR 

NRITE(«,«)'  INPUT  IPR-' 

READ(5,*)IPR 

NR1TE(6,*)'  INPUT  RREP-  k9/«*3’ 

READ(5,«)RHEP 
NRITE(6,*)’  INPUT  AMEX-' 

READ(6,*)AnEX 

-TO  INPUT  TCCIN  CHANGE  NRSPECS  SUBR  AND  CALL - 

NRITE(6,*)'  INPUT  TCC1N-' 

READ(S,*)TCCIN 
RETURN 
END 


SUBROUTINE  INLET! PIN, TIN, PIN, PO, TO , PO, VO , PHACH, PIDMAX) 
COnnON/LOSSE5/DIPLOS,CCLOSS.ABLOES,BPLOSS 
COHMON/CONST/CJ , C 
HIN-ENTHPN(PIN,TIN) 

PHIIN-ENTBPN( PIN.TIN) 

GLOCIN-GLOC( PIN,TIN,HIN) 

PTIN-PIN*( 1 .♦(GLOCIN-1. )/2. •PHACH**’. ) *  * < CLOCIN/( GLOCIN-1 . ) ) 
TTIN-TIN*(l.*(GLOCIN-l.)/2.*PMACB**2. ) 

FID-PIOHAX 

IP(FnACR.GE.l.O)  PID-PlDHAX*(l.-0.07S*(FnACB-l. )**1.35) 

DIPLOS-1 ,-PID 

P0-PTIN*(1-DIFL08) 

TO-TTIN 

PO-0.0 

BO-HIN 

PHIO-PRIIN 

VO-PHACH*S0RT(GLOCIN*BGA8(PIH)*TIN*G) 

lfFITE(9,10) 

WRITE(9,20>TIN,TO,PIH,PO,P1H,PO,VO 
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10  rOnHAT(5X, 'STATION  lOENTt  ENGINE  INLET,  STATIONS  0-1'/) 

20  rOAHAT(5X, 'T1-' .rs.o.lx, 'TO-' ,r5.0,2x, 'PI-' ,r5.1,2x, 'SO-' , 
>r5.i,2x.'ri-',r5.3,2x,'ro-',r4.s,2x,'vo-' ,P5.o/l 

SETURN 

END 


SUBROUTINE  COnPR( PR , PIN.TIN, PO, TO, TOlS , PIN, fO, IOC ) 

COMtION/ETAS/ETALC ,  ETAHC ,  ETACC ,  ETAHT ,  ETALT ,  BTAAB ,  VCOSPE 

COHNON/LOSSES/DI PLOS , CCLOSB , ABLOS8 , SPLOSS 

CONMON/NOFKl/NORRLC.NORKHC 

C0HP10N/C0NST/C3 ,  G 

IDC-lDC+1 

IP(IDC.EO.l)  GO  TO  30 
C  »*«»»»»COHPUTE  HPC 

HlN-ENTHrN(PIN,TIH) 

PHI 1N-ENTRPN( PIN.TIN) 

PHI01S-PHIIN+RGAS( PIN)/CJ*ALOO(PR) 

CALL  PIIIPNfPIN.PBlOIS.TOIS) 

H01S-ENTHPN(PIN,T0IS) 
nO-HIN4^  (  HOI  S-niH  )/BTARC 
NOFPHC-RO-BIN 
CALL  HPN(PIN,BO,TO) 

PHIO-ENIFPN(PIH,TO) 

PO-PIN 
PO-PIN*PR 

GBAF-GAVG(P0,TIN,T0IS,PniIN,PH10IS) 

TEMPHCl-ALOG( TO/TIN) 

TENPHC2-AL0G( PO/PIN) 

EHPC- ( GBAR- 1 . 0 )/GBAR«TEHPHC2/TEHPHC1 
WR1TE(«,* ) 'TO-' ,TO, 'TIN-' ,TIN, 'PO-' ,PO, 'PIN-' , 

>  P I N , ' ETAHC- ' , ETAHC , ' GBAR- ' , GBAR 
KF 1 TE ( 6 , » ) ' EHPC- ' , EHPC 
NRITE(9,10) 

NF1TE(9,20)TIN,TO,PIN,PO,PIH,FO, ETAHC, GBAR, EHPC 
to  POFNAr(5X, 'STAION  IDBNTl  HP  CONPRESSOR,  STATIONS  2-3'/) 

20  PORMAT(5X, 'T1-' ,P5.0,2X, 'TO-' ,P5.0,2X, 'PI-' ,P5.1,2X, 'PO-' , 
>P5.1,2X,'PI-',P«.3,2X,'PO-',P4.3/,5X,' ETAHC- ' , P6 . 4 , 2X , ' GBAR- ' , 
>P6.4/,SX, 'EHPC-' ,r6.4/) 

RETURN 

“••COMPUTE  LPC 
30  BlN-ENTHrN(PIN,TIN) 

PHIIN-ENTFPNdlN.TIN) 

PHIOIS-PHIIN+HGAS(FIN)/CJ^ALOG(PF) 

CALL  PHirN(PIN,PHlOIS,TOIS) 

HOIS-ENTHPN(FIN,TOIS) 

H0-HIN4(H0IS-HIN)/BTALC 
NORRLC-HO-HIN 
CALL  HrN( PIN.HO.TO) 

PHIO-ENTRPN(PIN,TO) 

PO-PIN 
PO-PIN^PR 

GBAR-GRVG(P0,TIN,T0IS,PHIIN,PHI0IS) 

NB1TC(6,‘)'T0-' ,T0, 'TIN-' ,T1N, 'PO-' ,FO, 'PIN-' , 

>PIN, 'ETALC-' ,ETALC, 'GBAR-' ,GBAR 
TEMPLC1-AL0G( TO/TIN) 

TEnPLC2-ALQG( PO/PIN) 

ELPC- ( CBAF-1 . 0 )/GBAR*TEnPLC2/TEHPLCl 
C  NFITE(6,»)'ELPC-',ELPC 

MFITC(9,40) 

NRITE(9,SO)TIN,TO.P1N,PO,PIN,PO,BTALC, GBAR. ELPC 
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40  rOFMAT(5X, 'STATION  IDEMTi  LP  COMPRESSOR,  STATIONS  1-2’/) 

50  P0RMAT(5X, 'T1-' .r5.0,2X, 'TO-' ,r5.0,2X, 'PI-' ,r5.1,2X, 'PO-’ , 
>r5.1,2x, 'PI-' ,r4.3,2x, 'fO-' ,r4.S/,5X, 'RTALC-' , PS . 4 , 2X, ' OBAR-' , 
>P6.4/,5X, 'ELPC-',P6.4/) 

RETURN 

END 

C  *•«  *******ft****************«***ft««*ft***«****«4«««4««*«««*««*4«« 

c  •••«kAVC  ROTOR  8UBROOTlNE***»»»*»»»»»«‘**«***‘«o«*****»«»*** 


SUBROUTINE  BOTOR( PRATl , PTl , PIN, PO, TIN, TO, PREP, RREP, TREE , PIN, PO, 
>EPR,TPP,AMEX,TCCIN,IDB1 ,CONriG) 

COHMON/LOSSES/DlPLOS,CCLOSS,ABI,OSS,BPLOSS 
COMMON/BURN/HPUE', ,  PLH  V ,  THAXCC  ,  TMAXAB  ,  TNAXMR 
COMMON/CONST/C J . G 
DUMMY-0. 

GGUESS-1 . 3 
PO-f IN 

TSO-TMAXCC/{l.+(GGOESS-l. )/2 . ‘AMEX* *2 . ) 

RCT-0 

10  RCT-RCTtl 

HSO-ENTHPN(PO,TSO) 

GLOCO-GLOC(PO,TSO,HSO) 

PO-PIN*TPR*(l.-CCLOSS) 

PSO-PO/( 1 .♦(GLOCO-1. )/2.«AMEX**2. )»*(GLOCO/(GLOCO-l. ) ) 
PREP-PSO/SPR 

C - CONVERT  DENSITY  PROM  kg/»‘?  to  lb«/ft*3 - 

TREr-PBEP/RGAS(PO)/RREP*2306.65e7 
STR-SPR**( (GLOCO-1. )/GLOCO) 

TSO-STR*TREP 

TO-TSO‘(l.+( GLOCO-1. )*AMEX**2 ./2 . ) 

C  WRITE( ',*) 'TREE-' , TREE, 'IMAXWR-’ ,TMAXWR, ’ RCT- ’ , KCT, ' PO- ' , PO 

IP( (ABS(TREP-TKAXWR) .LT.S.O) .OR. (RCT.CT.S) )  GOTO  20 
THAXHR-TREP 

C - COMPUTE  TEMP  INTO  THE  BURNER  (TCCIN) - 

TCCIN-TIN* ( PRATl *PT1/P1N)**( (GLOCO-1 . )/GLOCO) 

WRITE) 6, •) 'PRATl-' ,PRATl, '  PT1-' ,PTl 

WRITE(6,*)'TCCIN  -  ’, TCCIN 

CALL  BURNER ( DUMMY , TCCIN , DUMMY , TMAXWR , PIN , PO , IDBl , RODE , 1 ) 

GOTO  10 
20  CONPIG-O 

IDBl-lDBl+1 

WRITE) 6 , * ) ' 1DB1-' , IDBl 
WFITE)9,30) 

WRl ;i) 9, 40)TCCIN,TO,TREr,PIN,PO,PREr,riN,PO,GLOCO, IDBl 
30  PORMAT)5X, 'STATION  IDENTi  WAVE  ROTOR,  STATIONS  3-4’/) 

40  PORMAT) 5X, 'TCCIN-' ,P7 .1 , 2X, 'TO-' ,P7 . 1 , 2X, 'TREE-' , P7 . 1/,5X, 

> ’ PIN- ' , PS . 1 ,2X, 'PO-' , rS.l.ZX, 'PREP- ' , PS. J/, SX, 'PIN-' ,P6. 4,2X, 
>’PO-’ ,P6.4,2X,’GLOCO-' ,P5.3,2X, 'IDBl-' ,Il/) 

RETURN 
END 


SUBROUTINE  BURNER) FIN, TIN, PO, TO, PIN, PO, IDBl , RODE, CONFIG) 

INTEGER  CONFIG 

COMMON/ETAS/ETALC , ETAHC , ETACC , ETALT , ETAAB , VCOEPP 
COMMON/LOSSES/DIPLOS.CCLOSS, ABLOSS ,BPLOSS 
COMMON/BURN, 'HPUEL  ,  PLHV,  TMAXCC  ,  TMAXAB  ,  TMAXWR 
CHECX  COUNTER  OPERATION 

- CONPIG-1  GOTO  WAVE  ROTOR  CALC 

-  CONFIG  RESET  -  0  IN  ROTOR  SUBROUTINE 
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-  IDB  INCREASED  BY  1  IN  ROTOR  8UBR 

- CONFIG-O  AND  IDB-1  GOTO  RAIN  CC  CALC 

- CONriO-0  AND  IDB>1  GOTO  A/B  CALC 

IF ( CONFIG. eO.l)  GOTO  SO 
ID82-IDB2'fl 

NFITE(  6, • ) 'CONFIG-' .CONFIG, '  IDB2-’ ,IDB2 

if(iob2.bq.1)  goto  10 
GOTO  20 

- RAIN  BURNER  CALC - 

10  TO-TNAXCC 
ETAB-ETACC 

DELTAF-(A(T0)-A(TIN)+FIN‘(B<T0)-B(TIH) ) )/(HFOEL+ETAB»FLHV-B(TO) ) 
FO-DELTAF+FIN 
PO-PlN‘f 1 .-CCLOSS) 

HIN-ENTHrN( FIN, TIN) 

HO-ENTHFN( FO.TO) 

PHI1N-ENTFFN( FIN.TIN) 

PHIO-ENTFFN( FO.TO) 

NRtTF(9,60) 

HRITB(9,70)TIN,TO,PIN,PO,FIH,FO 
RETURN 

- AFTERBURNER  CALC - 

- A/B  ON 

20  IF(FODE.EO.O)  GOTO  30 
to-tra::ab 

ETAB-ETAAB 

DELTAr-(A(TO)-A(TlH)4FlN*(BlTO)-B(TlN) ) )/( BFUEL+ETAB* FLBV-B ( TO ) ) 
GOTO  40 

C - A/B  OF' 

30  TO-TIN 
DELTAF-0. 

40  FO-DELTAF+FlN 

PO-PIN*(l.-ABLOSS) 

BIN-FNTKrN(FlN,TlN) 

HO-ENTHFN) FO.TO) 

PHIIN-ENTBFN(FIN,TIN) 

PH10-ENiRFN( FO.TO) 

HF1IE(9,B0) 

KRlTEi9,70)TIN,TO,PIN,PO,FIN,FO 

RETURN 

(  --WAVE  ROTOR  CALC - 

50  TO-TMAXWR 
ETAB-ETACC 

DELTAF-{ A(TO)-A(TIN)+FIN*(B(TO)-B(TIN) ) )/( RFUEL+ETAB»FLBV-B ( TO ) ) 
FO-FIN+DELTAF 

60  FORRAT( 5X, 'STATION  IDENTl  RAIN  BURNER,  STATIONS  3-4'/) 

70  FOFnAT( 5X, 'TI-' ,F5.0,2X, 'TO-' ,F5.0,2X, 'PI-' ,F5.1,2X, 'PO-' ,F5.1, 
>2X,'FI-' ,F4.3,2X,'FO-' ,F4.3/) 

SO  )Om„H  5X, 'STATION  IDENTl  AFTERBURNER,  STATIONS  9-10'/) 

RETURN 

END 


SUBROUTINE  RIXCP(FIN1,FIN2,PIN1,TIN1 ,PIN2,TIN2,FO,PO,TO,VAR,IDH) 
COriMON/ETAS/ETALC ,  ETAHC ,  ETACC ,  ETABT ,  ETALT ,  ETAAB ,  VCOEFF 
COnnON/LOSSES/DI FLOS . CCLOSS .ABLOSS , BPLOS8 
iDH-irn+l 

IF(IDn.EO.l)  GOTO  10 
ZETA-VAR 

GOTO  20 

10  I  ETA- VAR/' (1. -VAR) 
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20  HINl«ENTHFN(riNl,TINl) 

RIN2-ENTHFN( FIN2,TIN2) 

FHl IN1-ENTRFN( FINl ,TIN1 ) 

FHIIN2-EHTRFN(FIN2.T1N2) 

B0-(ZETA*(1.4FIN1  )*HINl4(1.4^FIN2)*BIN2)/(EETA*(1.4^FINl)4(1.4FlN2) ) 
F0-(ZETA»FIN14F1N2)/(1,+1BIA) 

F0-FIN2 

CALL  BFN(FO,BO,TO) 

FB10-ENTBFN(F0,T0) 

IF(IDM.EO.l)  GOTO  30 
GOTO  40 

30  NRITE(9,S0) 

♦IRITE(9,60)T1N1,TIM2,TO,FO,FIN1,FIN2,FO 

RETURN 

40  NBITE(9,70) 

KRITB(9,60)TIN1,TIN2,TO,PO.FIN1,F1N2,FO 
50  F0RMAT( 5X, 'STATION  IDENTj  FRE  LPT  BIXING,  STATIONS  3,5-6’/) 

60  FORMAT ( 5X, 'Tl-' ,F5.0,2X, ’T2-' ,F5.0,2X, 'TO-' , F5 . 0 , 2X , ' FO- ' . F5 . 1 , 

>2X, 'F1-' ,F4 . 3,2X, 'F2-' ,F4. 3,2X, 'FO-' ,F4.3/) 

70  F0RMAT(5X, 'STATION  IDENTi  FRE  NOZZLE  BIXING,  STATIONS  9,10-11'/) 
RETURN 
END 


SUBROUTINE  TURBIN(PIN,TIN,FO,TO,IOI8,FIN,FO,BO,FRAT,FRLC,FRBT,FRLT 
>,XI , BETA, F ASET, IDT) 

COMMON/ETAS/ETALC , ETAHC , ETACC , ETABT, ETALT, BTAAS , VCOEFF 

COMMON/LOSSES/DI  FLOS  ,  CCLOSS ,  ABLOSS  ,  BPLOSS 

COMMON/NORR1/NORKLC,NORRHC 

COMMON/CONST/CJ ,G 

lDT-IDT+1 

BIN-ENTHFN(FIN,TIN) 

FHIIN-ENTRFN(FIN,TIN) 

IFdDT.EQ.ll  GOTO  10 
IF(RASET.EO.l)  GOTO  20 
GOTO  30 

C - HPT  CALCULATION - 

10  HO-HlN-WOFFHC/( (l.-XI)»(l.*FIN) ) 

CALL  HFN(FIN,HO,TO) 

PHIO-ENTRFN(FIN,TO) 

HOI S-HIN-(HIN-HO) /ETABT 
CALL  HFN(FIN.BOIS,T01S) 

PBIOIS-ENTPrN(FIN,T01S) 

FRHT-1./EXP(CJ/BGAS(FIN)»(PHI0IS-PHIIN) ) 

FO-FIN 

PO-PIN/PRHT 

GBAP-GAVG( FO,TIN,roiS,PHIIN,PHIOIS) 

TEMPHTI-ALOGIPO/PIN) 

TEHPHT2-AL0G(T0/T1N) 

EHPT-GBAR/(GBAR-1 . ) •TEBPHT2/TEMPHT1 
KRITE(6,* ) 'TO-' , TO, 'TIN-' ,TIN, 'PO-’ ,PO, 

>’PIN-' ,PIN, 'GBAR-' ,GBAR, 'ETABT-' , ETABT 
WR1TE(6,*)’EHPT-',EHPT 
NBITE(9,40) 

NRITE(9,50)TIN,TO,PIN,FO,FIN,FO,ETABT,GBAR,EHPT 
RETURN 

C - CALCULATE  OPTIHUH  BYPASS  RATIO 

20  FRLT-1 ./( ( 1 .-BPLOSS )/(l. -CCLOSS )*PRLC/PRAT*PRHT) 
FHIOIS-ENTBFN(FIN,TIN)+RGAS(FIH)/CJ*ALOG(l./FRLT) 

CALL  PHIFN(FIN,PHIOIS,TOI8) 

BOIS-ENTBFN( FIN,TOI8) 
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HO«RIN-ETALT*(RIN-tiOI8) 

CALL  lirN(  riN,RO,TO) 

PHIO-ENTRrN(riN,TO) 

■ETA-(l.+riN)*<HIN-RO)/(MORRLC)-l. 

C  NRITC( 6,*) 'TURBIN  BETA- '.BETA 

ro-riH 
PO-PIN/PRLT 

GBAR-GAVG(r0,TlN,T0IS,PRIlN,PR101S) 

TEHPLT1-AL0G(P0/P1N) 

TE«PLT2-AL0G(T0/T1N) 

ELPT-GBAR/(GBAR-1 . )»TBRPLT2/TBHPLT1 
NR1TE( 9,60) 

NrilTE(9,50)TIN,TO,PlN.PO,riN,PO,ETALT,GBAR,ELPT 

RETURN 

c - LPT  CALCULATION  BYPASS  PRESCRIBED - 

30  MO-fllN-( 1 .+BETA) *(WORRLC)/(l .+riN) 

CALL  MENU  IN,  HO,  TO) 

PK10-ENTRrN(riN,T0) 

R01S-H1N+(H0-HIN)/BTALT 
CALL  RrN(riN,ROIS,TOI5) 

PRIOIS-ENTRrN(PIN,TOIS) 

PRLT-I ./BXP(CJ/RGAS(riN)*(PB10IS-PBIIN) ) 

ro-riN 

PO-PlN/PRLT 

GBA(>-GAVG(r0,TlN,T0IS,PBIIH,PH10IS) 

TErPLTl-ALOG(PO/PlH) 

TEnPLT2-ALOG( TO/TIN) 

Ei,PT-GBAB/(GBAR-l.  )»TEHPLT2/TEMPLT1 
WBITE(6,»)'T0-' ,TO, 'TIN-' .TIN, 'PO-' ,PO, 'PIN-' ,PIN, 
>'GBAR-’,GBAR,'  ETALT- '  ,  BTALT 
NBITE(6,*)'ELPI-’,BLPT 
KE2TE(9,60) 

WHTE(  9,50)TIN,TO,PIN,PO,riN,rO,ETALT,CBAR,BLPT 
40  rORNAKSX, 'STATION  IDENTi  BP  TURBINE,  STATIONS  4-5'/) 

50  FORMAT) 5X, 'T1-' ,P5. 0, 2X, ' TO- ' , r5 . 0 , 2X , ' PI- ’ , r5 . 1 , 2X , ' PO- ' , 

>r5. 1  .2X, 'fl-' ,r4.3,2X, 'PO-' ,r4.3/,5X, 'ETA-' ,P6.4,3X, 'OBAR-' , 
>rf.  <1  ',5X,  'ET-'  ,r6.4/) 

60  rOHMAT) 5X, 'STATION  IDENTi  LP  TURBINE,  STATIONS  6-S'/) 
RETURN 
END 


SUBROUTINE  DUCT) PIN ,TIN, E IN , PO, TO, EO, HO) 

COMMON/LOSS  ES/DI ELOS , CCLOSS , ABL0S8 , BPLOSS 
PO-PIN*) 1 .-BPL06S) 

TO-TIN 
EO-0 .0 

BIN.ENTHEN(riN,TIN) 
p;u  it;-tHTp.rN)riN,TiH) 

BO-ENTHEN)EO,TO) 

PHIO-ENTREN)EO,TO) 

HBITE)9,10) 

KRITE)9,20)TIN,TO,PIN,PO,EIN,rO 
10  EORMRT)SX, 'STATION  IDENTi  BYPASS  DUCT,  STATIONS  2-7'/) 

20  rORnAT(RE, 'TI-' ,r5.0,2x, 'TO-' ,E5.0,2X,'PI-' , P5 . I , 2X, ' PO-' , 
>E5.1 ,2X, 'PI-' ,r4. J,2X, 'fO-' ,r4.3/) 

RETURN 

END 


SUBROUTINE  RIXCA) PIN2 . PINl , TIN2 , TINl , RIN2 , BlNl , PIN2 , PINl , BETA, CHIN 
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>i,piiiAx,TO,ro,po,cno) 

CO«MON/CONST/CJ,0 

FUNCTION  DEFINITIONS 

PHI IFNI  GANNA, AHACH)-AMACH«»2. •( 1. ♦ t OANNA-1. )/2.*AMACH«*2. )/| |1.+GA 
>NNA*ANACH*»2. )»»2. > 

PSIlFN(G,GANNA,P,ANACM)-|AnACB*S0HT(O»GAHNA/II)  )/( 1  .-f  (GAHKA-1 .  )/2  . * 
>ANACH**2. )*•( IGANNA+l. )/(2.«(GANHA-l. ) ) ) 

FOR  FAN  8TREAN  MACH  NUNBER 

CHECK  IF  BYPASS  GREATER  TRAN  0.0 
IF  (BETA. GT. 0.0)  THEN 
GOTO  5 

ELSE 

TO-TINl 
P0-PIN1»PINAE 
RO'HINl 
FO-FINl 
CNO-CNINI 
AI2AI1-0.0 

ENDir 
GO  TO  9 

S  GL0CI2-GL0C(F1N2,TIN2,R1N2) 

C  NFITEI 6, * ) 'AMACH8-' ,CH1n1 

GLOCIl-GLOCIFINl .TINl.BINl) 

FLAG1-(PIK2/P1N1*(I .♦(GLOCIl-1. )/2 . *CHIN1* *2 . ) ••(GLOCIl/IGLOCI 
>1-1 .) ) )**( (GLOCI2-1. )/GLOC12) 

C  KP1TE(6, * ) ■ FLAGl-* .FLAGl 

IF  ( BETA. GT. 0.0)  THEN 

IF(FLAGl.LE.l.O)  GOTO  «0 
ENDIF 

CHIN2-SOBT(2./(GLOCI2-1, )«({PIN2/PINl*(l,+(GLOCIl-l. )/2.»CHIN1»«2. 
>)*• I GLOCI 1/(GL0CI1-1. )))»»( (GLOCI2-1. )/3LOCl2)-l. ) ) 
IF(CNIN2.GE.0.e)  GOTO  40 


Ct*tff»t»f*»»»*t****t**t»**»»*S0lVt  FOR  HIKER  INLET  AREA  RATIO 

A12AI1-8ETA/(1.+FIN1)*CH1N1/CH1N2*S0RT( (GLOCIl*ROAS(riN2)«TIN2*(l. 
>+( GLOCI 1-1 . )/2.*CHlMl»»2. ) )/(GLOC12*RGAS(FINl )*TIN1*( 1.+ (GLOCI 2-1 . 
> )/2.*CN1N2*'2.  ) )  ) 

FOR  HIRER  TEHPERATORE  RATIO 


BIN2-RGAS(r]N2) 

RINl>RGAS( FlNl ) 

CPIN2-CP( riN2,TIN2) 

CPINl-CP(EINl,TINl) 

BETAPR-BETA/( 1 .+F1H1 ) 

CP0-(CP1N1+BETAPR*CPIN2)/|1 .+BETAPR) 

TO-TINl* ( CPIN1/CP0*( 1 .+BETAPR*CPIN2*TIN2/(CPIN1*T1N1 ) )/( 1 .+BETAPR) 

>  I 

TAUMIK-TO/TINl 

C*******«**********<*******«***<**sOLVE  FOR  HIRER  ERIT  HACH  NUNBER 

Fo-(B^T.^'r;n2  +  ^INl  )/(1.+beta) 

HO-ENTHFNI rO,TO) 

RO-RGAS( rO) 

CLOCO-GLOC(FO,TO,HO) 

PHIIK2-PHIlrN(GLOCI2,CHIN2)  • 

PHIlNl-PHIlrN(GLOCIl,CHlNl) 

PHIO-( ( ( 1 .+BETAPR)/(SQRT(1./PHI1N1 )+BETAPB*SQRT(1./PHI1N2)*SQRT(RI 
>N2*TIN2*GL0C11/!RIN1*TINI*0L0CI2) ) ) )**2. ) •OLOCI1*RO»TAUH1R/(GLOCO* 
>RIN1 ) 

CHO-SQRT( 2 . *PHIO/( (I .-2.*GLOCO*PHIO)*SORT( 1.-2.*(GL0C0+1 . )*PHI0) ) ) 

C**<***«**<<*><*«<a<*<>*>*><******<***««*«gQi.vE  FOR  HIRER  PRESSURE  RATIO 
POI-PINl*( 1 .♦BETAPR)*(1./(1.4AI2AI1 ) )*8QRT(TA0HIR)*PSI1FN(0,GL0CI1 
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>  ,RIN1  ,CniNl  )/PSIlFN(G,GLOCO,llO,CnO) 

Al-FSIlrN(G,GLOCIl.ItINl,CHIHl) 

A2-F5IirN(G,GLOCO,RO,CIIO> 

FIHIXI-FOI/FlNt 
FiniX-FiniXI*FIKAX 
FO-FiniX*FINi 
9  MRITCO.IO) 

NRITE(9,30)T1N1.T1N2,TO,FIN1,F1N2,FD 

MRITE(9.30)FlNl.FIN2,FO,CniHl,CniN2,CHO,AX2AIi 

RETURN 

10  rORHAKSX, 'STATION  IDENTt  CONST.  AREA  HXXER,  8TATXON8  7,8-9’/) 

20  rORnAT(5X, 'T1-' ,r5.0,2X, 'T2-' ,r5.0,2X,'TO-' ,F5.0,2X, 'F1-' ,F5.1,2X, 
>’F2-’ ,r5.1,2X,'FO-' ,F5.1/) 

30  rORKAT(SX, 'rl-',r4.3,2X, 'F2-',F<.3,2X, 'rO-',F4,3,2x,'Rl-',F4.2,2X, 
>'H2-’ ,F4.2.2X, 'NO-' .F4 .2 . 2X, * ARAT> ’ , P4 . 2/) 

40  NRITC(6,S0) 

50  FORnAT( SX, 'FAN  STREAM  MACH  NO  XS  GREATER  TBAN/BOUAL  TO  0.8’) 

RETURN 

60  NRXTE(6,70) 

70  rCRHATf 5X, 'FAN  STREAM  MACS  NO  XS  LESS  TRAN/BQUAL  TO  0.0’) 

RETURN 

END 


8URROUTINE  N0ZZLE(PXN.TXN.F0,T0,T0X8,PXN,PRLC,VEXIT,BXHACB,FS0, 
>irLAG) 

COMMON/ETAS/ETALC ,  ETAHC ,  ETAwC ,  ETAHT ,  ETALT ,  ETAAB ,  VCOEFF 
COHMON/LOSSES/DIFLOS.CCLOSS.ABLOSS.BPLOSS 
COMMON/CONST/C J , G 

PRN-1 , /( ( 1 . -DX  FLOS ) * ( 1 . -BPL08S ) * ( 1 . -ABLOSS ) * ( FRLC ) ) 

PRN-PSO/PIN 

HIN-ENTHrN(FIN,TXN) 

PRI1N>ENTRFN(FXN,TXN) 

C10CO»GI.OC(FXN,rXN,»IN) 

FRN-PSO/PXN 

XF(FRN.GB.l.O)  GOTO  30 
wrlte(6,«|'GLOCO-',GLOCO 
wtlte(6,«)’FSO-',PSO 
wtlte(6,*)’PBN-’,FRN 
PHI0XS-PHIIN+BGAS(FIN)/CJ*AL0G<FRN) 

CALL  FHIFN(FIN,FHX01S,TOXS) 

H0IS-ENTHFH(FXN,T0XB) 

FO-FIN*PRN 

ETAN-VCOErF*»2. 

HO-HIN-ETAN*(HIN-HOrS) 

CALL  HrN(FXN,HO,TO) 

VEXIT-SQRT( 2 . *G*CJ*(HlN-HO) ) 

AEXIT-SQFT(G*RGAS(FXN)*TO/(l.-RGAS(FIM)/CJ»rO/BO) ) 

BXHACH-VEX  X  T/AEXXT 
NRITE(9,10) 

NRXTE(9,20)TIN,TO,FXN,FO,FXN,VEXIT,EXHACB 
NRXTE(9,25)0L0C0,FRN,AEXXT 
GOTO  40 

10  FORnAT( SX, 'STATION  XDENTt  EXRAUST  NOZZLE,  STATIONS  11-12’/) 

20  FORHAT(5X, 'TI-' ,F5.0,2X, 'TO-' ,F5.0,2X, ’FI-’ ,F5.1,2X, ’FO-' ,F5.1, 
>2X, 'FI-',F4.3,2X, 'VO-’,F5.0,2X, •MO-’,F4.2) 

25  FORHAT(5X, 'GLOCO-' ,F6.4,2X,'PRN-’ ,F6.5,2X, ' AEXXT-' , F6 . 1/) 

30  XFLAG-l 


64 


no  n  n  o o 


40  RETURN 
END 


SUBROUTINE  FERE (AHACRB . PIN, VO. VEEIT, BETA, BPTRR, SEC, BXNRCR ) 

COMMON/CONST/C J , G 

SFrHR-( (1 .+EIN) •VEXIT-V0)/0 

SEC-EIN»3600./BFTHR 

KBITE(9,10) 

NR I TE ( 9 , 2  0 ) ARACRB , 8PTHR , SEC, BETA , EXHACR 
NRITEt9,25) 

10  EORHAT(5X, 'PEREORRANCE  PARAHETERS ' /) 

20  E0RMAT(5X,'AMACR8-’,E6.3/,5X,'BP.THRUST  ,E6.2,2X, 'SEC  -  ’,E5.3 
>  ,  4X , • BETA-  •  ,  E6 . 3 , 2X, • EXMACH- ’ , E6 . 3  ) 

2 5  EORMAT ( 5X , • ( LB E/LBM/S ) ' , BX . ' <  LBH/RR/LBE ) ' / ) 

RETURN 

END 


SUBROUTINE  HEN(ERATIO,Hex.TEHPEX> 

DUMMY-1 . 

ITERNO-0 

I HEX- 1  El X( HEX/1 00. )»1 00 
HONE-ELOAT( IHEX) 

HTWO-HONE+100 . 

CALL  TEMP ( HONE , HTWO , DUHHY ,DUHM* , TEMPI , TEMP2 ) 

10  ITEfnO-ITEBNO+1 

IE( ITERNO.GT.50)  GOTO  20 
HONE-ENTHEH!  EP.ATIO,  TEMPI  ) 

HTWO-ENTHFH(FRATIO,TEMP2) 
BLOPE-(HTWO-HONE)/(TEMP2-TEMP1) 
TEMPl-TEMPl+(HEX-HONE)/SLOPE 
H0NE-BNTHEN( ERATIO, TEMPI ) 

IE(ABS(HONE-HEX) .GT. 0.002)  GOTO  10 
TEMPEX-TEHPl 
GOTO  30 
20  NBITE(6,40) 

40  EORHAT( 5X, 'SECANT  HETHOO  ITERATION  DOES  NOT  CONVERGE'/) 
30  RETURN 

END 


SUBROUTINE  PHI rN< ERATIO, PHIEX.TBMPBX) 

DUMMY-1 . 

ITEPNO-0 

IPHIEX-IEIX(PHIEX*10.  ) 

PH  I  ONE- FLOAT ( I  PH  I  EX) *0.1 
PHITWO-PHIONE+0.1 

CALL  TEMP  I  DUMMY, DUMMY, PHIONE, PHI TWO, TEMPI, TEMPI) 
10  ITEFNO-lTERNO+1 

IF( ITERNO.GT.55)  GOTO  20 
PHIONE-ENTFFNIERATIO, TEMPI ) 

PHITWO-ENTRFN(  ERATIO, TEMF2  ) 

SLOPE-) PHITWO-PHIONE)/(TEHP2-TEMPl) 

TEMP 1-TEHP1*(PH1EX-PR TONE) /SLOPE 
PHIONE-ENTFFfJ)  ERATIO, tempi  ) 

IF(ABS( PHIONC-PBIEX) .CT. 0.000003)  GOTO  10 
TEMFEX-TEHFl 
GOTO  40 
20  NRITC(6,35) 
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35  rORHAT(5X, 'SECANT  HETBOD  ITERATION  DOES  NOT  CONVERGE'/) 
40  RETURN 
END 


FUNCTION  STATEMENTS 


FUNCTION  BNTnFN(FRATIO,TENF) 

IFCTEMP.LT.O.O)  GOTO  10 
Cl-0. 24062 
C2-0.017724E-03 
C3-0.038056E-06 
C4-0.012662E-09 
C5-0.0013012E-12 
Dl-0. 22091 
D2>0.5ie22E-03 
O3-0.19462E-06 
D4-0 .04S089E-09 
05-0. 0043275E-12 

ENTHrN-(l./(l.+FRATIO))»(Cl»TEMP-C2/2.*TEMP*‘2.+C3/3.*TEMP**3.-C4/ 
>4  .  «TEMP**4  .+05/5  .  »TEMP*»5.  )♦(  FRAT10/(  1  .♦FRATIO)  )  ‘  ( D1  •TEMP4-D2/2  . ‘TE 
>MP**2.-D1/3.»TEMP»»3.+D4/4.»TB«P»*4.-D5/5.*TEHP*«5. ) 

GOTO  25 

10  VRITEI 6,20)TEMP 

20  FORMAT! SX, 'NEGATIVE  TEMP  IN  ENTHFB  ENCOUNTERED- F6 . 1/) 

25  CONTINUE 
RETURN 
END 

FUNCTION  ENTRrN(FRAT10,TEHP) 

IF(TEMP.LE.O.O)  GOTO  10 
Cl-0. 24062 
C2-0.017724E-03 
C3-0.03e056E-06 
C4-0.012662e-09 
C5-0.0013012E-12 
Dl-0. 22091 
D2-O.51022E-O3 
D3-0.19462E-06 
D4-0.045089E-09 
D5-0.0043275E-12 

ENTRFN-!  1  ./(I  .-fFRATlO)  )*(Cl*ALOG(rEnP)-C2*rEMP+C3*0. 5‘TEMP»*2.-C4/ 
>3.*TEMP*»3 .♦C5/4.*TEHP**4. )+(FRATIO/( l.+FRATIO) ) ‘ ( Dl ‘ALOG! TEMP ) +02 
>‘TEMP-D3*0.5»TEMP»»2.+D4/3.»TEMP**3.-D5/4.»TEHP**4. ) 

RETURN 

10  NRITE(6,20)TEMP 

20  FORMAT! 5X, 'NEGATIVE/IERO  TEMP  IN  BKTRFN  ENCOUNTERED  -’,F6.1) 

RETURN 

END 

FUNCTION  A!T) 

C1-. 24062 

C2-1. 77246-5 

C3-3.8056E-8 

C4-1.2662E-11 

C5-1.3012E-15 

D1-. 22091 

D2-5.1822E-4 

D3-1.9462E-7 

O4-4.5089E-11 

D5-4.3275E-15 
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A-C1*T-C2/2.»T«»2.+C3/3.*T**3.-C4/4.»T»*4.^C5/5.*T**5. 

RETURN 

END 

C 

rUNCTION  B(T) 

C1-. 24062 
C2-1 .7724E-5 
C3-3.8C56E-8 
C4-1 .2662E-11 
C5-1.3012E-15 
D1-. 22091 
D2-5.ie:2E-4 
D3-1 .9462E-7 
D4-4.b089E-ll 
D5-4.3275E-15 

B-Dl*T+D2/2.*T»»2.-D3/3.»T‘*3.+t)4/4.«T**4.-D5/5.«T«*5. 

RETURN 

END 

C 

rUNCTlON  CPIfRATlO.TENP) 

Cl-0. 24062 
C2-0.017724E-03 
C3-0  C3S056E-06 
C4-0 . 012662E-09 
C5-0.0013012E-12 
n-r  22091 
D2-0.;i822E-03 
D3-0.19462E-06 
04-'’  ,046089E-09 
O5-0.0043275E-12 

CP-( 1 ./( 1 .4rRAriO) ) *(C1-C2*TEMP+C3‘TENP**2.-C4‘TEHP**3.+C5»TEHP‘*4 
> . l+(rRATIO/(1.4PRATlO) )»(Dl+D2*TEnP-D3»TEnP**2.+D4*TEMP**3.-D5*TEN 
>P“4.  ) 

BEa 

END 

C 

rUSCriON  RGRS(rRAT) 

W'.C'I  .( 1  ♦rPATl/(0. 034522+0. 035648*rRAT) 

RG;.'!  -1545 . 43/WHOL 
RETURN 
END 
C 

rUNCTION  GAVG( ERAT.TIN.TOIS, PHIIN.PHIOIS ) 

COflAON/CONST/CJ .  G 

c  HPIIEl 6, 5)rRAT,TJN,TOIS,PHIIN,PHI01S,CJ 

c  5  rOF«AT(5x,r3.1,5x,2rB.l,5x,2r7.2,5A,r6.2) 

GAVG-l ./(I .-(BGAS(rRAT)/Ca)»(ALOGlTOIS/TlM)/(PHIOIS-PHlIN) ) 1 
P"’T'-rn 
ENO 
C 

FUNCTION  GLOC( r LOG, TLOC, BLOC) 

COHAON/CONST/CJfG 

GLOr-1 ./( 1 .-RGAS(rLOC)/CJ»TLOC/HLOC) 

RETURN 

END 


eUBSOUTINE  TEMP ( HONE, HTWO, PHIONE, PH ITNO, TEMPI ,TENP2) 
jri (HONE.EO.0,00) .OB. (PHIONE.LE. 1 . 3) )  GOTO  10 
IF' (RONE.E0.100. 1 .OR. (rHlONB.EQ.1.4} )  GOTO  20 
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ir< (RONE. eQ.200. ) .OR. (PRIONB.eO.l -5) )  GOTO  30 
IP((RONE.ED.300.).CR.(PR10NE.ED.1.6))  GOTO  40 
ir(  (HONE.EQ.  400.  )  .OIt.(PRIONB.E0.1 .7)  )  GOTO  50 
IP((RONB.EQ.500.).OR.(PR10N£.B0.1.B))  GOTO  60 
ir((HONE.GE.600.).OR.(PR10HB.GB.1.9))  GOTO  70 
NRITE(  6  ,  DRONE ,  PHIONE 

1  rORnAT(SX, 'NOTHING  HRTCRBDtBONB-' , P13 . 9 , 5X, ' PRIONE-' , P13 . 9/) 
STOP 

10  TEMPl-0. 

TEMP2-450. 

RETURN 

20  TEMPl-300. 

TEMP2-eOO. 

RETURN 

30  TERPl-500. 

TEMP2-1250. 

RETURN 

40  TEMPl-750. 

TEHP2-1600. 

RETURN 

50  TBHPl-1100. 

TEHP2-2000. 

RETURN 

60  TEMPl-1500. 

TEHP2-2400. 

RETURN 

70  TEHPl-2200. 

TEMF2-3300. 

RETURN 

END 


PLOT  SUBROUTINES  USING  'OISSPLA'  SOFTWARE  -  PACILITT  DEPENDENT 


SUBROUTINE  PLOT ( XARAY , YARAY , lARA Y , BPR , J J , I 1 ) 

DIKENSION  XARAYdl)  .YARAYdl)  ,ZARAY(I1)  ,BPR(11) 

REAL  XARAY, YARAY, ZARAY, BPR 

INTEGER  11, JJ 

CALL  PHYSOR(0.5,0.5> 

CALL  AREA2DdO.  ,7,5) 

CALL  PRAME 

CALL  GRAr( 0. , 'SCALE' , 140 ., 0 .,’ SCALE' , 100. ) 

CALL  GRAF ( 50, , 'SCALE' ,120. ,0.6, 'SCALE' ,1.6) 

CALL  XNAMECSP. THRUST' ,9) 

CALL  YNA«E( 'TSFC' ,4 ) 

GOTOd0,20,30,40,50,10,20.30,40,50)  ,JJ 


10 

CALL 

SETCLR( 'YELLOW' ) 

CALL 

MARKER! IS) 

GOTO 

60 

20 

CALL 

EETCLRI 'CYAN' ) 

CALL 

HARKERdS) 

GOTO 

60 

30 

CALL 

SETCLRl 'MAGENTA' ) 

CALL 

HARKER(5) 

GOTO 

60 

40 

CALL 

6ETCLR( 'GREEN' ) 

CALL 

MARKER! 9) 

GOTO 

60 

50 

CALL 

SETCLR! 'RED' ) 

CALL 

MARKER! 2 ) 

60 

CALL 

CURVE ! XARAY , YARAY , I 1 . + 1 ) 

68 


uuuuuu  ouuuuuu  uuu  u  o  u  uuu 


CALL  ENDGR(O) 

RETURN 

END 


PLOT  ROUTINE  USING  GRAPRIT  SOFTWARE 


SUBROUTINE  GRSPLOTl(](lA,ylA,y2A,  IAN4  ,  ICPl ) 

PLOT  or  SPECIFIC  THRUST  AND  SFC  VS  INPUT  PARAMETER 
CHOICES  OF  INPUT  PARAMETERS  AREl 

1.  ByPASS  RATIO 

2.  PAN  PRESSURE  RATIO 

3.  COMPRESSOR  PRESSURE  RATIO 

4.  LPT  EXIT  MACS  NUMBER 

REAL  XlA(50),yiA(50),y2A(50> 

CKARACTEi'  XlLABEL*6  ,  ylLABEL*33  ,  y2LABCL*  33 
CHARACTER  FIG1*7 ,NUMBl«3 ,TOPLBL*10 


SPECIPy  X,y  LEFT  AND  T  RIGHT  LABELS - 

yiLA3CL-’At  Specific  Thruit  Ibf /( Ibsi/e ) ' 
y2LA8EL-'Bi  SFC  (Ibe  fue3/br>/lbf  thru»t’ 

CHOOSE  THE  X  PARAMETER  TO  PLOT 
IP  (IAN4.E0.1)  THEN 
Xt LABEL-'  BETA  ' 

ELSriF  (IAN4.EQ.2)  THEN 
X1LABEL-’  PRLC  ’ 

ELSFIF  (1AN4.EQ.3>  THEN 
X1LABEL-'  PFHC  ' 

ELSEir  (IAN4.E0.4)  THEN 
XI LABEL-'AMACHS' 

END  IF 

WBlTE(6,‘)'lNe”Ti  ' 

HRlTEte,*)'  0  -  AUTOMATICALLV  DETERMINE  XIMIN  AND  XIMAX' 
WF,Hi,6,*)'  1  -  INPUT  XIMIN  AND  XlKAX' 

REAO(5,*)IANSS 
IF  'IANS5.EQ.1)  THEN 

«RITE(6,*)'  INPUT  XIMIN  -l' 

PEAD(5,*)X1MIN 
VPITEie,*)'  INPUT  XIMAX  -l' 

READ*  5,« )X1HAX 
ELSE 

AUTCMATICALLy  DETERMINE  THE  MIN  AND  MAX  X  VALUE 

X1MIN-X1A( 1  I 

X1MAX-X1A{1) 

DO  J-l,ICPl 
X1TEMP-X1A( J) 

X2IE«F-XIA( J) 

IF  (TlMIN.GT.XlTEMP)  THEN 
XIMIN-XITEMP 
ENOir 

IF  (XIMAX. LT.X2TEHP)  THEN 
X1HAX-X2TEHF 
ENDIF 
END  DO 
ENDIF 

MRITE(6,* ) 'XIMIN-' .XIMIN, 'XlMAX-' , XIMAX 
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c - INPUT  THB  NUKBER  OP  THE  PIGURE  TO  PLOT  (PIG  1> - 

NRITE(6,*)'INPUT  THE  PIGURE  NUMBER  (!•  1  OR  2,..)' 
READ(S, ’ (A9)' )NunBl 
C  NRITE(6,»)’NUMBl-',NUHBl 

PIGl-'PIGURE  ' 

TOPLBL-PIG1//NUHB1 

C  NRITE(6, '(lX.All J'TOP  LABEL  .TOPLBL 

C—.—.— 

c - CHECR  THE  DATA  INPUT 

C  DO  11-1, ICPl 

C  MRITE(6,»)Il,XlA(Il),YlAtIl),T2A(Il) 

C  END  DO 

C - CONSTRUCT  THE  PLOT 

C 

NRlTE(6,*)‘STANDBy  POR  PLOT' 

C - OPEN  GRS,  OPEN  AND  ACTIVATE  NOP.RSTATIOH. 

CALL  OPNGRS< 1,2,1) 

C 

C  BET  UP  THE  COLOR  TABLE 


CALL  GSCRd, 0,0. 00, 0.00, 0.50) 
CALL  GSCRd, 1,1. 00, 1.00, 0.00) 
CALL  GSCRd, 2, 0.00, 1.00, 0.50) 
CALL  GSCRd, 3, 0.00, 1.00, 1.00) 
CALL  GSCRd, 4, 1.00, 0.00, 0.00) 
CALL  GSCRd, 5, 1.00, 0.00,1. 00) 
CALL  GSCRd, 6, 1.00, 0,00, 1.00) 
CALL  GSCRd,7,1.00,l,00,1.00) 
CALL  GSCRd, a, 0.00, 0.00, 0.00) 


C - SET  VIEWPOINT  AND  DEPINE  THE  RANGE  OP  USER’S  DATA  TO  BE  USED 

CALL  AGSETP  ( 'GRID/LEPT. ' ,  .15) 

CALL  AGSETP  ( 'GRID/RIGHT. ' ,  .B5) 

CALL  AGSETP  ( ' GRID/TOP . ’ ,  .85) 

CALL  AGSETP  ('GRID/BOTTOM.*,  .15) 

C - SET  UP  POR  FIRST  PLOT 

C -  DEFINE  THE  TOP,  BOTTOM,  LEFT  AND  RIGHT  LABELS 

C -  PUT  BLANRS  IN  FIRST  CURVE  LABELS 


CALL  AGSETC  { ' LABEL/NAHE . • ,  'T') 

CALL  AGSETI  (’LINE/NUMBER.',  100) 

CALL  AGSETC  ('LINE/TEXT.',  '  ') 

C 

CALL  AGSETC  (’LABEL/NAME.',  'B') 

CALL  AGSETI  ( ' LINE/NUHBER. ' ,-100 ) 

CALL  AGSETC  ( ' LINE/TEXT . ' ,  '  ') 

C 

CALL  AGSETC  ('LABEL/NAHE.',  'L') 

CALL  AGSETI  ( 'LINE/NUMBER. ' ,  100) 

CALL  AGSETC  ('LINE/TEXT.',  '  ') 

CALL  AGSETC  ('LABEL/NAHE.',  'R') 

CALL  AGSETI  (' LINE/NUMBER .' ,-100 ) 

CALL  AGSETC  ('LINE/TEXT.',  '  ') 

C - turn  OPP  RIGHT  AXIS 

CALL  AGSETP  ( 'AX I S/RIGHT/CONTROL .', 0 . 0 ) 

C - SET  THIN  AND  TMAX  POR  FIRST  PLOT 

CALL  AGSETP  ( 'X/HINIHUH. ' ,  XIHIN) 

CALL  AGSETP  ( 'X/HAXIHUM. ' ,  XIMAX) 

CALL  AGSETP  ( ' Y/HINIMUM . ' ,  0.0) 

CALL  AGSETP  ( ' Y/HAXIHUH. ' ,  120.0) 

C - 


on  n  onooo  on  noonn 


CALL  AGSTUP  (XIA, 1 , 1 , ICPl . 1 , Y lA, 1 , 1 , ICFl , 1 ) 

CALL  AGBACK 
CALL  GSPLCI  (7) 

CALL  DASHDC  3 , 1  > 

C - plot  curve  If  XIA  VS  ylA 

CALL  AGCURV  ( XIA, 1 , VIA, 1 , ICPl , -1 ) 


- SET  UP  roR  second  plot 

-  DEEINE  THE  TOP,  BOTTOM,  LEFT  AND  RIGHT  LABELS 

-  THE  LARGER  LINE  NUMBER  WRITES  ABOVE  PREVIOUS  LABELS 

CALL  AGSETC  ( ' LABEL/NAME . ' ,  'T') 

CALL  AGSETI  ( ' LINE/NUMBER . ' ,  100) 

CALL  AGSETC  ( ' LINE/TEXT . ' ,  TOPLBL) 

CALL  AGSETC  <’ LABEL/NAME ,  'B') 

CALL  AGSETI  {' LINE/NUMBER -1 10 ) 

CALL  AGSETC  (' LINE/TEXT XILABEL ) 

CALL  AGSETC  ('LABEL/NAME.’,  'L') 

CALL  AGSETI  ('LINE/NUMBER.',  110) 

CALL  AGSETC  (' LINE/TEXT. ', YILABEL ) 

CALL  AGSETC  ('LABEL/NAME.',  'R') 

CALL  AGSETI  | 'LINE/NUMBER -1 00 ) 

CALL  AGSETC  t  '  LINE/lEXI  .  '  ,  VILAbEL 
TURN  orr  THE  TOP,  BOTTOM  AND  LEFT  AXES 
CALL  AGSETF  ( • AXIS/TOP/CONTROL . * , 0 . 0 ) 

CALL  AGSETF  (’ AXI S/BOTTOM/CONTROL 0 . 0 ) 

CALL  AGSETF  < 'AXI S/LEFT/CONTROL 0 . 0 ) 

- TURN  ON  THE  RIGHT  AXIS 

CALL  AGSETF  ( ' AXIS/RIGHT/CONTROL . ' , 1 . 0 ) 

- CHANGE  the  YMIN  AND  THAX 

CALL  AGSETF  ( ' T/MINIMUM . ' ,  0.0) 

CALL  AGSETF  < ' T/MAXIMUM . ' ,  1.5) 

- TORN  ON  THE  RIGHT  AXIS  NUMERICS 

CALL  AGSETF  ( ' RIGHT/NUHERIC/TYPE . ’ , 1 . E6 ) 

CALL  AGSTUP  (XIA, 1 , 1 , ICPl , 1 , TZA, 1 , 1 , ICPl , 1 ) 

CALL  AGBACX 
CALL  GSPLCI  (4) 

CALL  DASHDC  ( ’ $$SFCS ’  "  "  , B , 1 ) 

- PLOT  CURVE  2l  XIA  VS  Y2A 

CALL  AGCURV  ( XIA, 1 , Y2A , 1 , ICPl , -2 ) 

- CALL  FRAME  TO  ADVANCE  THE  FRAME 

- CLOSE  GES,  CLOSE  AND  DEACTIVATE  WORRSTATION 

CALL  FRAME 
CALL  CLSGRS 

C—.... 

END 


71 


APPENDIX  B 


PROGRAM  WRCOMP 


A.  PROCEDURES 

1 .  WRCOMP 

a .  Log  on  to  the  VAX 

b.  At  the  DCL  prompt  type: 

$run  WRCOMPIA . runs  WRCOMP  program 

. select  options  from  screen 

2.  INTERFACING  WITH  ENGINE 

a.  Run  ENGINE  program  in  either  configuration. 

b.  $type  ENGINE1.DAT  or 

$print  ENGINEl.DAT/que=la210_l 

c.  From  the  ENGINE1.DAT  determine: 

1)  Pressure  out  of  HPC 

2)  Temperature  out  of  the  HPC 

3 )  Gamma  at  the  HPC 

4)  Fuel-air  ratio  (f)  from  the  burner 

d.  Run  WRCOMPIA  using  these  values  as  inputs  and  make 
a  guess  for  RREF. 

e.  If  min  doesn't  equal  mout  from  VJRCOMP,  then  change 
the  guess  for  RREF  and  iterate  until  they  are 
equal . 

f.  Run  ENGINEIA  in  the  wave  rotor  configuration  with 
the  correct  value  for  RREF. 
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B.  SAI'IPLE  RESULTS 


riLE  WRPERrl.OUT  PROM  KRCOMPlA.POR  PROGRAM: 

t  or  TIME  STEPS  K-  2000 

CFLNUH-  0.6000000 

I  or  BIEMANN  ITER.  (QSTOP)-  20 
GRID  CELL  WIDTH  DX-  0.01 
PTO  -  79.20000  pel 

PTOTIN-  546064.8  Pa 

RTIN  -  0.1956737  lbra/ft'3 

RTOTIN  -  3.134375  kg/m‘3 

PSEXITl  75.54148  psl 

PSEXIT  -  520840.2  Pa 

PRErl  -  188.8537  (psl) 

PREF  -  1302101.  (Pa) 

xre£  -  0.1600000  a 

xte£l  -  0.5903512  ft 

TIMING  AND  MASS  PLOW  BALANCE! 

INLET  PORT  OPFNS  AT: 

N-  741  TTO'^AL-  1.2576689 

EXHAl'ST  PORT  CLOSES  AT: 

N-1014  TTOTAL-  1.7417805 

INLET  PORT  CLOSES  AT: 

N-1681  TTCIAL-  3.1865280 

P(195)-0. 68590'  P(201)-0.684240 

Min-  0.561596  Mout-  0.626005 


riN  -0.0165  AMEXIT-0.9000 

RTIN-0.1957  RTOTIN-3.1344 

TIE  -  1964.0  TE  -  1707.4 


COMPUTER  RUN  TIME  -  835.5391  sees. 


WAVE  ROTOR  PERPORMANCE  MAP  ANALYSIS: 

SPR-  0.40  T''  -  1.60  TTR  -1.135C 

PTO-  79.20  TOTIN-  1056.00  GAHHA-1.371C 
PTE-126.72  PE  -  75.54  PREP  -ISo.BS 

TREF-  2187.9  RGUESS-O, 22520  RREP-  3.6073 


TIME-  0.0003768 

time-  0.0005218 

TIME-  0.0009547 

P(  203  )-0. 684240 
SWL-  1  SWR-  1 
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ononnonn 


C.  WRCOMP  PROGRAM  LISTING 


Sections  which  are  commented  out  with  a  "C"  contain 
features  that  have  been  disabled  (eg.  DISSPLA)  or  they  contain 
amplifying  remarks. 


c 


c 


c 


302  1 


3022 
C - 


c 


PROGRAM  WRCOMP  WITH  VAN  DER  CORPUT  SAMPLING  AND  SINGLE  TIME  STEP 
INTEGER  QPRINT,  QSTOP , SWL , SWR 
REAL  MASSIN.MASOUT 

DIMENSION  WNORM<12) ,ID1GT{12) ,XARRAY(100) 

DIMENSION  P(203),R(203),0(203),A(203),S(203),X(203) 

COMMON/SUB S/P, R,U, A, S,X 

COHHON/GLIMMl/PGLIM,RGLIM,UGLlM,PL,RL,UL,PR,RR,OR,AL,AR,GL,GR -EPS 
COMMON/GLIMM2/DT,DX,XI 
COHMON/rUNl/G, PA,RA,UA,RB,RMU 
COMHON/SAMPLE/WNORM, IDIGT 
COMMON  XARRAY.Nl 
CALL  COMPRS 
CALL  BLOWUP (1.5) 

CALL  PAGE( 11 . 0,0. 5 ) 

CALL  HWSCAL( ' SCREEN' ) 

DATA  K,SWL,oWR/2000,l,2/ 

DATA  N,CFLNUM,TTOTAL/0, 0.60, 0.0/ 

DATA  PTOTIN,RTOTlN,PSEX it/547292. 900, 3.250,497390.1/ 

DATA  PREF, RREr,XREE/1243475 .4,3.60,0.1800/ 

DATA  SWL,SWR/1,2/ 


OPEN (UNIT- 20, FILE-'WRPERFl .OUT' , STATUS- ' NEW ' ) 

WRITE( 6 , * ) ' INPUT  NUMBER  OF  TIME  STEPS  (2000)!' 

P.EAD(  5  ,  *  )K 

WR1TE( 20 , » ) ' t ILE  WRPERFl.OUT  FROM  WRCOHPlA.FOR  PROGRAM:' 
WRIT£( 20 , •  ) 

WRITE(20,*)'  f  OF  TIME  STEPS  K-  ',K 

WRITE( 6 , ♦ ) ' INPUT  CFL  NUMBER:  (RANGE  0.2  <  CFLNUM  <0.8)' 
READ( 5, ♦ )CFLNUH 

WRITE(20,»)'  CFLNUM-  ', CFLNUM 


WBITF( 6 , * ) ' INPUT  MAX  NUMBER  OF  RIEHANN  ITERATIONS  (20) 
READ( 5, ‘ )OSTOP 
WFITE( 20,3021 )OSTOP 

FOBHAT(5X,'l  OF  RIEMANN  ITER.  ( OSTOP ) -  ',13) 

WRITE( 6 , * ) ' INPUT  GRID  CELL  WIDTH  (0.01):' 

READ( 5 , * )DX 
WRITE( 20, 3022 )DX 

FORMAT( 5X, 'GRID  CELL  WIDTH  DX-  ' , F4 . 2 ) 

--INPUT  TOTAL  CONDITIONS  INTO  THE  WAVE  ROTOR - 

WRITE( 6 , * ) ' INPUT  TOTAL  PRESSURE  IN  (PTO)  (pel):' 

READ( 5 , * )PTO 

FTOTIN-PTO‘144 . *47 .88026 
WRITE(6,»)'  PTO  -',PTO,'psi' 

WRITE(6,*)'  PTOTIN-' ,PIOTIN, 'Pa' 

WRITE(20,*)'  PTO  -',PTO,'psi' 

WRITE( 20, * ) '  PTOTIN-' ,PTOTIN, ' PB' 

WPITE( 6 , * ) ' INPUT  TOTAL  TEMP  IN  (TTOTIN)  (R):' 

PEAD( 5, • )TTCTIN 

WRITE( 6 , • ) ’ INPUT  GAMMA  IN:' 

READ( 5, * )GAMMA 
G-GAMMA 

WRITE( 6 , « ) ' INPUT  FUEL-AIR  RATIO  IN  (FIN):' 

READ(5,*)riN 

MG  I N-(l.+FIN)/(0. 034522  +  0. 03564e*FIN) 
RGIN-1545.43/MGIN 
- density  In  Ibm/ff  3 


RTIN-PTO*144 . 0/(BGlR*tTbtlfr) 

-density  In  kq/n"3 

RTOT  iN-RTi  N*  0 . 4  5  3?9*^TT77irr30^  BT'*^rr 


C 


n  n  o  o 


WBITE(6,*)'  RTIN  -  ',RTIN,’  lbBi/£t'3' 

WR1TE(6,»)’  RTOTIN  -  '.RTOTIN,'  kg/IB*3' 

KRITE(20,*)'  RTIN  »  '.RTIN,'  lb«l/ft'3' 

WRITE<20,*1'  RTOTIN  -  RTOTIN,*  kg/«*3’ 

WRITE! 6, •)' INPUT  TOTAL  PRESSURE  RATIO  (TPR)!* 
READ(5,»)TPR 

WR1TE(6,« 1 ' INPUT  EXIT  HACH  (ABSOLUTE)*' 

READ! 5, * )AMEXIT 
-COMPUTE  EXIT  STATIC  PRESSURE 

PSEXITl-PTO*TPR/(l+(G-l. )/2 .*AMEXIT**2 . )**!G/(G-1 . ) ) 
PSEXIT-PSEXITl *144 .0*47.88026 


PSEXITl  -  ' .PSEXITl , ’pel’ 
PSEXIT  -  * .PSEXIT, ' Pa' 
PSEXITl  -  * .PSEXITl , 'psi * 
PSEXIT  -  * .PSEXIT, 'Pa' 


WP1TE(6, *) ' 

WPITE(6 , » ) ' 

WRITE! 20, * ) ' 

WRITE! 20, * ) ' 

-INPUT  REF  DATA - 

WRllE! 6 , * ) ' INPUT  STATIC  PRESSURE  RATIO  (SPR)i 
READ(5,*)SPR 

-COMPUTE  PREF - 

PREFl-PSEXITl/SPR 
PREF-PREF1*144 .0*47 .88026 
«RITE(6,*)'  PREFl  -  ' . PREFl , '  !psi)' 
WBITE(6,*)'  PREF  -  '.PREF,'  ! Pa )  ' 
WR1TE!20.*)'  PREFl  -  ', PREFl,’  (psl)' 

WPITE!20,*)‘  PREF  -  '.PREF,'  (Fa)’ 

WRITE! 6 ,*)’ INPUT  GUESS  FOR  RREF  !lbm/ft'3)!’ 
READ! 5, * )RGUESS 

RREF-ROUESS*0.45359*!1./0.3048)**3. 
WRITE(6,*)'  RGUESS  -  ’.RGUESS,’  lb»/ff3’ 

WPITE(C,*)’  RREF  -  ’,RREF,’  )tg/lil'3’ 

-INPUT  XREF  IN  METERS - 

WBITC!6 ,*)' INPUT  XREFi  (0.1800  m)’ 

READ! 5 , * )XREr 

-CONVERT  XREF  TO  FEET  FOR  CORRECT  DIMENSIONS  • 
XREFl-XREF/O . 3048 

WRITC(C,*)’  xref  -  ’.XREF,’  »' 

WRIlEii,*)'  xrefX  -  ’, XREFI,’  ft’ 

WPtTE(20.*)’  xtaf  -  ’.XREF,’  m’ 

WFI7E{7),*)’  xtefl  -  '.XREFI,’  ft’ 


G-1 . 35 
GL-1  .  35 
GR-l  .  35 

GL-GAMMA 
GP-GAHMA 
EPS-1 . E-06 
C  O5TOP-70 

Hl-0 

JCOUNT-0 
ECOUNT-O 
UEXHAX-O . 

MASSIN-0 , 

MASOUT-0 . 

C - START  COMPUTER  TIME  CLOCK - 

RINTT-0 . 0 

RCTIME.>SECNDS!RINTT) 

OPEN(UNIT-15,FILE-’WRCOMPl.DAT’ , STATUS- ’ NEW ’ ) 
C  DX-0.01 

AFEF-S0RT( PREF/RREF) 

TIMREF-XREF/AREF 
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n  n  n  n  n  n 


RMU-SORT( (G-1. )/(G+l . ) » 

X(1)--0.5*DX 
ZETA-MDP( 1 ) 

XII-DX*(WDP(0)-0.5) 

DO  25  1-2,203 

X(I)-X(1-1)+0.5»DX 
25  CONTINUE 

DO  35  1-1,100 

XARRAY(1)-X(1*2+1) 

35  CONTINUE 

INITIAL  DATA 
CALL  INlTl 

CALL  IN1T2L(PSEXIT,PREF,RREF) 

CALL  1N1T2R(PSEX1T,PREF,RREF) 

CALL  IN1T3L( PSINL.RINL) 

CALL  1N1T3R( PSINR,RINR) 

NONDIMENSrONALIZATION 
DO  30  1-1,203,2 
P(  I  )-P( I )/PREr 
R(  I  )-R( I )/RREF 
U( I )-U( I )/AREF 
A( I )-A( I )/AREF 
S(l)-ALOG(P(l)/R(I)»*G) 

30  CONTINUE 
C  CALL  PLOTl(F) 

DO  40  J-1,K 
N-N  +  1 

XII-DX*(WDP(0)-0.5) 

QPRlNT-N/50 

DT-100. 

DO  50  1-1,203,2 

DTT-CFLNUM*DX/(2. ‘AMAXlt ABS( A( I )+U( I ) ) ,ABS(A( 1 )-U(l )  )  )  ) 
DT-AMIN1(DTT,DT) 

50  CONTINUE 

TTOTAL-TTOTAL+DT 
TII1E-TTOTAL»TlKREr 
XI  — XII 

DO  60  1-1,201,2 
PL-P( I  ) 

RL«R( 1 ) 

UL»U( I  ) 

PR-P( 1  +  2  ) 

RR-R( 1+2 ) 

UB-U( 1+2 ) 

XITEHP-XI 

IF(I.EO.l)  XI-ABS(Xl) 

IF( ( I .EO. 201 ) .AND. (XI .GT.0.0) )  XI--XI 
CALL  GLIMM(QSTOP,PSTAR,USTAR, ASTAR) 

XI-XITEMP 
PI I+l )-PGLIM 
R(I  +  1  )-RGHri 
0( I+l )-UGLIM 
60  CONTINUE 

CALL  GG( SWL,SWR,N,TTOTAL, TIME, UEXHAX,PTOTIN,PREF, NASS IN. HASOUT) 
DO  70  1-1,201,2 

IF(XI.LT.0. )  GOTO  80 
P(l+2)-P(l+l> 

C  WRITE(6,*)'P<',I+2,’)-',P(I+2> 

R(I+2)-R(I+l) 

C  WRITE(6, * ) 'B( ' , 1+2, '  )-' ,R( 1  +  2) 
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n  n  n  rj 


80 


70 


40 


U< 1  +  2  )-U( I  +  l ) 

WPITE(6,*)'U(',I+2.')-',U(I+2) 
A( 1+2 ) -SORT tG*P{ 1+2 )/R( 1+2) ) 
WR1TE(6,»)'A( ' ,1+2,' )-' ,A(l+2) 
S( l+2)-AL0G(P( l+2)/R( 1+2 )**G) 
KRlTE(6,*)'S(',I+2,')-',S(I+2) 
GOTO  70 
P(l)-P(l+1) 

R(l)-R(l+1) 

U(l)-U(l+1) 

A(1)-SQRT(G«P<1)/R(I)) 

S(l)-ALOG(P(I)/R(I)»»G) 

CONTINUE 


IF(SWL.EQ.l)  CALL 
ir(SWL.EQ.2)  CALL 
ir(SWL.EQ.3)  CALL 
ir(SWL.EQ.4)  CALL 
ir(SWL.E0.5)  CALL 
ir(SWR.EO.l)  CALL 
IF(SWR.EQ.2)  CALL 
if(swr.eq.3)  call 
ir(swR.EQ.4)  call 
IF<SWR.E0.5)  call 
ir(SWL.E0.4)  MASSIN- 
IF(SWR.E0.2)  hasout- 
IF(N.EQ. (QPR1NT»50) ) 
CONTINUE 
call  ENDPL(O) 
call  DONEPL 


BCLl 

BCL2(PSEXI1,?REF) 

BCL3(PSINL,R1NL,PREF,RREF) 

BCL4 ( PTOTIN , RTOTIN, PREF , RREF ) 

BCL5 

BCRl 

BCR2(PSEX1T,PBEF) 

BCR3( PSINR, RINR , PREF, RREF) 
BCR4 (PTOTIN, RTOTIN, PREF, RREF ) 
BCR5 

MASS1N+R( 3 ) ‘U( 3 ) *DT 
HASOUT+R( 201)*U(201)*DT 
CALL  PLOT2(N,K) 


-  — CALCULATE  PRESSURES  AND  TEMPS  FROM  INPUT  INFO - 

PTE-PTO+TPR 

PE-FTE/( ( 1 .  +  (G-1 . )/2. *AMEXIT**2 .  ) *  * ( G/( G-1 . ) ) ) 

PREF-PE/SPR 

TREF-PREF*14  4 . /( RG1N*RGUESS  ) 

TE-TREF*(SPB**( (G-l. )/G) ) 

TTE-TE *(!.+( G-1. )/2.«AMEXIT**2. ) 

C - CALC  TOTAL  TEMP  RATIO  FROM  REF  DENSITY  GUESS - 

TTR1-TPB**( (G-1 . )/G) 

C 

WP.ITEl  20  ,  •  ) 

KRITE(  20, 3024  ) 

WRITE! 20 , 3025 )SPB,TPR,TTR1, FIN, AMEXIT 
WRITE! 20 , 3026 ) FTO, TTOTIN , G, RTIN, RTOTIN 
WRITE (20, 3027) PTE, PE,PBEr,TTE,TE 
WRITE! 2 0,3028 )TREF,RGOESS, RREF 

3024  FORMAT! 2X, 'WAVE  ROTOR  PERFORMANCE  MAP  ANALYSIS!') 

3025  FORMAT! 5X, ' SPR-' ,F6 .2, 2X, 'TPR  - ' , FS . 2 , 2X , ' TTR  -',F6.4,2X, 
>'riN  ,r6 . 4 , 2x, 'AMEXIT-' ,r6 .4 ) 

3026  FORMAT! 5X, ’PTO-' ,F6.2,2X, 'TOT IN-' ,F8.2,2X, 'GAMMA-' ,r6 . 4,2X, 
>'RTIN-' , re. 4, 2X,' RTOTIN-' ,r6.4) 

3027  FORMAT! 5X, 'PTE-' , re .2,2X, 'PE  - ' , r8 . 2 , 2X , ' PREF  - ' , r6 . 2 , 2X , 

>'TTE  -' ,r7.1,lX,'TE  -',F7.1) 

3028  FORMAT! 4X, 'TREF-' ,r7.1, IX, 'RGUESS-' ,r7. 5, 2X, 'RREF-  ',F6.4) 

RCTIME-SECNDS(RINTT)-RCTIME 

WRITE! 6, * ) 'COHFUTER  RUN  TIME  -  ',RCTIME,'  SeCB.' 

WRITE! 15, *) 'COMPUTER  RUN  TIME  -  ',RCTIME,'  tecs.' 

WRITE! 20, •  ) 

WRITE! 20, ») 'COMPUTER  RUN  TIME  -  ',RCTIME,'  sees.' 
CLOSE!UNIT-15) 
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CLOSE(UN1T-20) 

STOP 

END 


SUBROUTINE  GLIMM{ QSTOP , PSTAR, USTAR , ASTAR ) 

INTEGER  Q, QSTOP 
REAL  ML-MR.HLN.MRN 

COMHON/GLIfim/PGLIM,FGHH,UGLI»l,PL,RL,UL,PR,RR,UR,AL,AR,GL,GB,EPS 

C0MH0N/GLIMM2/DT,DX,XI 

DATA  Q,ML,f1R/0,100.  ,100./ 

PSTAR-0.5*(PL4PR) 

COEFL-SORT( PL*RL) 

C0EFP-SQRT( PR*RR) 

ALPHA-1 . 

C  BEGIN  GODUNOV  ITERATION 

30  0-0+1 

IF(PSTAR.LT.EPS)  PSTAR-EPS 
C  COMPUTE  NEXT  ITERATION  FOR  ML  AND  MR 

MLN-COEFL*PHI( PSTAR, PL) 

MRN-COEFR+PHI ( PSTAR, PR) 

DIFML-ABS(MLN-ML) 

DirMR-ABS(HRN-MR) 

ML-MLN 

MR-MRN 

C  COMPUTE  NEW  PSTAR 

PTIL-PSTAR 

PSTAR- ( UL-UP+PL/ML+PR/MR )/( 1 ./ML+1 ./MR) 

PSTAR-ALPHA*PSTAR+( 1 . -ALPHA) *PTIL 
IF(O.LE.QSTOP)  GOTO  10 
IF(ABS(PSTAR-PTIL) .LT.EPS)  GOTO  20 
C  COMPUTE  NEW  ALPHA 

ALPHA-0 . 5+ALPHA 

0-0 

jr( ( 1 .-ALPHA) .LT.EPS}  GOTO  20 
10  ir(DIFML.GE.EPS)  GOTO  30 
ir(DIFHR.GE.EPS)  GOTO  30 
C  END  OF  GODUNOV  ITERATION?  COMPUTE  USTAR 

20  USTAR-( PL-PR+ML*UL+HR*UR)/(ML+MR) 

C  BEGIN  SAMPLING  PROCEDURE 

IF  (XI ,LT.USTAP*DT)  GO  TO  40 
C  RIGHT  SIDE;  SELECT  CASE  OF  SHOCK  OR  EXPANSION 

IF  ( PSTAR.lt. PR)  GO  TO  50 
C  RIGHT  WAVE  IS  A  SHOCK  WAVE 

WR-UR+MR/RR 

IF  ( XI . LT.WR*DT)  go  TO  60 
C  RIGHT  or  FIGHT  SHOCK  CASE 

RGLIH-FR 
PGLIM-PP 
UGLIM-UR 
RETURN 

C  LEFT  OF  RIGHT  SHOCK  CASE 

60  RGLIM--MR/(USTAB-WR) 

PGLIM-PSTAR 

UGLIM-USTAR 

RETURN 

C  RIGHT  WAVE  IS  A  RAREFACTION  WAVE 

50  CONST-PR/FR**GR 

RSTAR-(FSTAR/CONST)»»(l./GH) 

ASTAR-S0RT(CR«PSTAR/RSTAR) 

AR-SQRT(GR*PR/RR) 
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IF  (XI .GE. (USTAR+ASTAR)*DT)  GO  TO  70 
C  LEFT  OF  RIGHT  FAN  CASE 

RGLIM-RSTAR 
UGLIH-USTAR 
PGLIH-PSTAR 
LTAG-3 
RETURN 

C  SELECT  RIGHT  OF  FAN  OR  IN  FAN 

70  IF  (XI.GE. (UR+AR)*DT)  GO  TO  60 
C  IN  RIGHT  FAN  CASE 

UGLIM-2 ./(GR+1 . )*(XI/DT-AR4(GR-1. )/2.*0R) 

RGLIH-( ( (AR+(GR-1. ) /2 . * ( UGLIM-UR ) )«*2. )/(GR* CONST) ) * • ( 1 . /( GR-1 . ) ) 
pglih-const*rglim*»gr 

RETURN 

C  RIGHT  OF  RIGHT  FAN  CASE 

80  RGLIH-RR 
PGLIH-PR 
UGLIM-UR 
RETURN 

C  LEFT  SIDE;  SELECT  CASE  OF  SHOCK  OR  RAREFACTION 

40  IF  (PSTAR.LT.FL)  GO  TO  90 
C  LEFT  WAVE  IS  A  SHOCK  WAVE 

WL-UL-ML/RL 

IF  (XI .GE.WL*DT)  GO  TO  100 
C  LEFT  OF  LEFT  SHOCK  CASE 

RGLIH-RL 
PGLIM-FL 
UGLIH-UL 
RETURN 

C  RIGHT  OF  LEFT  SHOCK  CASE 

100  RGLIM-r.L/(USTAR-WL) 

PGLIM-PSTAR 

UGLIH-USTAR 

RETURN 

C  LEFT  WAVE  IS  A  RAREFACTION  WAVE 

90  CONST-f  L/Hl‘' “OL 

RSTAR-(  PSTAF./CONST)  •«  (  1  ./GL) 

A S TA R -  3 OR T  G „  ‘  P S T AR/RS TAR  ) 

AL-SQRT(GL*PL/RL) 

IF  (XI  .LT.  ('.i<  TAR-ASTAR)»DT)  GO  TO  110 
C  Rir  'T  OF  LEFT  FAN  CASE 

RGLIH-RSTAR 
PGLin-PSTAP, 

UGLIH-USTAR 

RETURN 

C  SELECT  LEFT  OF  FAN  OR  IN  FAN  CASE 

110  IF  (XI .LT. (UL-AL)»DT)  GO  TO  120 
C  IN  LEFT  fan  CASE 

UGLi:i-2./(&^J^i.  )*(AL+(GL-1.  )/2  . ‘UL  +  XI/DT  ) 

RGLIH-( ( ( AL+(GL-1 . ) /2 . * ( UL-UGLIH ) )**2. )/(GL*CONST) )*•( 1 ./(GL-1 . ) ) 

PGLIH-CONST*RGLIM»*GL 

RETURN 

C  LEFT  OF  LEFT  FAN  CASE 

120  RGLIH-RL 
PGLIH-FL 
UGLIH-UL 
RETURN 
END 

. . . . . . **************** 

FUNCTION  PHI(y,Z) 
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REAL  RHU 

C0MH0N/FUN1/G . PA , RA , UA , RB , RHO 

EPS-1 . E-06 

PARAM-Y/Z 

IF  ( ABS( 1 .-PARAH) .GE.EPS)  GO  TO  10 
PHI-SQRT(G) 

RETURN 

10  IF  < PARAH. GE. 1 .  )  GO  TO  20 

PHI-(G-1 . )/2 .*( 1 .-PARAH)/(SQRT(G)*(1 .-PARAM**( < G-1 . ) /( 2 . *G ) ) ) ) 
RETURN 

20  PH I- SORT ( <G+1 . )/2.*PARAM+(G-1. )/2. ) 

RETURN 

END 

FUNCTION  PHIUPB) 

REAL  RHU 

COHHON/FUNl/G , PA , PA , UA , RB , RHU 

PHI1-(PB-PA)»SQRT( (l.-RMU**2. )/( RA* ( PB+RHU* * 2 . * PA ) ) ) 

RETURN 

END 


FUNCTION  PSI(PB) 

REAL  RHU 

COHMON/FUN 1/G , PA , RA , UA , RB , RHU 

PS  1 -SORT ( 1 .-RHU* *4 .  !/RKU**2./SQRT(RA)*PA**< 1 ./(2 . *G ) ) * ( PB* * ( ( G-1 . ) 
>/(2. ‘G) )-PA**( <G-1 . )/(2.*G) ) ) 

RETURN 

EHC 

SUBROUTINE  INlTl 

DIHENSION  P(203),R(203),U(203),A(203),S(203),X(203) 

COriMON/FUNl/G,PA,RA,UA,RB,RHU 

COMMON/SUBS /P,R,U, A, S,X 

DO  10  1-1,9, 2 

P(I)-830600.00 

B( I )-0.7132 

U( I )-644 . 4 

A( I )-SQRT(G*P( I )/R( 1  )  ) 

10  CONTINUE 

DO  20  1-11,203,2 
P( I )-101325.0 
R( I )-l . 22 
U{ I )-0 .0 

A(I)-SQRT(G*P(I)/R(1)) 

20  CONTINUE 
RETURN 
END 

Q*  . . 

SUBROUTINE  INIT2R( P5EX I T , PREF , RREF ) 

DIMENSION  P(203),R(203),U(203),A1203),S(203),X(203) 

COMHON/FUNl/G, PA,RA,UA,RB,RHU 

COMHON/SUBS/P,R,U,A,S,X 

DO  10  1-3,201,2 

P( I )-PBEr 

R( I  )-BREF 

U( I  )-0 . 0 

A(I)-S0RT(P(I)*G/B(1)) 

10  CONTINUE 
P( 1)-P( 3) 

R(l)-B(3) 

U(l)— U(3) 
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A{ 1 )-SQBT(G*P(l )/R(l ) ) 

P( 203 )-PSEXIT 
R(203)-n(201) 

PA-P(201  ) 

RA-R(201  ) 

OA-U(201) 

PB-P(203) 

PB-R(203» 

IF(PA.GT.PB)  GO  TO  20 
U(203)-UA-PHI1(PB) 

GO  TO  30 

20  U(203)-UA-PSI(PB) 

30  A(203»-SORT(G«P(203)/R(203) ) 

RETURN 

END 

Ct********»*****»********»»»***»»»»t,*********»*»»**»**»»*»»*i, 

SUBROUTINE  I NIT2L( PSEXIT ) 

DIMENSION  P(203),B(203),U(203>,A(203),S(203),X(203) 

COKMON/FUN 1/G , PA , RA , UA , RB , RMU 

COHMON/SUBS/P,F,0,A,S,X 

DO  10  1-3,201,2 

P(  I  >-285080. 0 

R(  I)-0.e97 

U(I)-0.0 

A( I )-SQRT(G*P(I )/R( 1 > ) 

10  CONTINUE 
RETURN 
END 


SUBROUTINE  I NIT3L( PSINL , RINL ) 

DIUENEIC-  P(203) ,R(203),U(203) ,A(203) ,S(203) ,X(203) 

COMHON/FON1/G,PA,RA,OA,RB,RMU 

COHMON/SUBS/P,R,0,A,S,X 

DO  10  1-3,201,2 

P( I >-2390000.0 

R(I  >-9.7H7 

U(1)-0.0 

A(  I  >-S  JRT(G*P(  I  >/R(  I ) 

10  CONTINUE 
P',  1  >-FrT  'IL 
R( 1 )-RIND 
PA-P(3) 

RA-R( 3  > 

UA-U( 3  ) 

PB-P( 1 ) 

U( 1 )-UA+PHIl ( PB  > 

A(1>-SQBT(G*P(1>/R(1)) 

P( 203 )-P( 201 > 

R(203>-R(201) 

U(  203  >— U(  201  > 

A(203>-SORT(G*P(203)/B(203) > 

RETURN 

END 


SUBROUTINE  1NIT3R( P5INR , RINR > 

DIMENSION  P(203> ,R(203> ,U(203> ,A(203> ,S(203> ,X(203> 
COMMON /SUBS/P, R,U, A, S,X 
COMMON/rUNl /G , PA , RA , UA , RB , RMU 
DO  10  1-3,201,2 
P( I >-2421667.5 
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R(I)-9.787 

U(I)-0.0 

A( I )-SQRT(G*P( I )/R(I ) ) 

10  CONTINUE 

P( 203 )-PSINR 
PA-P{ 201 ) 

RA-R( 201  ) 

UA-U( 201  ) 

PB-P(  203  ) 

U(203)-UA-PH11(PB) 

R( 203 )-RINR 

A(203)-SORT(G‘P(203)/R(203) ) 
P(l)-P(3) 

R(l)-R(3) 

U(l)— U(3) 

A( 1 )-SQRT(G*P(l )/R(l) ) 

RETURN 

END 


SUBROUTINE  BCLl 

DIMENSION  P(203),R(203),U(203),A(203),S(203),X(203) 

COMMON/SUBS/P,R,U,A,S,X 

COMMON/FUNl/G ,  PA ,  RA,  OA ,  RB  ,  RHU 

P(l)-P{3) 

R(l)-R( 3) 

U( 1 )--U( 3) 

A( 1 )-SQRT(G*P( 1 )/R(l ) ) 

RETURN 

END 

C‘  *  *  •  * 

SUBROUTINE  BCR3 ( PS INR , RI NR , PREF , RREP ) 

DIMENSION  P(203),R(203>,U(203),A(203),S(203),X(203) 

COMMON/SUBS/P, R,U, A, S,X 

COMMON/FUN 1 /G , PA , RA , UA , RB , RMU 

r(  203  »-PSINR/PREr 

B( 203 )-RINB/RREF 

PA-P( 201 ) 

RA-R( 201  ) 

UA-U( 201  ) 

PB-P(  203  ) 

U{203)-UA-PHI1(PB) 

A(  20  3  )-SQBT(G*P(203)/R(203)  ) 

RETURN 

END 

C*  *  * ♦ * 

SUBROUTINE  BCL3(PSINL,RINL,PREF,RPEF) 

DIMENSION  P(203),R(203),U(203),A(203),S(203),X(203) 

COHMON/SUBS/P,R,U,A,S,X 

COMHON/FUNl /G , PA , RA , OA , RB , RHU 

P( 1 )-PSINL/PREF 

R( 1 l-RINL/RREF 

PA-P( 3) 

RA-R( 3) 

UA-U( 3  ) 

PB-P( 1 ) 

U( 1 )-UA+PHIl ( PB> 

A( 1 )-SQRT(G»P( 1 )/R( 1 ) ) 

RETURN 

END 
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SUBROUTINE  BCRl 

DIMENSION  P{203|,R(203),0(203),A(203),S(203),X(203) 
COMMON/SUBS/P, R,0, A, S,X 
COMMON /rVNl/G , PA , RA , UA , RB , SHU 
P(2C3)-P(201) 

R(203)-R(201) 

U(203)— U(201) 

A(203)-SQRT(G*P{203)/R(203) J 

RETURN 

END 

subroutine  BCL2(PSEXIT,PREr) 

DIMENSION  P(203),R(203),0(203),A(203),S(203),X(203) 

COMHON/SUBS/P.R.U, A,S,X 

COMMON/FUNl/G,PA,RA,OA,RB,RMU 

P( I )-PSEXIT/PREr 

R(l)-R(3) 

U(l  )-U( 3) 

A(1  I-S0RT(G*P(1)/R(1) ) 

RETURN 

END 


SUBROUTINE  BCR2 ( PSEXIT , PREF ) 

DIMENSION  P(203),R(203).U(2C3).A(203),S(203),XU03) 

COMMON/SUBS/P,R,U, A,S,X 
COMMON/FUNl /G , PA , FA, UA, RB , RHU 
P(  203  )-PSEXIT/PREF 
R(203)-R(201) 

FA-P( 201 ) 

ra-r;201  ) 

UA-U( 201 ) 

PB-F( 203 ) 

RB-R( 203 ) 

IF(PA.GT.PB)  GO  TO  10 
U(  203  )-UA-PHIl(PB) 

GO  TO  20 

10  U(203)-UA-PSI(PB) 

20  A(203)-SORT(G*P(203)/B(203)) 

RETURN 

END 


FUNCTION  WDP(II) 

DIMENSION  WNORH( 12 ) , 1DIGT( 12) 

COMMON/SAMFLE/WNORM, IDIGT 

IF  (II.EQ.O)  GO  TO  10 

Ll-2 

L2-1 

DO  20  JJ-1,12 
IDIGT(00 )-0 

WNORH( JJ)-1 ./FLOAT! Ll** JO) 

20  CONTINUE 
WDP-0. 

RETURN 

10  DO  40  JO-1,12 
Ll-2 
L2-I 

RJ0-1DIGT( JJ) 

RJN-MOD( ( KJO  +  I ) ,Ll  ) 

IDIGT! JO )-KJN 
IF  !KJ0.LT.KJN)  GO  TO  50 
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40  CONTINUE 
50  SUM-0. 

DO  60  JJ-1,12 

KNEW-MOD( IDIGT( JJ ) *L2,  Ll ) 

SUM-SUM+  FLOAT ( KNfcW ) ‘WNORHC  J J ) 

60  CONTINUE 
WDP-SUH 
RETURN 
END 

C****************************************«********«******************** 

C  SUBROUTINE  PLOT2(N,K) 

C  DIMENSION  XORG( 4 ) ,yORG( 4 ) ,¥MAX( 4 ) , VMIN( 4) ,KNT( 4 ) ,iyNAM( 10) 

C  DIMENSION  PARRAV ( 100 ) ,RARRAy( 100 ) ,UARRAy( 100) ,SARRAy( 100 ) ,XARRAY( 1 

C  >00 ) 

C  DIMENSION  P(203),R(203),U(203),A(203),S(203),X(203) 

C  COHMON/SUBS/P,R,U,A,S,X 

C  COMMON  XARRAy 

C  DATA  XORG/0.5,4.75,0.5,4.75/ 

C  DATA  yORG/0.5.0.5,4.75,4.75/ 

C  DATA  yMAX/2.5b,2.5,2.5,2.5/ 

C  DATA  YMIN/O. 0,0.0, -1.0, -2.0/ 

C  DATA  KNT/1,4,6,9/ 

C  DATA  lyNAM/'PRES' , 'SURE' , ' DENS ' , ’ ITYf ' , ' VELO ’ , ' C I TY ' , ' $  ' 

C  > , ' ENTR’ , 'OPy$ '/ 

C  DO  200  1-1,100 

C  PARRAYC 1 )-P( 1*2  +  1  ) 

C  RARRAy( I )-R( 1 *2  +  1  ) 

C  UARBAy( I )-U( 1*2+1 ) 

C  SARRAy( I )-S( 1*2+1  ) 

C  200  CONTINUE 

C  DO  300  1-1,4 

C  CALL  PHySOR(XORG(l),YORG(l)) 

C  CALL  AREA2D( 3 . 5, 3 . 5) 

C  CALL  XNAME(  ’X'  ,1  ) 

C  CALL  YNAMEI iyNAM(RNT( I ) ) ,100) 

C  CALL  GRAr(0. , 'SCALE' ,1 .0, y«IN( I ), 'SCALE' ,yMAX( I ) > 

C  IF(I.EQ.l)  CALL  SETCLR( ' YELLOW' ) 

C  ir(I.EQ.2)  CALL  SETCLR( 'CYAN' ) 

C  ir{I.EQ.3)  CALL  SETCLR ( ' BED ' ) 

C  IF(I.EC.4)  CALL  SETCLRl 'MAGENTA' ) 

C  IFIN.EO.R)  CALL  SETCLR (' WHITE '  ) 

CC  IF{N.EQ.2398)  CALL  SETCLR( 'BLUE' ) 

C  IFd.EO.l)  CALL  CURVE  (  XARRAY ,  PARRAY ,  3  00 , 0  ) 

C  IF(I.EQ.2)  CALL  CURVE  ( XARRAY , RARRAY , 1 0 0 , 0 ) 

C  IF(I.EQ.3)  CALL  CURVE  ( XARRAY , UARRAY , 1 00 , 0 ) 

C  1F(I.EQ.4)  CALL  CURVE  ( 3CARRAY ,  SARRAY  ,  1 00 , 0  ) 

C  CALL  ENDGR(O) 

C  300  CONTINUE 
C  RETURN 

C  END 

. . ♦♦«*.*»»«***.*»*..*.**.*.********.***** 

C  SUBROUTINE  PLOTl(K) 

C  DIMENSION  XOBG( 4 ) , YORG( 4 ) , YMAX( 4 ) , YMIN( 4 ) 

C  DATA  XORG/0.5,4.75,0.5,4.75/ 

C  DATA  YORG/0.5,0.5,4.75,4.75/ 

C  DATA  YHAX/2.50,2.5,2.5,2.5/ 

C  DATA  YMlN/0. 0,0.0, -1.0, -2.0/ 

C  DO  10  ;-l , 4 

C  CALL  PHYbC'»(XORG(  I  )  ,  YORG(l )  ) 

C  CALL  AREA2D( i . 5, 3. 5) 
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nnnnnnnonnnn 


CALL  FFAME 

CALL  GRAF(0. , ’SCALE’ , 1 . 0 , YHI N ( I ) , ' SCALE ’ ,YMAX(I) ) 
CALL  ENDGR(O) 

10  CONTINUE 

CALL  PHYSOR(e.5,0.5) 

CALL  AREA2D(2.25,7.75) 

CALL  FRAME 

CALL  GRAF ( 0 . , ’ SCALE ’ , 1 . , 0 , ’SCALE ' , K ) 

CALL  ENDGR(O) 

RETURN 

END 


SUBROUTINE  BCL4  ( PTOTIN , RTOTIN , PREF , RREF ) 

INTEGER  QOUT 

DIMENSION  P( 203) ,R(203),U(203),A(203) ,S(203) ,X(203) 
DIMENSION  QR( 3  ) 

DIMENSION  XARRAY(IOO) 

COMMON/SUBS/P,R,U,A,S,X 
COMNOU/rUNl/ G, PA.RA.UA.RB.RMU 
COMMON  XARRAY.Nl 
Nl-NJ  +1 

QP( 3)-2./(G-l. )*A( 3)-U(3) 

QOUT-Nl /5 
PTOT-rrOTIN/PREF 
RTOT-PTOTIN/RREF 
ATOT-SORT(G*FTOT/BTOT) 

STOT'ALOG( PTOT/RTOT»»G) 

U( 1  !-U( 3  ) 

A( 1 )-S0RT( ATOT**2 .-(G-1 . )/2.*ABS(U( 1 ))»*2.) 
0P(l)-2./VC-l.)»A(l)-U(l) 

AMACH-U( 1 )/A( 1 ) 

I F( AKACH.lt. 0.0)  GO  TO  60 

P( 1 )-PTOT/l 1 .  +  (G-1 . )/2 .•AHACH»*2 . ) *  *  (  G/(  G- 1 .  ) ) 

R( 1 )-RTOT/( I .  +  (0-1 . )/2. ‘AMACH**! . )**{!. /{  G-1.)) 
S(l)-ALOO(P(l)/R(l)»*G) 

GO  TO  50 

60  P( 1 )-PT01/a .+(G-1 • )/2 . *ABS(AMACH) **2 . ) • * ( G/( G-1 . ) ) 
B( 1 )-RTOT/( 1 .  +  (G-1 . )/2.*ABS(AHACH)**2. ) *  * ( 1 ./( G-1 . ) ) 
S( 1 )-ALOG( r ; 1 )/B{l )**G) 

50  IF(N1  .E0.5*QOUT)  WPITEdO,!)  A(  1 )  ,  A(  3  )  ,  U  (  1  )  ,  U  (  3  )  ,  Nl 
1  rOPHATI  2.r,'(F9.6,5X,l5) 

:  TURN 
END 

•  *  * 

SUBROUTINE  BCR4  ( PTOTIN , FTOTI N , PREF , RREF ) 

INTEGER  QOL'T 

DIMENSION  P(203),R(203),U(203),A(203),S(203),X(203) 
DIMENSION  QL( 203  ) 

DIMENSION  XARRAYdOO) 

COrtMON/SUd.9/  P  ,  R  ,  U ,  A ,  S  ,  X 
COMHON/rUNl/G,PA,RA,UA,RB,RMU 
COMMON  .XARRAY.Nl 
N1-NIh 1 

QL{201)-2./(G-l.)*A(201)+U(201) 

QOUT-Nl/25 
PTOT-PTOTIN/PPEE 
RTOT-RTOTIN/PREF 
ATOT-S0RT( G‘FTOT/RTOT) 

STOT-AI.OG(  PTOT/RTOT»*G) 

U(203I-U(401) 
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20  3  |-SQRT(AT0T*»2  .-(G-1  .  )/2  . ‘ABSCyi' 203  )  )**2.  ) 
QL(2C3)-2./(G-1. ) »A( 203 » +0(203) 

AHACH-U(  203  )/A(  203  ) 

Xr(AMACH.LT.O.O)  GO  TO  60 

P(203)-PTOT/(1.+(G-1.)/2.*AMACH**2. )**(G/(G-1. ) ) 

P( 203)-BTOT/( 1 .+(G-1. )/2.*AMACH**2 .)**{! ./(G-1 . ) ) 
S(203)-ALOG(P(203)/R(203)**G) 

GO  TO  50 

60  P( 203 )-PTOT/l 1 .+(G-1 . )/2 .*ABS(AMACH) **2 . ) **(G/{G-1 . ) ) 
R( 203 )-RTOT/( 1 .+(G-1 . ) /2 . » ABS ( AMACH ) * • 2 . )*♦(! ./(G-1 . ) ) 
S(203)-ALOG(P(203)/R(203)»»G) 

50  tF(Nl . EQ. 25*OOUT)  WRITEllO,!)  P( 201 ) ,U( 201 ) ,Nl ,N 
1  rORHAT(2X,2rl3.6,5X,2l5) 

RETURN 

END 


SUBROUTINE  BCL5 

DIMENSION  P( 203) ,R(203) ,U(203) ,A(203 ) , S ( 20 3 ) , X( 203 ) 

COMMON/SUB S/P, R,U, A, S,X 

P(l)-P(3) 

P( 1 )-R( 3) 

U(l)-U(3) 

A(1 )-A( 3) 

RETURN 

END 


SUBROUTINE  BCP5 

DIMENSION  P(203),R(203),U(203),A(203),S(203),X(203) 
COMMON/SUB S/P,R,U, A, S,X 

P(  203  )-P' 201  ) 

R(  203  )-R( 201  ) 

U(203)-U(201) 

A(203)-A(201) 

RETURN 

END 


l**««**<«*«****<««**<**«**< 


SUBROUTINE  C” ( SWL, SWR, N , TTOTAL , T1 ME , UEXMAX , PTOTIN , PREF , MASS IN , MASO 
‘ur ) 


INTEGER  SWL, SWR 
P.  M,  MASSIN,HASOUT 

DIMENSION  P(203),R(203),U(203),A(203),St203),X(203) 
COMMON/SUBS/F , F , U , A , S , X 
IF( (SKL.EQ.l) .AND.(SKR.E0.2))  GO  TO  10 
IF ( 1 SWL . EQ. 4 ) , AND, ( SWF . EQ. 2 ) )  GO  TO  30 
IF<  (  SWL.EQ .4 ) . AND. (SWR.EQ.l ) )  GO  TO  50 
I F( ( SWL , EQ . 1 ) .AND. ( SWR.EQ.l ) )  RETURN 
10  rWAI,t.»r(2) 

ir(PWALL.LE. (PTOTIN/PREF) )  GO  TO  20 
RETURN 
20  SWL-4 

WRITE) 6, 74  ) 

WRITE(6,75)  N.TTOTAl .time 
WRITE( IS , 74  ? 

WRITE(15,75)  N, TTOTAL, TIME 
WRITE( 20  .  •  ) 

KRITE( 20 ,♦) 'TIMING  AND  MASS  FLOW  BALANCE:' 

WRITE(  20,74  ) 

KP,ITE(20.75)  N, TTOTAL, TIME 
RETURN 

30  UEX1T-U(202) 
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IF(UEXHAX.LT.OEXIT)  UEXMAX-UEXIT 
ir(UEXIT.LT.UEXMAX/2. )  GO  TO  40 
RETURN 
40  SWR-1 

WRITE(6,76) 

WRITE(6,75)  N.TTOTAL.TIME 
WRITE( 15,76  ) 

WRITE(15,75)  N.TTOTAL.TIME 
WRITE (20, 76) 

WR1TE(20,75)  N,TT0TAL,TIME 
RETURN 

50  P1SH0K-P(2) 

C  IF(N.EQ.1590)  GO  TO  60 

IF(PlSHOK.GT.PTOTIN/PREF)  GO  TO  60 
RETURN 
60  SWL-1 

WR1TE(6,77) 

WRITE(15,77) 

WP.ITE(  20,77  ) 

wr:te(6,75)  n,ttotal,time 
write(15,75)  n,ttotal,tihe 
WRITE{20,75)  n,ttotal,tihe 
WRITE(6,78)  P(199),P(201),P(203) 

WRITE (15. 78)  P(199) ,P(201) ,P(203) 

WRITE(20,78)  P(199),P(201),P(203) 

WRITE( 6 , 73 )  MASSIN,MASOUT,SWL,SWR 
WRITE( 15 , 73 )  HASSIN,HASOUT,SWL,SWR 
WRITE( 20, 73 )  HASSlN,HASOUT,SWL, SWR 

73  FORMAT (5X, 'Min-' , F9 . 6 , 2X , ' Mou t- ' , F9 . 6 , 5X , 'SWL-' ,I4,2X, 

74  rORHAT(  5x,  '  INLET  PORT  OPENS  AT:') 

75  FORMAT  (  5x,  'N-'  , 14 , 5X, ' TTOTAL- ' ,F14.7,2X,'TIME-',ri4.7) 
78  FORMAT cx, 'P(199)-',r8.6,2X,'P(201)-',Fb.6,2x,'P(203)- 

76  rORMAT( 5X, ' EXHAUST  PORT  CLOSES  AT:’) 

77  rORMAT( 5X, ' INLET  PORT  CLOSES  AT;') 

RETURN 

END 


’SWR-' ,14 

' ,re.6) 


87 


APPENDIX  C 


COMPUTER  SYSTEM  AIDS 


A.  TRANSFERRING  FILES 

1.  FROM  MAINFRAME  TO  AERO  PC  LAB 

The  following  procedures  use  SIMPC  to  download 
programs  from  the  NPS  mainframe  to  the  Aero  Department's  PC 
Lab.  The  modem  and  SIMPC  hooked  up  to  the  IBM  PC-XT-2  in  the 
Aero  PC  lab.  The  following  boldfaced  commands  outline  the 
procedures . 

C\  <dir>  >cd\SIMPC . change  dir  to  SIMPC  dir 

C\SIMPC>SIMPC  ACCESS. SIM. . .execute  SIMPC  macro  ACCESS. SIM 

•  •  • 

«  «  « 

LOGON:  9812P  ....logon  to  mainframe  (user  number) 

PASSWORD:  _  ....enter  your  mainframe  password 


<enter>  ....executes  your  PROFILE. EXEC 


FORSIMPC  GET  ENGINE  FORTRAN  Al..file  to  copy  from 

mainframe 


A: ENGINE. FOR 


FORSIMPC  PRINT 
to 

LOGOFF 
<ALT>  <G> 

<ALT>  <Q> 


ENGINE  FORTRAN  Al.. optional  if  you  wanted 

send  your  file  to  the  current  printer. 

logoff  the  mainframe 

hangs  up  the  modem 

quits  SIMPC  and  returns  to  DOS 
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The  procedure  above  uses  the  macro  ACCESS. SIM  which 
was  written  specifically  to  link  between  the  Aero  PC  lab  nd 
the  NPS  mainframe. 

2.  FROM  AERO  P.C.  LAB  TO  VAX  SYSTEM 

The  following  procedure  uses  PROCOMM,  an  executable 

modem  link,  and  ALT  0,  a  modem  macro  written  specifically  to 

link  the  Aero  PC  lab  to  the  Aero  VAX  lab.  KERMIT  is  the 

actual  file  transferring  program. 

C\  <dir>  >cd\COM. .. change  to  COM  directory. 

FROCCKK  .executable  program 
<ALT>  0  .macro  modem  link 
ATDT2953 
CONNECT 

USERNAME:  GUEST  .gouge  username 
PASSWORD:  _  .password  will  not  appear 

Are  you  logged  on  from  a  PC?  Y 
Do  you  want  to  connect  to  WASP?  Y 

Username:  _  .typical  VAX  logon 

Password: 


C\>KERMIT  ..executes  file  transfer  program 
to  send  a  file: 

Receive  >RCOMPl . FOR. .. Program  you  wish  to  call 
F3:  WR!.  P.FOR  ...  file  you  wish  to  send  but  this  file 
must  be  in  the  PROCOMM  dir  on  the  C  (hard)  drive! 

LOGOFF  ...logoff  the  VAX 

<ALT>  X  ...exits  PROCOMM 


If  you  have  any  difficulties  FIO  is  the  help  menu! 
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B.  USING  "SIGHT"  UTILITY  ON  THE  VAX 


SIGHT  is  a  utility  drawing  program  available  on  the  Aero 
Department  VAX  computer.  It  was  used  to  create  pop~up  menus 
for  the  configurations  for  the  ENGINE  program.  To  run  SIGHT, 
at  the  DCL  prompt  enter 
$SIGHT 

A  pop-up  menu  provides  simple  instructions.  Save  your  drawing 

as  a  <fn>.UIS  file  and  then  use  the  VAX  "render"  command  to 

create  your  output  file  as  a  <fn>.REN  file.  An  example  for 

sending  creating  and  sending  an  output  file  to  the  laser 

(ln03)  printer  is: 

$RENDER  <fn>.UIS/DEVICE_TYPE=ln03 
$PRINT/PASSALL/QUE=lnC3  fn.REN 

SIGHT  can  be  used  to  generate  a  picture  from  inside  a  program. 
This  can  be  done  by  inserting  the  following  statement  in  a 
(FORTRAN)  program. 

CALL  LIB$SPAWN  ('RENDER  <fn>.UIS') 

This  was  done  for  the  ENGINE  program. 
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C.  GRAFkit  PROCEDUREf 

PLOT3.FOR  is  a  GRAFkit  plotting  program  which  was 
incorporated  in  the  ENGINEIA.FOR  program.  PLOT 3 . FOR  will  be 
available  as  an  engine  performance  plotting  program  on  the  VAX 
network.  It  can  accept  up  to  50  separate  data  inputs  for  a 
selected  X-axis,  a  Yl-axis  for  specific  thrust,  and  a  Y2-axis 
for  specific  fuel  consumption.  The  instructions  for  running 
PLOT 3 . FOR  are  included  in  the  comments  section  of  the  program 
listing. 


PLOT3 


)roaram  listinc 


C - Tins  PLOTTING  PROGRAM  OSES  GRAFEIT  SOFTWARE - 

C - EdSL'PE  THE  FOLLOWING  SEOO'NCE  IS  ENTERED  PRIOR  TO 

C  EXECUTION  OF  PLOTS: 

c  ligon  thp  VAX 

C  SGrSETUP 

C  $UIS 

C  $nUN  PLOTS 

C - IF  YOU  WANT  YOUR  RESULTS  SENT  TO  EITHER  THE  LASER  OR 

C  HIGH  SPEED  PRINTERS: 

C  logon  the  VAX 

C  $GESETUP 

C  $ln03 

C  SOEFINE  GR_OOT  PLOT.DAT  ("PLOT. DAT"  IS  YOUR) 

C  $ln03s  ep  (OUTPUT  fn.ft  ) 

C  SPUN  PLOTS 

C . to  Bend  output  to  the  selected  printer: 

C  $FRINT/PASSALL/OUE-LNOS  PLOT.DAT 


c - 

PROGRAM  PLOTS 

REAL  Xl(50),Yl(50),Y2(50) 

CHARACTER  XI LABEL* 1 0 , Yl LABEL*  3 3 , Y2 LABEL* S S 


CHARACTER  F I G1 *  7 , NUMBl * 3 , TOPLBL * 1 0 

C - .............................................. 

c - SPECIFY  X.Y  LEFT  AND  Y  RIGHT  LABELS - 

C 

YlLABEL.'A:  Specific  Thrust  Ibf /( Ibis/s  ) ' 
YiLABEL-’B:  SFC  (Ibm  fuel/ht)/lbf  thrust’ 
c - CHOOSE  THE  X  PARAMETER  TO  PLOT 


MRITE(  6,  *) ’CHOOSE  YOUR  X  AXIS  LABELd.e.  1):’ 
WRITE! 6,*)'  1.  BYPASS  RATIO' 

WRITE(6.*)'  2.  FAN  PRESSURE  RATIO' 

WRITE(6,*)'  3.  COMPRESSOR  RATIO' 

WRITEie,*)'  4.  LPT  EXIT  MACH  NUMBER’ 

WRITE (6,*)'  5.  OTHER  (MAX  10  CHARACTERS)’ 

READ( S, * ) lANSl 
IF  (lANse.EQ.i)  then 
X1LABEL-’  BETA  ' 

ELSEIF  (IAN54.E0.2)  THEN 
XILABEL-'  PRLC  ’ 

ELSEIF  (IANS4.E0.3)  THEN 
XILABEL-'  PRHC  ' 

ELSEIF  (1ANS4.E0.4)  THEN 
XILABEL-'AHACHS' 

ELSEIF  (IANS4.E0.5)  THEN 
W*  **E(C,*J'  INtUl  A  AaIS  LACE':’ 
fiEAD(5,’(A10)' )X1LABEL 
END  IF 
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uuuuu  u  u  u  uuuuuuuuou 


- CHOOSE  THE  NUMBER  OE  DATA  RTS  TO  PLOT  NAX-50 - 

NRITE(6,*)’INPUT  NUMBER  OF  DATA  POINTS  TO  PLOT  (MAX-50)’ 
BEAD(5,»)  ICPl 
- INPUT  X  AXIS  VALUES: 

NRITE( 6 , * ) ’ INPUT  VALUES  FOR  '.XlLABEL 
DO  I-1,ICP1 

NR1TE(6,1000)’X(*,I,’) 

REAO(5,«)Xl(I) 

END  DO 


NRITE(«,»)'INPUT»  ' 

HRITE(6,*)’  0  -  AUTOMATICALLY  DETERMINE  XlHIN  AMD  XINAX' 

KRITE(6,«)’  1  -  INPUT  XlMIN  AND  XIMAX’ 

REAO(5,*)IANS5 
IF  (IANS5.E0.1)  THEM 

NRITEIS,*)'  INPUT  XIMIN  -t' 

P.EAD(5,«  )X1HIN 
NR1TE(6,*)’  INPUT  XIMAX  -i' 

READ(5,*)X1MAX 

-  AUTOMATICALLY  DETERMINE  THE  MIN  AND  MAX  X  VALUE 

ELSE 

XIMIN-Xl ( 1 ) 

XIMAX-Xl ( I ) 

DO  J-1,ICP1 
XlTEMP-Xl ( 0 ) 

X2TEMP-Xl(a) 

IF  (XIMIN. GT.XlTEMP)  THEN 
XlMIN. XITEMP 
END  IF 

IF  (X1MAX.LT.X2TEMP)  THEN 
X1MAX-X2TEMP 
ENDIF 
END  DO 
ENDIF 

- INPUT  VALUES  FOR  YI  AND  Y2  ARRAYS 

NRITE(6,»)’ INPUT  VALUES  FOR  '.VILABEL 
DO  J-l,ICPl 

NRITE(6,1000) 'Ylt • , J. ' ) 

READ(5,»)Y1(J) 

END  DO 

NRITE(6,»)'lNPOt  VALUES  FOR  ' ,Y2LABEL 
DO  E-l,ICPl 

KRITE(6, 1000) ' Y2( ’ ,K, ' ) 

READ(5<«)Y2(K) 

END  DO 

lOOC  FORMAK  1X,A3.  I2.A3) 


INPUT  THE  NUMBER  OF  THE  FIGURE  TO  PLOT  (FIG  1) - 

MBITE(6,*) • INPUT  THE  FIGURE  NUMBER  (ie  1  OR  2...)' 

READ( S, ' (A9) ' )NUMB1 
NRITE( 6 . • ) ’NUHBl-' ,NUMBl 
FIGl-'FIGURE  ' 

TOPLBL-FIG1//NUMB1 

NRITE(6,M1X,A11)’ )'TOP  LABEL  -'.TOPLBL 


CONSTRUCT  THE  PLOT 


WRITE! 6,*) 'STANDBY  FOR  PLOT' 

OPEN  GKS,  OPEN  AND  ACTIVATE  WORXSTATION. 
CALL  OFNGKSd.a.l) 

SET  OP  THE  COLOR  TABLE 

CALL  CSCRd, 0,0. 00, 0.00, 0.50) 

CALL  GSCRd, 1,1. 00, 1.00, 0.00) 

CALL  GSCRd, 2, 0.00, 1.00, 0.50) 

CALL  GSCRd, 3, 0.00, 1.00, 1.00) 

CALL  GSCRd, 4, 1.00, 0.00, 0.00) 

CALL  GSCRd, 5, 1.00, 0.00, 1.00) 

CALL  GSCRd, 6,1 .00,0.00,1.00) 

CALL  GSCRd, 7, 1.00, 1.00, 1.00) 
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CALL  GSCRd, 8,0.00, 0.00, 0.00) 

C 

C - SET  VIEWPOINT  AND  DEFINE  THE  RANGE  OF  USER'S  DATA  TO  BE  USED 

CALL  AGSETF  ( ' GRID/LEFT . ' ,  .15) 

CALL  AGSETF  ( 'CRIO/RIGBT. • ,  .85) 

CALL  AGSETF  { • GRID/TOP . ' ,  .85) 

CALL  AGSETF  ( 'GRID/BOTTON. ' ,  .15) 

C - SET  UP  FOR  FIRST  PLOT 

C -  DEFINE  TRE  TOP,  BOTTOH,  LEFT  AND  RIGHT  LABELS 

C -  PUT  BLANES  IN  FIRST  CURVE  LABELS 

CALL  AGSETC  ('LABEL/NAME.',  'T') 

CALL  AGSETI  ( 'LINE/NUMBER. ' ,  100) 

CALL  AGSETC  (' LINE/TEXT. ’ ,  '  ') 

C 

CALL  AGSETC  ( ' LABEL/HAME . ' ,  'B') 

CALL  AGSETI  ( 'LINE/NUMBER. ' ,-100) 

CALL  AGSETC  ( ’ LIHE/TEXT . ' ,  '  *) 

C 

CALL  AGSETC  ('LABEL/NAME.',  'L') 

CALL  AGSETI  ( ' LINE/NUMBER . ' ,  100) 

CALL  AGSETC  ('LINE/TEXT.',  '  ') 

CALL  AGSETC  ('LABEL/NAME.',  'R') 

CALL  AGSETI  ( 'LINE/HURBER. ' ,-100 ) 

CALL  AGSETC  ('LINE/TEXT.',  '  ') 

C - turn  OFF  RIGHT  AXIS 

CALL  AGSETF  ( ' AXIS/RIGHT/CONTROL . ' , 0 . 0 ) 

C - SET  YMIN  AND  YHAX  FOR  FIRST  PLOT 

CALL  AGSETF  ( ' X/HININUM . ' ,  XIMIN) 

CALL  AGSETF  ( 'X/MAXINUH. ' ,  XIMAX) 

CALL  AGSETF  < ' Y/MINIMUN . ’ ,  0.0) 

CALL  AGSETF  t  '  Y/NAXlNUfl. '  ,  120.0) 

C - 

CALL  AGSTUP  ( XI , 1 , 1 , ICPl , 1 , VI , 1 , 1 , ICPl , 1 ) 

CALL  AGBACK 
CALL  GSPLCI  (7) 

CALL  DASHDC  (  •  $$$$$S$$SSS  "  "  "  "  "  $SSS$$$$$S ’ , 3 , 1 ) 
c - plot  curve  It  xl  VS  yi 

CALL  AGCURV  (XI , 1 , Yl , 1 , ICPl , -1 ) 

c 

c-  — — 

c - SET  UP  FOR  SECOND  PLOT 

C -  DEFINE  THE  TOP,  BOTTOH,  LEFT  AND  FIGHT  LABELS 

c - the  larger  line  number  writes  above  previous  labels 

CALL  AGSETC  ( ' LABEL/NAHE . ' ,  'T’) 

CALL  AGSETI  ( ’ LINE/NUHBER . ' ,  100) 

CALL  AGSETC  (' LINE/TEXT .', TOPLBL ) 

C 

CALL  AGSETC  ('LABEL/NAHE.',  'B') 

CALL  AGSETI  ( ' LIHE/HUHBER . ' , -1 1 0 ) 

CALL  AGSETC  (' LINE/TEXT. ' ,X1LABEL ) 

C 

CALL  AGSETC  ('LABEL/NAHE.',-  'L') 

CALL  AGSETI  ( ’ LTNE/NUHBER . ' ,  110) 

CALL  AGSETC  ( 'LINE/TEXT. ' ,T1LABEL) 

CALL  AGSETC  ('LABEL/NAHE.',  'B') 

CALL  AGSETI  ( ' LIHE/NUHBER. ' ,  -100) 

CALL  AGSETC  ( 'LINE/TEXT. ' ,Y2LABEL) 

C  TURN  OFF  THE  TOP,  BOTTOH  )WD  LEFT  AXES 
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f 


c— 

c— 

c— 

c 

c— 

c- 

c- 


CKLL  AGSETr  » 

■’“^LrAGEEt"  r  SlGHT/NUnERlC/TY PE  .  •  ,  1  •  E6  , 

OLL  AGSTUP  <X1,1.1.ICP1,1.Y2.1.1>JCP1,1I 
CftLL  AGBACE 

CALL  GSPLCI  .  .  1) 

call  dashoc 

—  plot  curve  2I  El  ys  Y2 

KORESTATION 


C 


END 
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